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.98.9% 382.7 / ‑2.9 88intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (46/41);  minor base T (1/0);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_pae_pgi/4234519‑4234808
                                                                                                                                                   |                                                                                                                                               
gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                                >  2:583127/1‑148 (MQ=255)
    atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCCTCGCCCCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑atcatc                                                                                                                                           >  2:946953/1‑149 (MQ=255)
     tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                          >  1:581933/1‑149 (MQ=255)
      ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                               <  1:455426/144‑3 (MQ=255)
      ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                               >  2:455426/1‑142 (MQ=255)
       taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGTCCACCCCGAGAATAGATGCATTGAAGTGGCGCTTGCATTG‑ATCAtcgtc                                                                                                                                        >  2:291919/1‑149 (MQ=255)
                gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAg                                                                                                                               >  1:631012/1‑149 (MQ=255)
                gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAg                                                                                                                               <  1:291568/149‑1 (MQ=255)
                   aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAAATGGCAGAATTTAGTGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGcc                                                                                                                            >  1:786948/1‑149 (MQ=255)
                    tgtaCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTATGCCg                                                                                                                           <  2:1163180/146‑1 (MQ=255)
                    tGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                           >  1:1256839/1‑149 (MQ=255)
                     gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCgg                                                                                                                          >  1:878168/1‑149 (MQ=255)
                      cACCCCAAGGGCGTTTTATCCTTTTCACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGa                                                                                                                         <  1:198193/149‑1 (MQ=255)
                       cccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                        <  2:1371450/148‑1 (MQ=255)
                       aCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                        >  1:1486351/1‑149 (MQ=255)
                       aCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                        >  1:494757/1‑149 (MQ=255)
                         cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGTGGCGGTTGCAGTG‑GTTATCGTCGTTCTGTTGGCCGGATaa                                                                                                                      >  2:531711/1‑149 (MQ=255)
                          ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  1:1543929/149‑1 (MQ=255)
                           cAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                    <  2:929196/149‑1 (MQ=255)
                               ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                >  2:299493/1‑149 (MQ=255)
                               ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATCGAAGAGGCGCTAGCAGTG‑ATCATCGTCGAGATGTAGGCAGGCGAAGGAGtt                                                                                                                >  2:1170575/1‑149 (MQ=255)
                                  gTTTTATTCTTTTTACGAAATTTGGGCATCACCTCTGTATAACGGGCGGAGATTAGAGGCGTAAATCCCCTGGCCACCCGGAGAATAGATGCATTGAAGAGGCGCATGCAGTG‑ATCATCGTCGCTATGTAGGCCGGATAAGGCGTTCAc                                                                                                             <  1:1228729/149‑1 (MQ=255)
                                         ccTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTCCCGTG‑ATCCTCGTCGATATGTAGGCCGGCTTAGGCGTCCCggccgcct                                                                                                      >  2:1368587/1‑147 (MQ=255)
                                           tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTCGGCCGGATAAGGCGTTTTCGCCGCCTcc                                                                                                    >  2:193105/1‑149 (MQ=255)
                                           tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                    >  2:464226/1‑149 (MQ=255)
                                            ttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                   <  1:464226/149‑1 (MQ=255)
                                                aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               >  1:596407/1‑149 (MQ=255)
                                                 cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGAGCTTGGCAT‑TATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              >  1:780014/1‑149 (MQ=38)
                                                  gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGGCGATATGTAGGCCGGAAAAGGCGGTCAAGCCGCATCCGGCAAcc                                                                                             >  2:539947/1‑149 (MQ=255)
                                                        tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAGGAGGAGCTTGCAGTG‑ATGATCGTCGATATGTAGAACGGATACGGCGTTCACGCCGCATCGGGCAACCGCTGcc                                                                                       >  2:587924/1‑149 (MQ=255)
                                                             aTCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                  >  2:257741/1‑149 (MQ=255)
                                                               cGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACCCGTTCCCTGATTCg                                                                                >  2:1477432/1‑149 (MQ=39)
                                                                 cccctGTGAAACTGGCAGAGATCAGGGGCGTAAATCCCCTGGCCACCCCGGGAATAGATGCATTGAAGGGGCGCTTGCTGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  2:213256/146‑1 (MQ=37)
                                                                 ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  2:631012/149‑1 (MQ=255)
                                                                   tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCTCCTGGCCACCCCGAGAATAGATGCATTGCAGAGGCGCTTCCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGAGTTCACGCTGCATCCGGCAACAGAATCCTGATGCGAcga                                                                            >  2:1586385/1‑148 (MQ=255)
                                                                   tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           <  1:434029/118‑1 (MQ=255)
                                                                   tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           >  2:434029/1‑118 (MQ=255)
                                                                       tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        >  1:440979/1‑149 (MQ=255)
                                                                        aaaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                                                                       >  1:272337/1‑149 (MQ=255)
                                                                          aaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                     <  1:1076804/149‑1 (MQ=255)
                                                                           aCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                    <  2:952622/149‑1 (MQ=255)
                                                                           aCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                    <  2:1486351/149‑1 (MQ=255)
                                                                            cTGGCAGAGATTAGAGGCGTAAACCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  1:400086/149‑1 (MQ=255)
                                                                             tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                  >  2:309026/1‑149 (MQ=255)
                                                                               gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCAACCCATACTTGTTACTAACCGTTCGCGTCttt                                                                >  2:1420822/1‑148 (MQ=37)
                                                                               gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:911959/1‑149 (MQ=255)
                                                                                cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               <  2:347044/149‑1 (MQ=255)
                                                                                  gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGACTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:1019763/1‑149 (MQ=255)
                                                                                   agaTTAGGGGCGTAAACCCCCTGGCCACCCGGGAAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  2:1549136/149‑1 (MQ=37)
                                                                                   agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  1:58951/149‑1 (MQ=255)
                                                                                        agagGCGTAAATCCTCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTa                                                      <  1:5031/149‑1 (MQ=255)
                                                                                         gatgCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcg                                                      <  1:1368587/146‑5 (MQ=37)
                                                                                          agGCGTAAATCCCCTGGCACCCCCGTTAAGAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACa                                                    >  1:1596340/1‑149 (MQ=25)
                                                                                           ggCGGAAATGCCCGGGCCACCTCGAGAATAGAGGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAg                                                   <  1:587924/149‑1 (MQ=25)
                                                                                            gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAgg                                                  >  1:509264/1‑149 (MQ=255)
                                                                                            gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAgg                                                  >  2:1199104/1‑149 (MQ=255)
                                                                                                aaaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGa                                              <  2:1497544/149‑1 (MQ=255)
                                                                                                    ccccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGcc                                          >  1:1423187/1‑149 (MQ=255)
                                                                                                    ccccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGcc                                          <  1:1468722/149‑1 (MQ=255)
                                                                                                     cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCg                                         >  1:1545408/1‑149 (MQ=255)
                                                                                                      ccTGGCGACCCCGAGAATAGATGCATTGAAGCGGCGTTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGa                                        <  1:946953/149‑1 (MQ=255)
                                                                                                      ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCAATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATCCGACGCGGGCGCATCTTAGCAACA‑CGACAGGTCGATGccag                                        >  2:1567286/1‑147 (MQ=37)
                                                                                                       cTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGat                                       <  2:581933/149‑1 (MQ=255)
                                                                                                        tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGata                                      <  1:309026/149‑1 (MQ=255)
                                                                                                               cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCAGGGTAAGGCGTTCAAGTCGCATCCGGCAACCGATGCCTGATTCGACGCCGTCTCGTCTTATCATTC‑CCACAGGTCGATGCCGATATGTACAt                               >  2:381641/1‑149 (MQ=255)
                                                                                                                ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATc                              >  2:1106899/1‑149 (MQ=255)
                                                                                                                ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATc                              <  1:1335419/149‑1 (MQ=255)
                                                                                                                ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATc                              <  2:440979/149‑1 (MQ=255)
                                                                                                                   agaATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTa                           <  2:770361/149‑1 (MQ=255)
                                                                                                                     aaTAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGAAAAGCACCTCGTact                         >  2:528587/1‑147 (MQ=25)
                                                                                                                      aTAGATGCATTGAAGAGGCTCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTc                        <  1:175341/149‑1 (MQ=255)
                                                                                                                       tAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCg                       <  1:1106899/149‑1 (MQ=255)
                                                                                                                       tAGATGCATTGAAGAGGCGCTTGCAGTG‑AGCGTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATCTACATCGTATTCg                       <  1:1321480/149‑1 (MQ=255)
                                                                                                                        aGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCgg                      <  1:1192819/149‑1 (MQ=255)
                                                                                                                        aGATGCATTGAAGAGGCGCTTGCAGTG‑AGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCgg                      <  2:290698/149‑1 (MQ=255)
                                                                                                                          aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                    >  2:773518/1‑149 (MQ=255)
                                                                                                                           tGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                   <  1:747230/149‑1 (MQ=255)
                                                                                                                            gCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCact                  >  2:409896/1‑147 (MQ=255)
                                                                                                                            gCAGTGAAGAGGCGCTTGCAGTG‑ATCATCGTCTATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                  <  1:1090579/149‑1 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCATTTa                >  1:655218/1‑149 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                >  2:752237/1‑149 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                >  1:1379298/1‑149 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CGACAGGTCGATGCCGATATGTACATCGGATTCGTCACTTa                >  2:1372253/1‑149 (MQ=37)
                                                                                                                              aTTGAAGAGGCACTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGGTATGTACATCGTATTCGGCAATTa                >  1:1013956/1‑149 (MQ=255)
                                                                                                                               ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTCAGCCGTGTTAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGGCGATGCCGATATGTACATCGTATTCGGCAATTaa               >  1:1589040/1‑149 (MQ=37)
                                                                                                                               ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa               <  2:1423187/149‑1 (MQ=255)
                                                                                                                                  aaGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc            <  2:1379298/149‑1 (MQ=255)
                                                                                                                                     agGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa         <  2:272337/149‑1 (MQ=255)
                                                                                                                                       gcgcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  2:647294/97‑1 (MQ=255)
                                                                                                                                       gcgcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            >  1:647294/1‑97 (MQ=255)
                                                                                                                                       gcgcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc       >  1:1475471/1‑149 (MQ=255)
                                                                                                                                            tGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGCTTAATCAGGCTCTATA‑CCAAACATCTATATGTACATCTTATTCGGCAATTAATACATAGCACgat  <  1:381641/149‑1 (MQ=21)
                                                                                                                                                   |                                                                                                                                               
GTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGC‑CTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_pae_pgi/4234519‑4234808

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