Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,172 T→A 7.6% intergenic (‑106/+489) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1720TA 7.6% 604.1 / 8.1 230intergenic (‑106/+489)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (51/151);  new base A (12/9);  total (63/167)
Fisher's exact test for biased strand distribution p-value = 4.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360023‑360321
                                                                                                                                                     |                                                                                                                                                     
gCCAGGGGCTTACCATCCCGCGCCACCATCCCGGGCTGGAGCTCCGTGGATACGCTGTCAAACATGAGAATGCCGATCCCCAAGTCCTGGATTTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                                       <  1:2370998/150‑1 (MQ=21)
  cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTGGCGTCGCAACCATCAAATcc                                                                                                                                                     >  1:4216966/1‑150 (MQ=38)
  cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                     <  2:1697394/150‑1 (MQ=255)
  cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                     >  1:3586821/1‑150 (MQ=255)
    caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                   <  2:2979274‑M1/73‑1 (MQ=255)
     cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                  <  2:3071305/150‑1 (MQ=255)
       cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGTCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCAGTTGCTTAGCAACgttcaaataacctac                                                                                                                                                <  2:3642857‑M1/76‑8 (MQ=255)
       cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATCCGCTcc                                                                                                                                                <  2:1453350‑M1/76‑1 (MQ=255)
        cTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                               <  2:3889060/150‑1 (MQ=255)
         ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCCACGTTCAAAAACGCTcccc                                                                                                                                              >  1:1304646/1‑149 (MQ=37)
         acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCCAATCCGCACCcg                                                                                                                                              <  2:2555649‑M1/78‑1 (MQ=255)
          caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCAACGTTCAaatacgcagcccg                                                                                                                                             <  1:1181588‑M1/79‑7 (MQ=255)
          caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTGGCCACCTTCAAAACCGCTCCcgg                                                                                                                                             <  2:67226‑M1/79‑1 (MQ=255)
            aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGGAACGTTCAAAGCCGGTCgcggcg                                                                                                                                           <  2:4265703‑M1/81‑1 (MQ=255)
              caacagataaaacgaaaggcccagtctttcgactgagcctgtcgttttagttgttgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGGGGCACCGGGCAAATCCGCGCCCGGCGGa                                                                                                                                         <  2:3032259‑M1/83‑1 (MQ=255)
              aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                         <  2:4134265/150‑1 (MQ=255)
               aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                        <  2:1766945‑M1/84‑1 (MQ=255)
                     ccccccaTCCCGGGCCGGAGCTCCGATTATAAGCTGACAAACATTATAATTCCGACCCCAATTTCCTGGCTGTTTCTGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                  <  1:3424279/146‑1 (MQ=17)
                       caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCGTCAAATACGCGCCCCGCGGATTTGtacaa                                                                                                                                >  1:2796548/1‑146 (MQ=25)
                        accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               <  1:4297579/150‑1 (MQ=255)
                           aTCCTGTGCTGGAGCTCCGATGATATGCTGACAAACATGAGAATTCCGACCCCCAATTCCTGGCAGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                            <  1:2318901/150‑1 (MQ=18)
                               ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                                    <  1:2499876‑M1/72‑1 (MQ=255)
                               ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                                    >  2:2499876‑M1/51‑122 (MQ=255)
                               ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                        <  2:4021666‑M1/100‑1 (MQ=255)
                               cgTGCCGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                        <  2:3611716/149‑1 (MQ=255)
                                gcccagtcttgcgactgcgactttggttttatttgttgccgggggaattAATTCCTGGATGTTTATGGCGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                       >  1:2716082‑M1/50‑150 (MQ=255)
                                 cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCAATGCCGCTCCCGGCGGATTGGTCCTACTCAGGagag                                                                                                                      <  2:2624265‑M1/102‑1 (MQ=255)
                                   cAGGAGCTCGTATGATACTCTGTCATAAATGCTACGACCGTTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                    <  2:599122/150‑1 (MQ=14)
                                    agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                   <  2:4260219‑M1/105‑1 (MQ=255)
                                     ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                  <  2:4162673/150‑1 (MQ=255)
                                            cGATGATAAGCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                           <  2:2228374/150‑1 (MQ=255)
                                             gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGGAACCGTTAAATACCCCACCGGCGGATTTGTGCTAAGCCGGAGTGCGGTAACCGaca                                                                                                          >  1:433083/1‑150 (MQ=18)
                                               tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACTTGAAATACGCAACCCGCGGATTTGTGCAAATCAAGGGAGCGGTGACCCacaaa                                                                                                        >  2:370002/1‑150 (MQ=18)
                                               tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCAACCGTCAAAACAGCGCCCCGCGGGTTTGTACGACACAGGCGAGCGGTGAccccccaa                                                                                                        >  2:370052/1‑144 (MQ=18)
                                               tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAATCCGCGCCCGGCGGCTTTGTCCGACGCAGGAGCGCGTTCACCGacaaa                                                                                                        >  2:804647/1‑150 (MQ=18)
                                                gatACGCTTTAAAACATGATAAGTCCGATCCCAAATTCCTGGCTGTTTAGGGGGGGCGCCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                       <  1:457494/150‑1 (MQ=14)
                                                gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCAAAACCGCGCCCCGCGGAATTGTCCAACTCAGGAGAGCGTTCACCCacaaac                                                                                                       <  1:2862227‑M1/117‑1 (MQ=255)
                                                  tAAGCTTTAAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                     <  1:4138830/150‑1 (MQ=255)
                                                  tAAGCTCACAAACATGAGAATACCGACCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                                     <  1:1951313/150‑1 (MQ=255)
                                                   cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCAAaacaac                                                                                                    <  2:1570209‑M1/120‑1 (MQ=255)
                                                   aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACCTTCAAATACGCGCCCCGCGGATTTGTGCAACAAAGGAGAGCGGTAAACCacaaacaaa                                                                                                    >  1:3521524/1‑149 (MQ=17)
                                                     gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                  <  2:92621/150‑1 (MQ=255)
                                                     gCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                  <  1:1674896/150‑1 (MQ=255)
                                                      cTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                 <  2:2164747/150‑1 (MQ=255)
                                                         gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                              <  2:2472886‑M1/126‑1 (MQ=255)
                                                         aaaaaaATCAGAATTCCTACCCCAAATACCTGGAATTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                              <  1:3059258/148‑1 (MQ=11)
                                                          cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGGTGAAAAACCCTCCCCCCCGGGTTGTAGGACGAAGGAGTGCGTTGACCCACAAAAAACAaaaaaa                                                                                             >  1:906227/1‑150 (MQ=12)
                                                            aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                           <  2:3485365/150‑1 (MQ=255)
                                                            aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTAAAAAACGCGACCCGCGGAGTTTTGCAACTAAGGAGAGCGTTAAACGACCAACAACACAAAAAAc                                                                                           >  1:2253553/1‑150 (MQ=12)
                                                            aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                           <  2:3822479/150‑1 (MQ=255)
                                                             tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                          >  1:4260219‑M1/21‑150 (MQ=255)
                                                             aCATGAGAATACCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                          <  1:3846884/150‑1 (MQ=255)
                                                              cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACACATAAAACGa                                                                                         >  2:1523453/1‑150 (MQ=255)
                                                              atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTATTGTCCTACTCAGGGGAGCGTTGACCGACACACAACAGTTAAAACGa                                                                                         <  2:3736905‑M1/131‑1 (MQ=255)
                                                               tttgttgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGGTCAAATACGCTCCCGGCGGAGTTTGCCTACTCAGGTGAGCGGTCACCGACAAACAACAGATAAAacgta                                                                                        <  1:1434770‑M1/132‑3 (MQ=255)
                                                                 tgttgcctggggaattATTTCGTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                      >  2:3657854‑M1/17‑150 (MQ=255)
                                                                    tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGGTGAAATACCCTCCCCGCGGCTTTGTGCTAATCAGGAGAGCGTTCACCCACAAACAAAAGATAAAAACAAAGGcc                                                                                   <  2:35423‑M1/137‑1 (MQ=255)
                                                                      ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                 <  2:1959328/150‑1 (MQ=255)
                                                                       tCCGATCCCCATTTCCTGGCTTTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                <  1:4165862/150‑1 (MQ=11)
                                                                       ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                <  2:311023‑M1/140‑1 (MQ=255)
                                                                        tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCAAATCCGCACCCGGCGGTTTTGTCCCTCTCCGGAGAGCGTTCACCCCCAAACAAAAGATAAATCCAATGGCCCAGt                                                                               <  1:969356‑M1/141‑1 (MQ=255)
                                                                        tggggaattAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTAATCAGGAGAGGGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGt                                                                               <  1:2463217‑M1/141‑1 (MQ=255)
                                                                         ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                               >  1:1051097‑M1/9‑133 (MQ=255)
                                                                         ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                               <  2:1051458‑M1/125‑1 (MQ=255)
                                                                          ggtaattAATTCGTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTTTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             >  1:3337770‑M1/8‑150 (MQ=255)
                                                                          gggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             >  1:3978897‑M1/8‑150 (MQ=255)
                                                                          gATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             <  2:1105129/150‑1 (MQ=255)
                                                                           aTCCCCATTTCGTGGCTGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            <  1:2666175/150‑1 (MQ=11)
                                                                           aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            <  2:2512589/150‑1 (MQ=38)
                                                                             ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          <  1:628891/150‑1 (MQ=35)
                                                                             ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          <  1:3563201/150‑1 (MQ=35)
                                                                             ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          <  2:2715504/150‑1 (MQ=35)
                                                                              cccAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCACGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                         <  1:1784193/150‑1 (MQ=25)
                                                                              attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACACCAGCTAAAACGCAAGGCCCAGTCTTTCg                                                                         <  2:4199050‑M1/147‑1 (MQ=255)
                                                                                tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAATCCGCGCCCCGCGGCTTTGGCCTACACAGGAGAGCGTTCACCCACCAACAACAGATAAAACGAAAGGCCCAGTCATTCGAc                                                                       <  2:418177‑M1/149‑1 (MQ=255)
                                                                                cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                       <  1:2882341/150‑1 (MQ=33)
                                                                                 atttCCTGGCTTTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                      <  1:843478/148‑1 (MQ=9)
                                                                                 aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                      <  2:1819434/150‑1 (MQ=32)
                                                                                  aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                     >  2:2412019/1‑150 (MQ=32)
                                                                                  aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                     >  2:532174/1‑150 (MQ=32)
                                                                                                             tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                          <  1:2418664/150‑1 (MQ=11)
                                                                                                              gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                         <  2:2932360/150‑1 (MQ=31)
                                                                                                                cccccACCCTCCGGGCCGGTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                       <  1:2718110/147‑1 (MQ=12)
                                                                                                                ccGCCACCCTCCGGGCCGTAGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                       <  2:2909998/150‑1 (MQ=17)
                                                                                                                 cGCCACCCGCCGGGCCGTTGCATCGCAACGTACAACTCCGCGCCCGGCGGATTGGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  2:1903363/150‑1 (MQ=11)
                                                                                                                  gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGACTGGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  2:2094991/150‑1 (MQ=11)
                                                                                                                  gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  1:204530/150‑1 (MQ=34)
                                                                                                                  gCCACCCGCCGGGCCGTTGCTTCGCAACGTTCACAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGGCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  2:4216966/150‑1 (MQ=11)
                                                                                                                      cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                                 <  1:2979274/150‑1 (MQ=38)
                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                    >  1:3932847/1‑134 (MQ=35)
                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                    >  1:3339450/1‑134 (MQ=35)
                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                    <  2:3932847/134‑1 (MQ=35)
                                                                                                                                   tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                    <  2:3339450/134‑1 (MQ=35)
                                                                                                                                      ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGTAGAc                 >  1:292809/1‑150 (MQ=255)
                                                                                                                                       tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACg                <  2:1635946/150‑1 (MQ=255)
                                                                                                                                        cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGGCGTTGTAGTTGATGCCTGGGGCCGTGGGAATTACTTCTTGAAGACGa               >  2:1393012/1‑150 (MQ=14)
                                                                                                                                         gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa              >  1:4094779/1‑150 (MQ=255)
                                                                                                                                                  aaaTCCGCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACCGAGAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg     <  2:3199101/150‑1 (MQ=255)
                                                                                                                                                   aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa               <  1:136131/139‑1 (MQ=255)
                                                                                                                                                     tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa  >  1:1769602/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360023‑360321

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: