Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,819,892 C→A 94.2% 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,8920CA94.2% 552.4 / 24.6 191intergenic (‑10/+320)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (9/2);  new base A (98/82);  total (107/84)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 216 total aligned reads displayed.

GACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAA  >  CP009273/1819747‑1820036
                                                                                                                                                 |                                                                                                                                                
 aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtactt                                                                                                                                               <  1:216859/148‑1 (MQ=255)
 aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtactt                                                                                                                                               >  2:508161/1‑148 (MQ=255)
 aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtactt                                                                                                                                               <  2:3321772/148‑1 (MQ=255)
  ccGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttc                                                                                                                                              <  1:1636930/148‑1 (MQ=255)
   cGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttcc                                                                                                                                             <  2:1933861/148‑1 (MQ=255)
   cGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttcc                                                                                                                                             <  2:1475030/148‑1 (MQ=255)
    gAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttcct                                                                                                                                            <  1:3097457/148‑1 (MQ=255)
    gAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttcat                                                                                                                                            <  2:2449296/148‑3 (MQ=255)
     aaCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATacttcctt                                                                                                                                           >  1:1082967/1‑148 (MQ=255)
     aaCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATacttcctt                                                                                                                                           <  1:3184823/148‑1 (MQ=255)
     aaCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATacttcctt                                                                                                                                           >  2:2227438/1‑148 (MQ=255)
      aCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCttt                                                                                                                                          >  1:2723654/1‑148 (MQ=255)
      aCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCttt                                                                                                                                          <  1:199250/148‑1 (MQ=255)
       cTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTc                                                                                                                                         <  1:776771/148‑1 (MQ=255)
       cTTGCCTTTGTGGTGCATCATCGGATTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTc                                                                                                                                         >  2:2530111/1‑148 (MQ=255)
        ttGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCa                                                                                                                                        >  1:2487469/1‑148 (MQ=255)
         tGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAg                                                                                                                                       <  2:2917424/148‑1 (MQ=255)
              ttGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttatt                                                                                                                                  >  2:1239349/1‑148 (MQ=255)
                gtggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttattta                                                                                                                                >  2:2129776/1‑148 (MQ=255)
                 tggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttatttat                                                                                                                               <  2:2757753/148‑1 (MQ=255)
                  ggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTtatttatt                                                                                                                              >  2:790416/1‑148 (MQ=255)
                      catcatCGGCTTCGCGTCAGCCGGTGCTGCGGTTGTGGTGTCTGCGGCATGGTGCATGTTAGCCGCGCCAAGAGACATGGTAGAGGCATCAAACAGTGCAGTTAATTTACGCATTTTCTTTATACTTCCTTTCAGTTATTTATTActg                                                                                                                          <  2:1772563/148‑3 (MQ=255)
                        tcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCt                                                                                                                        <  1:3371141/148‑1 (MQ=255)
                              gCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTa                                                                                                                  <  2:2397952/148‑1 (MQ=255)
                               cTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTaa                                                                                                                 >  2:301999/1‑148 (MQ=255)
                               cTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTaa                                                                                                                 >  2:1164361/1‑148 (MQ=255)
                                 tCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGt                                                                                                               <  2:2097663/148‑1 (MQ=255)
                                           cGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgt                                                                                                     <  2:2283585/148‑1 (MQ=255)
                                           cGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgt                                                                                                     <  2:2429913/148‑1 (MQ=255)
                                             gTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgtgt                                                                                                   <  2:3237737/148‑1 (MQ=255)
                                               gctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTg                                                                                                 <  1:1837080/148‑1 (MQ=255)
                                                  gcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACg                                                                                              <  1:275338/148‑1 (MQ=255)
                                                     gTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACgaga                                                                                           <  1:1339171/148‑1 (MQ=255)
                                                         tggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTaa                                                                                       >  1:223037/1‑148 (MQ=255)
                                                           gtgtCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACt                                                                                     <  2:1960311/148‑1 (MQ=255)
                                                            tgtCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACtt                                                                                    >  2:487488/1‑148 (MQ=255)
                                                             gtCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACttt                                                                                   >  1:968134/1‑148 (MQ=255)
                                                               cTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTAc                                                                                 >  2:1302036/1‑148 (MQ=255)
                                                                tGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACt                                                                                >  2:2152830/1‑148 (MQ=255)
                                                                tGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACt                                                                                >  2:1121743/1‑148 (MQ=255)
                                                                 gCGACATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTg                                                                               >  2:426363/1‑148 (MQ=255)
                                                                    gCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACt                                                                            <  1:1112294/148‑1 (MQ=255)
                                                                    gCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACt                                                                            >  1:727150/1‑148 (MQ=255)
                                                                     cATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACtt                                                                           <  1:1426548/148‑1 (MQ=255)
                                                                      aTGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACttt                                                                          <  1:1050931/148‑1 (MQ=255)
                                                                              ggTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCaa                                                                  >  1:350099/1‑148 (MQ=255)
                                                                               gTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAAt                                                                 >  1:3194740/1‑148 (MQ=255)
                                                                               gTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAAt                                                                 >  1:1974382/1‑148 (MQ=255)
                                                                                   gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                             >  2:404710/1‑148 (MQ=255)
                                                                                   gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                             >  2:405767/1‑148 (MQ=255)
                                                                                    ccGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                            <  1:1471679/148‑1 (MQ=255)
                                                                                    ccGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                            >  1:2935747/1‑148 (MQ=255)
                                                                                    ccGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATc                                                            >  1:1002805/1‑148 (MQ=255)
                                                                                         ccAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGc                                                       >  1:2269899/1‑148 (MQ=255)
                                                                                          cAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCa                                                      <  1:391456/148‑1 (MQ=255)
                                                                                           aaGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCaa                                                     >  1:2618645/1‑148 (MQ=255)
                                                                                           aaGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCaa                                                     <  2:2682496/148‑1 (MQ=255)
                                                                                             gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                   >  2:2478245/1‑148 (MQ=255)
                                                                                             gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                   >  2:2650963/1‑148 (MQ=255)
                                                                                             gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                   >  2:3438748/1‑148 (MQ=255)
                                                                                              agCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACgg                                                  >  2:3117753/1‑148 (MQ=255)
                                                                                                      tAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACaa                                          <  2:171635/148‑1 (MQ=255)
                                                                                                      tAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACaa                                          <  2:1509612/148‑1 (MQ=255)
                                                                                                      tAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACaa                                          >  1:3440673/1‑148 (MQ=255)
                                                                                                       agagGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAAt                                         >  1:1128295/1‑148 (MQ=255)
                                                                                                          ggCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGaa                                      <  1:1491759/148‑1 (MQ=255)
                                                                                                             aacaaacaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgt                                   <  2:1737417/148‑1 (MQ=255)
                                                                                                             aacaaacaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAgt                                   <  2:3348985/148‑1 (MQ=255)
                                                                                                                aaacaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTa                                <  2:2240653/148‑1 (MQ=255)
                                                                                                                 aacaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaa                               >  1:50679/1‑148 (MQ=255)
                                                                                                                 aacaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaa                               >  2:871203/1‑148 (MQ=255)
                                                                                                                  acaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaa                              >  1:2336027/1‑148 (MQ=255)
                                                                                                                  acaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaa                              >  1:1203498/1‑148 (MQ=255)
                                                                                                                  acaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaa                              >  1:1155065/1‑148 (MQ=255)
                                                                                                                  acaGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaa                              >  1:1761014/1‑148 (MQ=255)
                                                                                                                    aGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCNTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaa                            >  1:1784576/1‑148 (MQ=255)
                                                                                                                      tGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAac                          >  1:440172/1‑148 (MQ=255)
                                                                                                                       gCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAaca                         <  1:3342189/148‑1 (MQ=255)
                                                                                                                       gCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAaca                         <  1:458651/148‑1 (MQ=255)
                                                                                                                        cAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAacac                        >  1:1591415/1‑148 (MQ=255)
                                                                                                                          gTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtt                      <  2:3227884/148‑1 (MQ=255)
                                                                                                                          gTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtt                      <  2:654887/148‑1 (MQ=255)
                                                                                                                            tAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtttt                    >  1:2670861/1‑148 (MQ=255)
                                                                                                                            tAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtttt                    >  1:1826594/1‑148 (MQ=255)
                                                                                                                            tAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtttt                    >  1:611839/1‑148 (MQ=255)
                                                                                                                               tttACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgcg                 >  2:1807480/1‑148 (MQ=255)
                                                                                                                                ttACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgctc                >  2:2331219/1‑146 (MQ=255)
                                                                                                                                ttACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgcgc                >  2:371889/1‑148 (MQ=255)
                                                                                                                                ttACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgcgc                >  2:57619/1‑148 (MQ=255)
                                                                                                                                 tACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGcc               <  1:2427357/148‑1 (MQ=255)
                                                                                                                                      atatTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTa          <  2:3217022/148‑1 (MQ=255)
                                                                                                                                      atatTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTa          >  1:3058298/1‑148 (MQ=255)
                                                                                                                                         ttCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATgg       <  1:2112886/148‑1 (MQ=255)
                                                                                                                                         ttCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATgg       <  1:1121743/148‑1 (MQ=255)
                                                                                                                                          tCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga      <  2:17043/148‑1 (MQ=255)
                                                                                                                                          tCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga      <  2:16988/148‑1 (MQ=255)
                                                                                                                                          tCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga      >  2:1551007/1‑148 (MQ=255)
                                                                                                                                           cTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGgag     >  2:595499/1‑148 (MQ=255)
                                                                                                                                            tatataCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGGGTGCCGTGTTGAATAGATTAACTTTACTGAATTTAGCGTAAATTAATCAGAGCAAAGGTAAAAAAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGgaga    <  1:2530111/148‑1 (MQ=255)
                                                                                                                                            tatataCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGgaga    >  1:3029660/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAA  >  CP009273/1819747‑1820036

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: