Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,114 A→T 5.5% intergenic (‑48/+547) 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,1140AT 5.5% 1372.4 / 7.2 465intergenic (‑48/+547)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (182/243);  new base T (7/31);  total (189/276)
Fisher's exact test for biased strand distribution p-value = 3.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 474 total aligned reads displayed.

AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGA  >  CP009273/359965‑360255
                                                                                                                                                     |                                                                                                                                             
ggccgttgcttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCa                                                                                                                                               >  1:4641211‑M1/140‑150 (MQ=255)
  ccgttgcttcgcaacgttcaaatccgctcccggcggatttgtcctacacaggagagcgttcaccgacaaacaacagataaaaccaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGt                                                                                                                                             >  1:3689504‑M1/138‑150 (MQ=255)
    gttgcttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagatacaacgacaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGttt                                                                                                                                           >  2:1731266‑M1/136‑150 (MQ=255)
    aGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGttt                                                                                                                                           >  2:603716/1‑150 (MQ=255)
     ttgcttcgcaacgttcacctccgctcccggcggatttgtcctccgcaggagagcgttcaccgaccaacaccagataaaacgcaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTa                                                                                                                                          >  2:2612131‑M1/135‑150 (MQ=255)
     gTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTa                                                                                                                                          >  2:4601031/1‑150 (MQ=255)
     gTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTa                                                                                                                                          >  2:1878375/1‑150 (MQ=255)
     gTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTa                                                                                                                                          <  2:2443510/150‑1 (MQ=255)
         ttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggc                                                                                                                                      <  1:1401996‑M1/20‑1 (MQ=255)
         aatcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggc                                                                                                                                      >  2:465049/1‑150 (MQ=255)
          tcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggcg                                                                                                                                     >  2:4490599‑M1/130‑150 (MQ=255)
          atcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcg                                                                                                                                     <  1:1967404/150‑1 (MQ=255)
          atcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcg                                                                                                                                     >  1:2193102/1‑150 (MQ=255)
           tcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcgg                                                                                                                                    >  2:4088868/1‑150 (MQ=255)
            caaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcggg                                                                                                                                   <  1:707898/150‑1 (MQ=255)
             caacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggcgggc                                                                                                                                  >  2:3851310‑M1/127‑150 (MQ=255)
             aaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcgggc                                                                                                                                  <  1:4073923/150‑1 (MQ=255)
              aacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggcgggcg                                                                                                                                 <  1:4781063‑M1/25‑1 (MQ=255)
               cTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                >  2:1671684/1‑150 (MQ=255)
               cTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                >  2:1324646/1‑150 (MQ=255)
                tGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTc                                                                                                                               <  1:562382/150‑1 (MQ=255)
                cgttcaaatccgctcccggcggatttgtcctactccggagagcgttcaccgacaaacaacagataaaacgacaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTc                                                                                                                               >  1:2864862‑M1/124‑150 (MQ=255)
                   ttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTg                                                                                                                            >  1:797649/1‑150 (MQ=255)
                   ttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGATTATGGCGGGCGTCCTg                                                                                                                            >  1:243995/1‑150 (MQ=255)
                   tcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTg                                                                                                                            >  2:2193102‑M1/121‑150 (MQ=255)
                    caaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcggtttgtttgatgccgggggtattAATTGCTGGGAggttatggcgggcggccggc                                                                                                                           <  1:1251347‑M1/31‑7 (MQ=255)
                     aaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaccgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGcc                                                                                                                          >  2:4433216‑M1/119‑150 (MQ=255)
                     aCCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGcc                                                                                                                          >  1:3000465/1‑150 (MQ=255)
                      ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                         <  1:3426175/150‑1 (MQ=255)
                         tgtgGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGcc                                                                                                                      >  1:3940616/1‑150 (MQ=255)
                          gtgGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCa                                                                                                                     <  1:1511190/150‑1 (MQ=255)
                            ggAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                   >  2:2820705/1‑150 (MQ=255)
                             tcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                                                                                  <  1:3469635‑M1/40‑1 (MQ=255)
                              cccggcggatttgtcctactcaggagagcgttcaccgacaaacaccagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                 >  1:4700608‑M1/110‑150 (MQ=255)
                                cGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGATTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                               <  2:3746793/150‑1 (MQ=255)
                                 ggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCAGCCACCCTCCg                                                                                                              <  2:1068435‑M1/44‑1 (MQ=255)
                                 gACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                              <  2:2889897/150‑1 (MQ=255)
                                    aTCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                           >  1:3172766/1‑150 (MQ=255)
                                     tCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                          >  2:1374977/1‑150 (MQ=255)
                                          gtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTtgg                                                                                                     <  1:1979464‑M1/53‑2 (MQ=255)
                                          agGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                     >  2:3928707/1‑150 (MQ=255)
                                             cACTTCACCGATTGCCTGGGGCTTACCATCCAGCGCCACCACCCAGTGCTGGAGCTCCGTTGATAAGCTGACAAACATGAGAATTCCGATCCCCAATTCCGGGCTGTTTAGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                  <  1:3918782/150‑1 (MQ=37)
                                              aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                 >  2:104043/1‑150 (MQ=255)
                                                ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                               >  2:3218760/1‑150 (MQ=255)
                                                 tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCGCCGGGCCGTTGCTTCGCaa                                                                                              <  2:2881124‑M1/60‑1 (MQ=255)
                                                    ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                           >  1:4528452/1‑150 (MQ=255)
                                                    ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                           >  2:1452767/1‑150 (MQ=255)
                                                      agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                         <  2:198381‑M1/65‑1 (MQ=255)
                                                       cgtGCCCGCGGCTTCCCAGCCCGCGCCACCAGCCCGTGCTGGCGCTCCGATGATACGCTGACAAACATGAGAATTCCGACCCCCACTTCCTGGCTGTTTTGGGCGGGCGGCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                        <  1:1064795/148‑1 (MQ=21)
                                                       cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                        <  2:2732882/150‑1 (MQ=255)
                                                            cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                   >  2:2302672/1‑150 (MQ=255)
                                                             aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCg                                                                                  >  2:4009583/1‑150 (MQ=255)
                                                              caccgacaaacaacagataaaacgaaagtcccagtctttcgactgagcatttcgttttagttgatgcgtggggcattAATTAATGGAAGTTTTTGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                 <  1:1093475‑M1/73‑1 (MQ=255)
                                                               cggataaccaccCCGCCCCACCATCCCGTGCTGGAGCTCCTTTTATACGCTTACAAACATGAGAATTCCGATCCCAATTTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                <  1:1017592/144‑1 (MQ=16)
                                                                  cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                             <  2:4286002/150‑1 (MQ=255)
                                                                   ttACCATCCCGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                            <  1:1150123/150‑1 (MQ=255)
                                                                   ttACCACCCCGCGCCACAATCCCGTGCAGGTCCTCCGATGATAAGCTGTCAAACATGAGAATGCCGATCCCCAATTCCTGGAAGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                            <  2:1280747/150‑1 (MQ=21)
                                                                   acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCCACCGTCAAAACCGCgcccc                                                                            <  1:4156618‑M1/78‑2 (MQ=255)
                                                                     aaacaacagataaaacgaaaggcccagtctttcgactgagtctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                          <  2:1913330‑M1/80‑1 (MQ=255)
                                                                     aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTa                                                                                                                                          <  1:649942‑M1/16‑1 (MQ=255)
                                                                       cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACGTTCAAATACGCAGCccgcgg                                                                        >  1:889113/1‑150 (MQ=37)
                                                                        caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggcgggcg                                                                                                                                 <  1:4220241‑M1/25‑1 (MQ=255)
                                                                               gcaaccaTCCAGTGCTGGAGCTCCTTTGATAAGCTTAAAAACATGATAATTCCGATCCACAATTCCTGGCTTTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                <  1:322433/147‑1 (MQ=18)
                                                                                accccgaatggcccagtctttcgactgcgcctttcgttttatttgaggcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                               >  2:2415326‑M1/60‑150 (MQ=255)
                                                                                  accaTCCCGCCCAGGAGCTCCGCTGAGAAGCTGTAAAACATTCGAATTCGGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                             <  2:251992/150‑1 (MQ=17)
                                                                                   ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                            <  2:1133345/150‑1 (MQ=255)
                                                                                       ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                        >  2:2976981/1‑150 (MQ=255)
                                                                                         aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                      <  2:1576497/150‑1 (MQ=255)
                                                                                              tgtccgtcgacggagcctttcgttttatttgctgcccggggcaatAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                 >  2:4380229‑M1/46‑150 (MQ=255)
                                                                                                   cTCCGATGATCAGCTGTCAAACATGAGAAGTCCGAGCCCCATTTCCTGGCAGTTGATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGGTCAcc                                            <  2:3380974/150‑1 (MQ=17)
                                                                                                     cgacggcgcattacgttttttttgatgccgggggaattACTCCCTGGCAGTTTATGGCGGGCGTCCCCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGa                                          >  1:2058021‑M1/39‑150 (MQ=255)
                                                                                                        atgatAATCTGTCAAAAATGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                       <  2:4470655/150‑1 (MQ=255)
                                                                                                         tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                      <  1:378169/150‑1 (MQ=255)
                                                                                                         tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                      <  2:3309842/150‑1 (MQ=255)
                                                                                                         tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAACTTCAAATACGCACCCCGCGGATTTGTCCTACACACGAGAGcgttcaccaaaaaa                                      <  2:475382‑M1/116‑7 (MQ=255)
                                                                                                          gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTGAAAATCCGCTCCCGGCGGATTGGTCCTACTCAGGAGAGCGTTCACCGacaaac                                     <  2:1513505‑M1/117‑1 (MQ=255)
                                                                                                          gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCCTCGCAAAGATAAAAAACGCACCCCCCGGATTTGTACTACTCAGGGGAGCGTTCACCCacaaac                                     <  2:51460‑M1/117‑1 (MQ=255)
                                                                                                           agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCACATCCGCTCCCGGCGGATTTGTGCTACGCAGGAGAGCGGTCACCGacacaca                                    <  1:3545149‑M1/118‑1 (MQ=255)
                                                                                                               ggtGTCAAACATGATAATGCCGATCCCCAATTCGGGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                <  1:2568861/148‑1 (MQ=14)
                                                                                                                   tttttatttgatgcgtgggtaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                            >  2:3689678‑M1/25‑150 (MQ=255)
                                                                                                                     tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                          <  2:3544033‑M1/128‑1 (MQ=255)
                                                                                                                      aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                         <  2:3430413/150‑1 (MQ=255)
                                                                                                                        cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGTTGGCAAAGTTAAAATACGCTCCCGGCGGATTTTTACTAAACAGGATAGCGTTCAACGACAAAAAACAGATtaaaaaa                       >  1:4320285/1‑147 (MQ=17)
                                                                                                                        cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATCCGCGCCCGGCCGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGa                       >  1:1918705/1‑150 (MQ=17)
                                                                                                                        cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAAACCGCTCCCCGCGGAATTGTCCGACTCAGGAGAGCGGTCCCCCACCCACAACCGATAAAAccc                       >  2:421248/1‑148 (MQ=17)
                                                                                                                         aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAAACCGCGCCCGGCGGAGTGGTGCGACTCAGGAGAGCGGTCACCGACAAACAACAGCTaaaaacaa                      >  1:4379294/1‑146 (MQ=14)
                                                                                                                            agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                   <  2:727502/150‑1 (MQ=255)
                                                                                                                               tttCCGATCCCCATTTCCGTGCTGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                <  1:1748445/149‑1 (MQ=11)
                                                                                                                                ttCCGATCCCCAATTCCTGTATTTTTTGGGGGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa               <  1:1671684/150‑1 (MQ=11)
                                                                                                                                cctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                          >  1:2335457‑M1/12‑139 (MQ=255)
                                                                                                                                   cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc            <  2:1809671/150‑1 (MQ=39)
                                                                                                                                    gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATACGCACCCGGCGGCTTTTTCCTACTAAGGAGAGCGTTCACCGACAAAAAACAGATAAAAAGAAAGGCCCAGTCt           <  2:463998‑M1/143‑1 (MQ=255)
                                                                                                                                      gaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTCCCAACGTTCAAATACCCTCCCCGCGGATGTGTACTACACAGGGGGGCGGTCACCGACCAACCACAAATAAAAACAAAAGCCCACTCttt         >  1:3055974/6‑150 (MQ=2)
                                                                                                                                       ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc        <  2:1156874/150‑1 (MQ=35)
                                                                                                                                        cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGTTCAAAACCCCTCCCGGCGGCTTTGTCCTACACAGGAGAGCGGTCACCGCCACAAAACACATAAAAAGAAATGCCCACTCTTTCg       >  1:2685095/1‑150 (MQ=11)
                                                                                                                                        attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg       <  2:3165148‑M1/147‑1 (MQ=255)
                                                                                                                                         ccAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa      <  2:895076/150‑1 (MQ=17)
                                                                                                                                          cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc     <  2:655548/150‑1 (MQ=33)
                                                                                                                                           aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAACAGATAAAACGAAATGCCCACTCTTTCGACt    >  2:2078956/1‑150 (MQ=255)
                                                                                                                                            agtCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg   <  2:3573983/148‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                                             
AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGA  >  CP009273/359965‑360255

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: