Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,114 A→T 6.4% 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 6.4% 1137.6 / 10.2 385intergenic (‑48/+547)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (159/191);  new base T (5/29);  total (164/221)
Fisher's exact test for biased strand distribution p-value = 4.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAG  >  CP009273/359965‑360256
                                                                                                                                                     |                                                                                                                                              
   cgttgcttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgacaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGtt                                                                                                                                             >  2:3012042‑M1/137‑150 (MQ=255)
   cAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGtt                                                                                                                                             >  2:1552241/1‑150 (MQ=255)
    gttgcttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGttt                                                                                                                                            >  1:3351608‑M1/136‑150 (MQ=255)
    aGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGttt                                                                                                                                            <  2:1164312/150‑1 (MQ=255)
       tcaatcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATg                                                                                                                                         <  2:130442/150‑1 (MQ=255)
       gcttcgcaacgttcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagccgttcgttttatttgatggctggggtattAATTACTGGCAGTTTATg                                                                                                                                         <  1:432781‑M1/18‑1 (MQ=255)
        caatcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATgg                                                                                                                                        >  1:3721728/1‑150 (MQ=255)
          tcgcaacggtcaaatccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagatacaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATggcg                                                                                                                                      >  2:4482159‑M1/130‑150 (MQ=255)
          atcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcg                                                                                                                                      <  2:1477592/150‑1 (MQ=255)
              aCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcgggcg                                                                                                                                  >  1:1233568/1‑150 (MQ=255)
                 gTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTcc                                                                                                                               >  2:2890643/1‑150 (MQ=255)
                  tttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCt                                                                                                                              <  2:4550117/150‑1 (MQ=255)
                      actccgctcccggcggatttgtcctactccggagagcggtcaccgacacccaccagctcacacgaccggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                          >  1:2927017‑M1/118‑150 (MQ=255)
                       cTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCg                                                                                                                         >  2:542312/1‑150 (MQ=255)
                       atccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTc                                                                                                                                >  1:2806757‑M1/117‑143 (MQ=255)
                          cgctcccggcggatttgtcctcctcaggagaccgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcatggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                                                                                      >  1:4875986‑M1/114‑150 (MQ=255)
                          cgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                          >  2:4065202‑M1/114‑146 (MQ=255)
                           gctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCt                                                                                                                              <  2:4621608‑M1/29‑1 (MQ=255)
                            ggAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                    >  2:2164209/1‑150 (MQ=255)
                             gAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                                                                                   >  1:2688285/1‑150 (MQ=255)
                              cccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                  >  1:1128408‑M1/110‑150 (MQ=255)
                               cgggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtatttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                                                                                 >  1:3466447‑M1/109‑150 (MQ=255)
                               ccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                                                                                 >  2:2930995‑M1/109‑150 (MQ=255)
                                 gACATCCAGAGGCACTTCACCGCTTGCAAGCGGCTTACAATCCAGCGCCACAATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                               >  1:4778547/1‑150 (MQ=255)
                                   cggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATg                                                                                                                                         <  1:3551175‑M1/18‑1 (MQ=255)
                                   cATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                                             >  2:2436062/1‑150 (MQ=255)
                                    ggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCGCCGGGc                                                                                                            <  1:4221109‑M1/47‑1 (MQ=255)
                                     tCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                           >  2:3852315/1‑150 (MQ=255)
                                      ccAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTTAAACATGAGAATTACGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                                                                                          >  2:1792843/1‑150 (MQ=255)
                                      ccAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                                                                                          >  1:3047389/1‑150 (MQ=255)
                                        ttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                        <  2:2294506‑M1/51‑1 (MQ=255)
                                        agagGCACTTCACCGCTTGCCAGCGGCTTACCATCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                        <  2:3242198/150‑1 (MQ=255)
                                         tgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTg                                                                                                       <  1:3141248‑M1/52‑1 (MQ=255)
                                           ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                     <  2:3607038/150‑1 (MQ=255)
                                           ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                     >  1:810683/1‑150 (MQ=255)
                                              aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCg                                                                                                  >  2:278836/1‑150 (MQ=255)
                                                 tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                               <  2:1127441‑M1/60‑1 (MQ=255)
                                                  cACCGCTTGCCAGCGGCTTACCAGCCCCCCCCACCATCCCGTGCTGGAGCTCCGATTATACGCTGTCAAACATGATAATTCCGATCCCCAATTCCGGGCTTTTTAGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                              <  1:3938162/150‑1 (MQ=25)
                                                   aCCGCTTGCCAGCGGATTACCACCCCGCGCCACAAGCCAGTGCAGGAGCTCCTATTATATGCTGAAAAACATGATATTTCCGATCCCAAATTCCTGGCTTTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                             <  1:2722271/150‑1 (MQ=21)
                                                     cGCTTGCCAGCGGATTACCATCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAACCATGAGAATACCGATCCCCAATTCGTGGCAGTTTATGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                           <  1:1105340/150‑1 (MQ=37)
                                                      agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGATCGCAACgtta                                                                                          <  2:3249770‑M1/65‑2 (MQ=255)
                                                       gagcgttcaccgacaaccaacataaaaaactaaaggcccagtctttcgcctgagcctttcgttttattttatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                         >  2:1425771‑M1/85‑150 (MQ=255)
                                                           gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                     <  2:2768109‑M1/70‑1 (MQ=255)
                                                               accgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                 >  2:3915599‑M1/77‑150 (MQ=255)
                                                                ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                           >  2:2807305‑M1/76‑123 (MQ=255)
                                                                 cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCGCCGGGCCGGTGGGTCGCCACCATCAAATACGCTcc                                                                               <  1:2589827‑M1/76‑1 (MQ=255)
                                                                    caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                            <  1:3243439‑M1/79‑1 (MQ=255)
                                                                           cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                     <  2:3195400/150‑1 (MQ=255)
                                                                           cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                     <  1:2425648/150‑1 (MQ=255)
                                                                           cAGCGCCACCATCCAGAGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                     <  2:4726415/150‑1 (MQ=255)
                                                                             gcgcCACCATCCAGTCCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                   <  2:3365072/150‑1 (MQ=255)
                                                                                 caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCAAGTTAAAATCCGCGCCCGGCGGGTTTGTccaa                                                               >  1:1424823/1‑148 (MQ=25)
                                                                                  accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATAAGCTACCCGCGGATTTGtactac                                                              >  2:591431/1‑150 (MQ=25)
                                                                                   ccaTCCTGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                             <  1:4751483/150‑1 (MQ=255)
                                                                                   ccaTCCAGCGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                             <  2:4263878/150‑1 (MQ=255)
                                                                                     aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCTATCCCCAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                           <  1:4413347/150‑1 (MQ=255)
                                                                                      tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGt                                                                                                                                 <  1:1975431/79‑1 (MQ=255)
                                                                                      tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCGTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                          <  2:16174/150‑1 (MQ=255)
                                                                                      aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                          <  1:1937616‑M1/33‑1 (MQ=255)
                                                                                         ggcgccgtctttcgactgagcctttggttttattgtctgcgcggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                       >  2:3721728‑M1/51‑150 (MQ=255)
                                                                                           tGCTGGTGCGCCGTTGATAGGCTTTCAAAAATTATATTCCCGAACCCCAATTCGTGGCTGTTTATGGCGGGGGGCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                     <  1:4144299/150‑1 (MQ=14)
                                                                                           cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                     >  2:1945431‑M1/49‑150 (MQ=255)
                                                                                            ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAACACCCCACCCCGCGGCTTTGTGCTACACAGGAGAGc                                                    <  1:1672364‑M1/103‑1 (MQ=255)
                                                                                                gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                <  2:4311231/150‑1 (MQ=255)
                                                                                                  tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                              <  2:805187‑M1/109‑1 (MQ=255)
                                                                                                  tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                              <  2:1884422‑M1/109‑1 (MQ=255)
                                                                                                  gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                              <  1:1559056/150‑1 (MQ=255)
                                                                                                   ctcgattgcgcttgtcgatttatttgaggcctgggtaattAATTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                             >  1:2977297‑M1/41‑150 (MQ=255)
                                                                                                       gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCGCGGCGGATTTGTCCTACTCAGGAGAGCGTGCACCGaca                                         >  2:4875986/1‑150 (MQ=255)
                                                                                                       gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                         >  2:2694115/1‑150 (MQ=255)
                                                                                                           agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                     <  2:4353218‑M1/118‑1 (MQ=255)
                                                                                                            gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                    <  2:2594431‑M1/119‑1 (MQ=255)
                                                                                                              agcttaaaaaCATTCGATTTCCTATCCCCATTTCCTGGCTGTTTAGGGGGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACGAacaaca                                  <  1:915556/143‑1 (MQ=11)
                                                                                                                 tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                               <  1:1466350/150‑1 (MQ=255)
                                                                                                                 tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGGGAAAATCCGCACCCGGCGGATTTGTCCGACAAAGGCGAGCGTTGACCGACAAAAAAAAGAt                               >  1:753565/1‑150 (MQ=17)
                                                                                                                  cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCAGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATCCGCTACCGGCGGATTTGTCAAACTAAGGAGAGCGTTCACCGACAAACAACAGATa                              <  2:3502475‑M1/125‑1 (MQ=255)
                                                                                                                   gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                             >  2:2728098‑M1/25‑150 (MQ=255)
                                                                                                                   gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                             <  2:4777758‑M1/126‑1 (MQ=255)
                                                                                                                     aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                           <  2:4241469/150‑1 (MQ=255)
                                                                                                                      aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                          <  2:3061135/150‑1 (MQ=255)
                                                                                                                       ccATGTGAATCCCGACCCCCATTTCGTGGCTTTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                         <  1:56437/149‑1 (MQ=11)
                                                                                                                       acttGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                         <  2:3300546/147‑1 (MQ=255)
                                                                                                                        cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGCCAAAGGTCAAATACGCACCCCGCGGATTTTTGCTAATGAGGAGATCGTTGACCGAAAAAAAACAGAAAAAAAGa                        >  1:1024547/1‑150 (MQ=17)
                                                                                                                         aTGAGAATTCCGATCCCCAATTCCTGGCTGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                       <  1:2348150/150‑1 (MQ=255)
                                                                                                                              aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                  >  2:3951504/1‑150 (MQ=255)
                                                                                                                                ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                <  1:2778931/150‑1 (MQ=255)
                                                                                                                                   cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc             <  2:3139851/150‑1 (MQ=39)
                                                                                                                                     aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTAAAAAACGCTCCCGGCGGCTTTGTGCTAAACAGGAGAGCGGTACCCCACAAACAACAGATAAAAAGAAAGGCCCAGTCtt           >  2:1231585/1‑150 (MQ=11)
                                                                                                                                     aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt           <  1:2164209/150‑1 (MQ=38)
                                                                                                                                      tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCATCGCAAACAGAAAAACCGCTACCCGCGGATTTTTCCTAAACAGGGGAGCGTTCCCCGACAAACAAAACATAAAACGAAAGGCCCACTCatt          >  1:1540843/1‑150 (MQ=11)
                                                                                                                                      gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGCTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt          <  2:3530962‑M1/145‑1 (MQ=255)
                                                                                                                                       atttAATTCCTGGCTGTTTCGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc         >  2:3808396‑M1/5‑150 (MQ=255)
                                                                                                                                        cccATTTCCTGGGTGTTTTGGGGGGGGGGCCGGCCCGCCACCCCCCGGGGCGTTGCTTCGGAAAGGTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg        <  1:4552647/150‑1 (MQ=11)
                                                                                                                                          tacttCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCGACCCTCCGGGCCGTTGCGTCGCTAGGTGCAAATCCGCTGCCGGCGGTTTTGTCCTACTCTGGAGAGCGTTCACGGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGAc      <  2:253501/147‑1 (MQ=11)
                                                                                                                                          cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc      <  1:1874143/150‑1 (MQ=33)
                                                                                                                                          cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc      <  2:1866234/150‑1 (MQ=33)
                                                                                                                                           aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt     <  2:943467/150‑1 (MQ=32)
                                                                                                                                           aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCAAAAACGCTCCCGGCGGATTTGTCCGACTCAGGAGCGCGGTCACCCACCAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACt     >  2:1273096/1‑150 (MQ=11)
                                                                                                                                            aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                              <  1:402078/60‑1 (MQ=35)
                                                                                                                                             ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAAAGGTCAAAAACGCTCCCGGGGGGTTTGTGCGGAACAGGCGCGCGGTCACCCCCCAACAACCGAAAAAAAGAAAGGCCCCCTCTTTGGACTg  >  2:2018372/1‑149 (MQ=2)
                                                                                                                                                     |                                                                                                                                              
AGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAG  >  CP009273/359965‑360256

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: