Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,819,892 C→A 6.6% intergenic (‑10/+320) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,819,8920CA 6.6% 678.6 / 16.6 181intergenic (‑10/+320)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (93/76);  new base A (6/6);  total (99/82)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 209 total aligned reads displayed.

CGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAA  >  CP009273/1819745‑1820039
                                                                                                                                                   |                                                                                                                                                   
cGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATGtc                                                                                                                                                     <  1:2433199/148‑1 (MQ=255)
cGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtc                                                                                                                                                     <  1:222977/148‑1 (MQ=255)
cGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtc                                                                                                                                                     <  2:411323/148‑1 (MQ=255)
cGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtc                                                                                                                                                     <  1:2987748/148‑1 (MQ=255)
 ggACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtcc                                                                                                                                                    <  2:68635/148‑1 (MQ=255)
  gACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtcct                                                                                                                                                   <  2:362152/148‑1 (MQ=255)
   aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtcctt                                                                                                                                                  >  2:1854558/1‑148 (MQ=255)
    ccGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtccttc                                                                                                                                                 >  1:1694795/1‑148 (MQ=255)
    ccGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttc                                                                                                                                                 >  1:3279176/1‑148 (MQ=255)
      gAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtccttcct                                                                                                                                               <  1:287348/148‑1 (MQ=255)
        aCTTGCCTTTGTGGTGCATCATGGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCttt                                                                                                                                             <  1:2007661/148‑1 (MQ=255)
        aCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCttt                                                                                                                                             <  1:496215/148‑1 (MQ=255)
        aCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCttt                                                                                                                                             >  1:2855090/1‑148 (MQ=255)
              cTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGtta                                                                                                                                       <  2:477231/148‑1 (MQ=255)
               tttGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttat                                                                                                                                      <  1:765381/148‑1 (MQ=255)
                  gtggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttattta                                                                                                                                   <  2:3480980/148‑1 (MQ=255)
                  gtggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttattta                                                                                                                                   >  2:1608058/1‑148 (MQ=255)
                   tggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttatttat                                                                                                                                  <  2:772066/148‑1 (MQ=255)
                    ggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTtatttatt                                                                                                                                 <  2:1495347/148‑1 (MQ=255)
                        catcatCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACgg                                                                                                                             <  2:3519563/148‑1 (MQ=255)
                        catcatCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACgg                                                                                                                             <  2:1347157/148‑1 (MQ=255)
                         atcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGc                                                                                                                            <  2:628663/148‑1 (MQ=255)
                         atcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGc                                                                                                                            >  1:2576781/1‑148 (MQ=255)
                            atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCttt                                                                                                                         >  1:313311/1‑148 (MQ=255)
                            atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCttt                                                                                                                         >  1:313431/1‑148 (MQ=255)
                            atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCttt                                                                                                                         <  2:1846577/148‑1 (MQ=255)
                             tcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTc                                                                                                                        >  1:1930294/1‑148 (MQ=255)
                                 cTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTaa                                                                                                                    >  2:2489884/1‑148 (MQ=255)
                                       gTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCg                                                                                                              <  1:1678771/148‑1 (MQ=255)
                                        tCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGt                                                                                                             >  1:3484062/1‑148 (MQ=255)
                                        tCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGt                                                                                                             <  1:1308946/148‑1 (MQ=255)
                                         cAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTg                                                                                                            <  2:3538571/148‑1 (MQ=255)
                                            ccGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCg                                                                                                         <  1:594985/148‑1 (MQ=255)
                                              ggTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgtg                                                                                                       <  2:131124/148‑1 (MQ=255)
                                                  ctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGa                                                                                                   <  2:107610/148‑1 (MQ=255)
                                                    gcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACg                                                                                                 <  1:2533116/148‑1 (MQ=255)
                                                     cGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACga                                                                                                <  1:24753/148‑1 (MQ=255)
                                                       gTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACgaga                                                                                              <  2:1201617/148‑1 (MQ=255)
                                                        tAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                             >  1:2401952/1‑148 (MQ=255)
                                                        tAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                             <  2:2258188/148‑1 (MQ=255)
                                                        tAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                             <  1:691075/148‑1 (MQ=255)
                                                        tAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAt                                                                                             >  1:1796108/1‑148 (MQ=255)
                                                         aGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGAtt                                                                                            >  1:2645759/1‑148 (MQ=255)
                                                           tggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTaa                                                                                          <  1:3464601/148‑1 (MQ=255)
                                                             gtgtCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACt                                                                                        <  1:1273370/148‑1 (MQ=255)
                                                                         tGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTa                                                                            >  1:2215393/1‑148 (MQ=255)
                                                                                ggTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCaa                                                                     >  2:3212640/1‑148 (MQ=255)
                                                                                   tAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTa                                                                  <  2:155309/148‑1 (MQ=255)
                                                                                     gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTACTTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                                <  2:1825433/148‑1 (MQ=255)
                                                                                     gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                                <  1:857595/148‑1 (MQ=255)
                                                                                      ccGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                               >  1:784639/1‑148 (MQ=255)
                                                                                      ccGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                               <  2:749162/148‑1 (MQ=255)
                                                                                       cgcgCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCa                                                              >  2:943461/1‑148 (MQ=255)
                                                                                        gcgcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCag                                                             <  1:1608058/148‑1 (MQ=255)
                                                                                        gcgcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCag                                                             <  2:1654545/148‑1 (MQ=255)
                                                                                        gcgcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTAAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCag                                                             <  2:1717891/148‑1 (MQ=255)
                                                                                            cAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCa                                                         >  2:1973550/1‑148 (MQ=255)
                                                                                            cAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCa                                                         >  1:2810958/1‑148 (MQ=255)
                                                                                               gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                      >  2:1353330/1‑148 (MQ=255)
                                                                                                agCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACgg                                                     <  2:2706758/148‑1 (MQ=255)
                                                                                                 gCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGt                                                    >  2:2691182/1‑148 (MQ=255)
                                                                                                 gCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGt                                                    <  1:2139420/148‑1 (MQ=255)
                                                                                                 gCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGt                                                    <  1:1617693/148‑1 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                                >  1:3164164/1‑148 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                                >  2:583802/1‑148 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                                >  2:239025/1‑148 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                                <  2:2297923/148‑1 (MQ=255)
                                                                                                      ggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGATATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAAc                                               >  2:2773425/1‑148 (MQ=255)
                                                                                                       gtcgagGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                              <  2:3181283/145‑1 (MQ=255)
                                                                                                       gTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                              <  1:1195337/148‑1 (MQ=255)
                                                                                                       gTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                              <  2:2944834/148‑1 (MQ=255)
                                                                                                       gTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                              <  1:2416426/148‑1 (MQ=255)
                                                                                                             gCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGaaa                                        >  1:2989951/1‑148 (MQ=255)
                                                                                                             gCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGaaa                                        >  1:3080706/1‑148 (MQ=255)
                                                                                                              caacaaACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAg                                       >  2:2732257/1‑148 (MQ=255)
                                                                                                              caacaaACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAg                                       <  2:1243985/148‑1 (MQ=255)
                                                                                                               aacaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgt                                      <  1:1020216/148‑1 (MQ=255)
                                                                                                                acaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgtg                                     >  2:1940429/1‑148 (MQ=255)
                                                                                                                acaaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgtg                                     >  2:1786756/1‑148 (MQ=255)
                                                                                                                  aaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTa                                   >  2:904874/1‑148 (MQ=255)
                                                                                                                      aGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaa                               <  1:1209091/148‑1 (MQ=255)
                                                                                                                       gTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaaa                              >  2:2178138/1‑148 (MQ=255)
                                                                                                                        tGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAac                             <  2:2120414/148‑1 (MQ=255)
                                                                                                                               aaTTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTACAACAATGAAAGTGTAAAAAACACTTTTg                      >  2:441494/1‑148 (MQ=255)
                                                                                                                               aaTTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTg                      <  1:1407598/148‑1 (MQ=255)
                                                                                                                                     cGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTACAACAATGAAAGTGTAAAAAACACTTTTGCGCCac                >  1:3275340/1‑147 (MQ=255)
                                                                                                                                        atatTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTa             <  1:1484629/148‑1 (MQ=255)
                                                                                                                                         tatTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGTGCCAATTAt            >  2:2948709/1‑148 (MQ=255)
                                                                                                                                          atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg           >  1:410168/1‑148 (MQ=255)
                                                                                                                                          atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg           >  2:1246944/1‑148 (MQ=255)
                                                                                                                                          atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg           >  1:389377/1‑148 (MQ=255)
                                                                                                                                           ttCTATATCCTTCCTTTCCGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATgg          <  2:2393194/148‑1 (MQ=255)
                                                                                                                                           ttCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATgg          <  2:2392195/148‑1 (MQ=255)
                                                                                                                                            tCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga         >  2:576070/1‑148 (MQ=255)
                                                                                                                                            tCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga         >  2:1007281/1‑148 (MQ=255)
                                                                                                                                            tCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga         >  1:2076087/1‑148 (MQ=255)
                                                                                                                                                tatCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaa     >  1:1569307/1‑148 (MQ=255)
                                                                                                                                                 atCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaaa    <  2:784639/148‑1 (MQ=255)
                                                                                                                                                 atCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaaa    >  1:2514295/1‑148 (MQ=255)
                                                                                                                                                  tccttcctTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaaaaa   >  2:38873/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAAAA  >  CP009273/1819745‑1820039

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