Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.99.4% 766.7 / ‑0.2 164intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (73/88);  total (75/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 172 total aligned reads displayed.

GGGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_pae_pgi/4234517‑4234807
                                                                                                                                                     |                                                                                                                                             
gggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                               <  1:1257607/150‑1 (MQ=255)
 ggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                              <  1:1707119/150‑2 (MQ=255)
  gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                             <  2:1104602/150‑3 (MQ=255)
  gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                             <  3:789273/150‑3 (MQ=255)
  gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                             <  4:818630/150‑3 (MQ=255)
  gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                             <  3:789289/150‑3 (MQ=255)
         taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtcg                                                                                                                                     >  1:929602/1‑150 (MQ=255)
         taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtcg                                                                                                                                     <  4:866960/150‑1 (MQ=255)
          aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACCGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                    >  2:206124/1‑150 (MQ=255)
           cGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCAGCCCAAACAAAAAATAGACACATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                   <  4:55552/150‑1 (MQ=255)
            gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGata                                                                                                                                  >  3:1160971/1‑150 (MQ=255)
              aGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATg                                                                                                                                <  4:1137578/150‑1 (MQ=255)
              aGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATg                                                                                                                                <  3:699160/150‑1 (MQ=255)
                 agagATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAg                                                                                                                             <  3:656755/150‑1 (MQ=255)
                    gaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGcc                                                                                                                          <  3:618359/150‑1 (MQ=255)
                       gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCg                                                                                                                         <  1:1923180/148‑1 (MQ=255)
                        cACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAt                                                                                                                      >  3:996765/1‑150 (MQ=255)
                         aCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATa                                                                                                                     >  4:1099947/1‑150 (MQ=255)
                          ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                    <  3:16671/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                   <  1:1802191/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTATACAAGATGTGGGCATCGCCTTTGTAGAACTGGCAGGGATTAGAGGCGTAAATGCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                   <  3:424680/150‑1 (MQ=255)
                             cAAGGGCGTTTTATCCTTTTTACGAGCTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                 >  3:271718/1‑150 (MQ=255)
                             cAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                 <  3:192326/150‑1 (MQ=255)
                              aaGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                <  3:133773/150‑1 (MQ=255)
                                 ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                             >  4:315189/1‑150 (MQ=255)
                                  gCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                            <  2:1817798/150‑1 (MQ=255)
                                      tttGTCCTTTTTACGATATGGGGTCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        <  2:1966672/150‑1 (MQ=255)
                                      tttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        <  4:19656/150‑1 (MQ=255)
                                            cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTa                                                                                                                              >  1:554639/1‑122 (MQ=255)
                                            cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTa                                                                                                                              <  2:554639/122‑1 (MQ=255)
                                            cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                  >  1:1423022/1‑150 (MQ=255)
                                             tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                 >  1:401201/1‑150 (MQ=255)
                                             tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                 >  3:597804/1‑150 (MQ=255)
                                                 tACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                             <  1:272515/150‑1 (MQ=255)
                                                 tACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                             >  1:1764321/1‑150 (MQ=255)
                                                  aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        >  2:1269284/1‑138 (MQ=255)
                                                   cttgaTGTGGTCATCGCCTCTGTATAACTGGCAGTTATTTGAGGCGTAAATCCCCTGGCCCCCCCTCGAATAGATTCATTTTAGAGGCGCTTGCCGTCATCATCGGCGATATGTTGGCCGGCTAAGTCGTTCTCGCCGCATCCGGCAAcc                                                                                           >  4:537673/4‑150 (MQ=255)
                                                     agaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                         <  1:1605411/150‑1 (MQ=255)
                                                            ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCg                                                                                                                         <  2:450876/111‑1 (MQ=255)
                                                            ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                  >  4:418585/1‑150 (MQ=255)
                                                            ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGGGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCTCAACCGGCAACCGATGCCTGa                                                                                  >  2:788385/1‑150 (MQ=255)
                                                             gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGCGAATAGATGCATTCAAGAGCCGCTTGCAGTGCTCATCGTCGTTATGTATGCCGTATAAGGCGTTCTCGCCGCGTCCGGCAACCGATGCCTGAt                                                                                 <  4:1051495/150‑1 (MQ=255)
                                                                 cGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                             <  2:1379301/150‑1 (MQ=255)
                                                                  gCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                            >  1:1324715/1‑150 (MQ=255)
                                                                    ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                  >  3:481483/1‑110 (MQ=255)
                                                                    ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                  <  4:481483/110‑1 (MQ=255)
                                                                    ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        <  2:1479659/120‑1 (MQ=255)
                                                                    ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        >  1:1479659/1‑120 (MQ=255)
                                                                      ctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATCCCTTTAATAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGAcgcga                                                                        <  1:1021139/150‑2 (MQ=38)
                                                                       tGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                      >  4:574353/1‑135 (MQ=255)
                                                                         tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                                                                     >  1:2120/1‑150 (MQ=255)
                                                                         tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGCCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                                                                     >  4:646563/1‑150 (MQ=255)
                                                                            aaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGAcgcgaccgcgt                                                                  <  3:148193/150‑8 (MQ=39)
                                                                                 gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                    >  1:961349/1‑107 (MQ=255)
                                                                                 gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                    <  2:961349/111‑5 (MQ=255)
                                                                                  cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                            >  3:873183/1‑150 (MQ=255)
                                                                                  cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                            >  2:1647417/1‑150 (MQ=255)
                                                                                     agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                             <  2:319134/98‑1 (MQ=255)
                                                                                     agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                             >  1:319134/1‑98 (MQ=255)
                                                                                     agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                         <  2:929602/150‑1 (MQ=255)
                                                                                     agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                         >  4:800864/1‑150 (MQ=255)
                                                                                      gaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                        <  3:1127763/150‑1 (MQ=255)
                                                                                       aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                       <  1:116163/150‑1 (MQ=255)
                                                                                       aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATAAGACGCGACCGCGTCGCATCAGGca                                                       <  1:865793/150‑2 (MQ=37)
                                                                                         tAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatc                                                     <  3:799946/150‑4 (MQ=38)
                                                                                           gagGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                   <  3:800864/150‑1 (MQ=255)
                                                                                              gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                <  3:1092853/150‑1 (MQ=255)
                                                                                                    aTCCCCTGGCCCCCCCGAGAATAGATGCATCGAAGAGGCGCGGGCAGCGATCACCCGCCGTATGCCGGCCGGATAACCCGTTCACGCCCCATCCGGCAACCGCTGCCTGACAGGGCGCCGCCGCTGCTTCTCGGGCCTCCAGGGCGAtct                                          >  4:77591/1‑148 (MQ=255)
                                                                                                     tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                         >  1:1587185/1‑150 (MQ=255)
                                                                                                      ccccTGGCCACCCCGGGAATAGATGCATTGGAGAGGCGCTTGCAGTGAACATCGTCGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                        <  2:1605954/150‑1 (MQ=255)
                                                                                                       cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                       >  2:547552/1‑150 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                      >  1:987376/1‑150 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                  <  3:536144/150‑1 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                  <  2:1587185/150‑1 (MQ=255)
                                                                                                             ccACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTa                                 <  4:265683/150‑1 (MQ=255)
                                                                                                                ccccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                              <  3:315189/150‑1 (MQ=255)
                                                                                                                 cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                             <  2:1900623/150‑1 (MQ=255)
                                                                                                                  ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                      <  2:697821/140‑1 (MQ=255)
                                                                                                                   cGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                           >  4:222940/1‑150 (MQ=255)
                                                                                                                        aTAGAGGCATTGAAGAGGCACTTGCAGTGATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                      <  4:732106/150‑1 (MQ=255)
                                                                                                                          aGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                   >  1:1961534/1‑55 (MQ=255)
                                                                                                                          aGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                   <  2:1961534/55‑1 (MQ=255)
                                                                                                                          aGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                    <  2:1324715/150‑1 (MQ=255)
                                                                                                                           gATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                   <  1:547552/150‑1 (MQ=255)
                                                                                                                           gATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                   <  1:168336/150‑1 (MQ=255)
                                                                                                                            aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATAAGACGCGACCGCGTCGCATCAGGcat                                                      >  2:1225759/1‑112 (MQ=21)
                                                                                                                            aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATAAGACGCGACCGCGTCGCATCAGGcat                                                      <  1:1225759/114‑3 (MQ=21)
                                                                                                                              gCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTCCATCGTATTCGGCAAtt                >  2:465143/1‑150 (MQ=255)
                                                                                                                                 ttGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt             <  1:52140/150‑1 (MQ=255)
                                                                                                                                  tGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa            <  2:1764321/150‑1 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        <  2:1164368/56‑1 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                   <  1:1040873/61‑1 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                   >  2:1040873/1‑61 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                    >  2:1565704/1‑124 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                    <  1:1565704/124‑1 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACa          >  4:976618/1‑150 (MQ=255)
                                                                                                                                      gagGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa        <  1:299424/150‑1 (MQ=255)
                                                                                                                                       agGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAgg                                                         <  1:1197826/100‑1 (MQ=21)
                                                                                                                                         gcgcTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa     >  2:522095/1‑150 (MQ=255)
                                                                                                                                            cTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAACCGATGCCTGAGGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACga  >  1:1276223/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GGGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_pae_pgi/4234517‑4234807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: