Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6830TA22.2% 394.8 / 48.7 188intergenic (+139/+234)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (35/98);  new base A (38/0);  total (88/100)
Fisher's exact test for biased strand distribution p-value = 5.18e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 190 total aligned reads displayed.

ATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGC  >  CP009273/2222533‑2222832
                                                                                                                                                       |                                                                                                                                                     
   atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                                                    <  1:1052399/151‑1 (MQ=255)
   atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCAGcc                                                                                                                                                    >  2:4454518/1‑150 (MQ=255)
   atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCa                                                                                                                                                    >  1:4051445/1‑151 (MQ=255)
     atTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGCCCCTACCCCccgccccccgcccc                                                                                                                                                  >  2:2395776/1‑151 (MQ=255)
     atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGCGGCCCAACGCCCCCCCCCCCAGCaac                                                                                                                                                  >  1:528710/1‑149 (MQ=255)
     atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                                                                                                                  <  1:4443792/151‑1 (MQ=255)
     atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                                                                                                                  <  1:2148486/151‑1 (MQ=255)
          aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                                                             <  1:3431213/151‑1 (MQ=255)
          aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                                                             <  1:1164726/151‑1 (MQ=255)
            aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                                                                                                                           <  1:3233687/151‑1 (MQ=255)
              ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctc                                                                                                                                         <  2:3988331/151‑1 (MQ=255)
               taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAAAGCCCCCCCCCCCCGCccccccccgcccg                                                                                                                                        >  1:5042248/1‑144 (MQ=255)
               taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                                                        <  1:554537/151‑1 (MQ=255)
                    aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAACGCCCCCCCCCCCCCCACCCCCCCGCGCTGGCgg                                                                                                                                   >  1:4166457/1‑150 (MQ=255)
                       gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACa                                                                                                                                <  1:2436990/151‑1 (MQ=255)
                         catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGCCCAAACCCCCCCCCCCCCCCCCCCCCCCGCACTGGCGAACaaa                                                                                                                              >  1:1050246/1‑151 (MQ=255)
                         catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGCCCAAACCCCCCCCCCCCAGCACCCCCCCGCACAGGGGAACaaa                                                                                                                              >  1:1000017/1‑151 (MQ=255)
                          atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGGCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGGCAATCGCCGCCGCCCCCGGCACCCCCGGGCGCAGGGGAACaaaa                                                                                                                             >  2:4735555/1‑151 (MQ=25)
                             aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                                                          <  1:4292467/151‑1 (MQ=255)
                               gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGCGGCCCAAAACCCCCCCCCCCCGCACCCCCCCGCGCAGGGGGACaaaaaaacc                                                                                                                        >  1:3240363/1‑146 (MQ=25)
                               gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAAAGCCGCCCCCCCCGGCCCCCCCCGGCACTGGCGAACaaaaaagcc                                                                                                                        >  2:1407673/1‑146 (MQ=25)
                                tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccg                                                                                                                       <  1:2750404/151‑1 (MQ=255)
                                 ccGAAAATAACGATAATATTTCAGTATATCAACCTGTTGCGCTCCTTTCTTGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccgc                                                                                                                      <  1:2845716/151‑1 (MQ=18)
                                 ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATAGCCGCCGCCCCCTGCAACCCCCGGCACTGGCGAAAAAAATagccgc                                                                                                                      >  2:1088909/1‑151 (MQ=25)
                                   gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                                                    <  2:3368794/151‑1 (MQ=255)
                                   gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAAAGCCCCCGCCCCCAGCACCCCCCGGCTATGGCGCAAAAAATagccgcag                                                                                                                    >  2:1123718/1‑149 (MQ=25)
                                     aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGGCAAACGCCGCCGCCCCCAGCCCCCCCCGGCTCTGGCGAACAAAAAAGCCgcggcg                                                                                                                  >  2:1498769/1‑151 (MQ=21)
                                     aaaTAACGAGAATAGTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCATTGCTTCAAATAAGAACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCCAACAAAAACGCCGCTgcg                                                                                                                  <  2:2125772/151‑1 (MQ=21)
                                      aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCCGGCTCATGCGAAAAAAAACGCCCCAgcgc                                                                                                                 >  2:1873870/1‑151 (MQ=21)
                                        tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCgg                                                                                                               <  2:4650719/151‑1 (MQ=255)
                                         aaCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCCCCCCCAGCACCCCCGGGCTCTGGCGAACAAAAAAGCCGCTGCGCGGt                                                                                                              >  2:3277474/1‑151 (MQ=21)
                                           cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGc                                                                                                            <  2:4083463/151‑1 (MQ=255)
                                            gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAAACCCCCCCCCCCCCGCCCCCCCCCGCGAAGGCGAAAAAAAAACCCCCAGCGCGGTggc                                                                                                           >  1:4908885/1‑149 (MQ=255)
                                            gagaAGATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCAGGCGAACAAAATCGCCCCGGCGCGGTGcc                                                                                                           >  2:3586958/1‑151 (MQ=255)
                                              caaTATTTCCGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCCTAGAAGCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACACAATCGCCGCTGCGCGGTGCCCt                                                                                                         <  2:3486209/150‑1 (MQ=18)
                                               aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAAAACCCCCCCCCCCCCGCACCCCCCCGCTCATGCGAACAAAAAAGCCCCGGCGCgggggccac                                                                                                        >  1:5382088/1‑144 (MQ=21)
                                               aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                                                        <  2:3960557/151‑1 (MQ=255)
                                               aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCGCCCCCAGCACCCCCCGGCTCATGCGAACAAAATCGCCCCTGCGCGGTGCCCTc                                                                                                        >  2:3199813/1‑151 (MQ=21)
                                                atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCCGCCCAACGCCCCCCCCCCCAACACCCCCCGGCTCGTGCGAACAAAAAAGCCCCAGCCCCGTGcccacc                                                                                                       >  1:2672099/1‑148 (MQ=18)
                                                  agttCAGTCTCTCATCCTGTTGCGCTCCTGGCATGTGCATTGCTTCATATAAGCACTGGCGCAAGGCGCGC‑GCTGGGGGCGGCCAATCCCGGCCGCCCCCTTGAGCCCCGGGCTGAGGGGGACCAAATGGCGGCTGCGCGGGGCGCGCggg                                                                                                     >  2:4956268/3‑150 (MQ=255)
                                                   tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCCCCCCCGCCCCCAGCACCCCCCGGCTAATGCGAAAAAAAACGCCCCAGCGCGGTGCCCACGGcg                                                                                                    >  2:1177834/1‑150 (MQ=18)
                                                    ttCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACCCCGCCGCCCCCTGCACCCCCGGGCTATGGCGAACAAAATAGCCCCAGCGCGGTGCCCTCGGCtt                                                                                                   >  2:4323278/1‑151 (MQ=18)
                                                      cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                                                 <  1:4121933/151‑1 (MQ=255)
                                                       aGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCCCCCCCGGCACCCCCGGGCACTGGCGAACAAAATAGCCGCGGCGCGGTGCCCACGGCTTATc                                                                                                >  2:3114279/1‑151 (MQ=18)
                                                        gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAATAGCCGCTGAGCGGTGCCCACGGCTTAtac                                                                                               >  2:364541/1‑149 (MQ=18)
                                                          ctctcATCCTGTTGCGCTCCTGTCATTTGCTGTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGTGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                                             <  2:4186093/151‑1 (MQ=255)
                                                             tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                                                          <  1:1968715/151‑1 (MQ=255)
                                                              cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                                                         <  1:1757221/151‑1 (MQ=255)
                                                                tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                                                       <  1:430830/151‑1 (MQ=255)
                                                                   tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAAAC‑CCCCCCCCCCCCGCCCCCCCCCGCCAAGGCGGACAAAAAAGCCCCAGCGCGGTGGCCCCCGCGGATACCCTACGGcaaa                                                                                    >  1:3364197/1‑148 (MQ=255)
                                                                   tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GTAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                                                    <  1:74326/151‑1 (MQ=255)
                                                                   tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                                                    <  2:3532371/151‑1 (MQ=255)
                                                                     ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                                                  <  2:706343/151‑1 (MQ=255)
                                                                        cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGCGGCCCAAACCCCCCCCCCCCACCACCCCCCGGCGCAGGCGGACAAAAAAGCCCCCGCGCGGGGCCCACCGCGGATACCCTACCGCAACCCGGt                                                                               >  1:1337184/1‑151 (MQ=255)
                                                                          cTCCGGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                                             <  1:3683496/151‑1 (MQ=255)
                                                                              tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                                                         <  1:250842/151‑1 (MQ=255)
                                                                               gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcgc                                                                        <  1:1204396/151‑1 (MQ=255)
                                                                               gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCGCCCCCACCACCCCCCGGCTCAGGCGAAAAAAAACACCCCAGCGCGGTGCCCAAGCCATATCACGTAAGGCAAACCGGTAGGgcgc                                                                        >  2:5364780/1‑151 (MQ=12)
                                                                                tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATAGCCGCCGCCCCCAGAACCCCCCGGCACAGGCCAACAAAAAAGCCCCAGCGCGGTGCCCACCGCTTAAACCATACGGCAACCCGGTCGGGcgcg                                                                       >  2:3499967/1‑151 (MQ=17)
                                                                                 cATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                                                      <  1:956413/151‑1 (MQ=255)
                                                                                   tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                                    <  1:3663085/151‑1 (MQ=255)
                                                                                       cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCCAACGCCCCCCCCCCCAGCACCCCCCGGCGCGGGCGAAAAAAAACCCCCCGGCGCGGGGGCCGCGGCGTAAACCAGAACGCGACCCGGGCGGGCCCGAGGGaaa                                                                >  2:2333786/1‑150 (MQ=255)
                                                                                         ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                              <  2:3296370/151‑1 (MQ=255)
                                                                                         ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                              <  2:4145639/151‑1 (MQ=255)
                                                                                         ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCCCCCCCCCCAGCAACCCCCGGCGCTGGCGAACAAAAAACCCCCACCCCGGGGCCCCCGCCGTAAACCCTAACGCGACCGCGTCGGGCGCGGCGTAACaa                                                              >  2:2922233/1‑150 (MQ=9)
                                                                                         ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCCGCCAAAAGCCCCCCCCCCCATCACCCCCCGGCTCTGGCCAACAAAAAAGCCGCATCGCGGTGCCCACGGCTTAAACCATACCGCTACCCGGTCGGGCGCGAGGTAACaa                                                              >  2:1223359/1‑150 (MQ=9)
                                                                                           gCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCCATCGCCCCCCCCCCCAGCCCCCCCCGGCGCAGGCGAAAAAAAAACCCCCGGCGCGGTGGCCACGGCGTAAACCATACGGCGACCCGGGCGGGGGCGAGGGAACaaac                                                            >  2:930381/1‑148 (MQ=255)
                                                                                           gCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                                            <  1:2697645/151‑1 (MQ=40)
                                                                                                   taATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                                                    <  1:4726538/151‑1 (MQ=34)
                                                                                                    tgccACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATAGCCCGCGCCCCCGGCACCCCCCGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCGCGGCGTATACCTTAACGCGACCGGGTCGGGCGCGAGGGAACATACCATTaaaaa                                                   >  2:2972609/4‑150 (MQ=1)
                                                                                                       cACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                                                <  1:161148/151‑1 (MQ=31)
                                                                                                       cACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                                                <  1:160346/151‑1 (MQ=31)
                                                                                                        aCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                               <  2:1368041/151‑1 (MQ=31)
                                                                                                          tcgcgcAAGGACCGC‑GCCGGGGGCGCCCAATCGCCGCCGCCCCCCGCACCCCCGGGCACAGCCCAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAACCCCCACGGCAACCGGGTCGGGCGCGAGGTAACAACCCTGTAAAACGCGccc                                             <  2:3786449/149‑1 (MQ=11)
                                                                                                           tgcgcAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCCGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACTGCTTCCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                                            <  2:2110267/150‑1 (MQ=9)
                                                                                                              gcAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg                                         <  1:3819449/151‑1 (MQ=31)
                                                                                                                 aGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc                                      <  1:2367438/151‑1 (MQ=31)
                                                                                                                 aGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc                                      <  1:4687813/151‑1 (MQ=31)
                                                                                                                 aGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc                                      <  1:4687770/151‑1 (MQ=31)
                                                                                                                    aGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                                   <  1:2068496/151‑1 (MQ=31)
                                                                                                                       gc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc                                >  1:596846/1‑151 (MQ=31)
                                                                                                                       gc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGAACCACCGCGCTCTGGCCATCATAATAGGCGCTCCGCTGGGCCCGCGGCTTATCGCTTACGGCAACAGGGGCGGGCGCGAGGTAACATCCCGGAAACACGCGCCCTCCCCCCACATcc                                >  2:4046705/1‑151 (MQ=2)
                                                                                                                          gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt                              <  1:427488/151‑1 (MQ=31)
                                                                                                                            aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                            >  2:5136042/1‑151 (MQ=31)
                                                                                                                             gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                           <  2:1877636/151‑1 (MQ=31)
                                                                                                                                  cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg                      <  2:5041013/151‑1 (MQ=31)
                                                                                                                                      cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                  >  1:2523900/1‑151 (MQ=31)
                                                                                                                                        aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg                <  1:2285202/151‑1 (MQ=31)
                                                                                                                                        aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg                >  2:1508419/1‑151 (MQ=31)
                                                                                                                                         tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGc               >  1:3464782/1‑151 (MQ=31)
                                                                                                                                           gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCGTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAACCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCt             <  1:632850/151‑1 (MQ=14)
                                                                                                                                            ccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCtt            <  2:2493198/151‑1 (MQ=31)
                                                                                                                                            ccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCtt            <  2:2229741/151‑1 (MQ=31)
                                                                                                                                             cgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTc           >  1:351766/1‑151 (MQ=31)
                                                                                                                                                cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg        >  1:1192799/1‑151 (MQ=31)
                                                                                                                                                  cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa      >  2:5177639/1‑151 (MQ=31)
                                                                                                                                                    cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc    >  2:4851411/1‑151 (MQ=31)
                                                                                                                                                    cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc    >  1:2824957/1‑151 (MQ=31)
                                                                                                                                                    cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAAc    <  2:859869/151‑1 (MQ=31)
                                                                                                                                                     ccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACg   <  1:4261899/151‑1 (MQ=31)
                                                                                                                                                       |                                                                                                                                                     
ATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGC  >  CP009273/2222533‑2222832

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