Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,147 G→C 33.4% intergenic (‑81/+514) 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,1470GC33.4% 156.6 / 34.2 79intergenic (‑81/+514)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (32/20);  new base C (24/3);  total (56/23)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAA  >  CP009273/359998‑360291
                                                                                                                                                     |                                                                                                                                                  
gACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                    >  1:1385733/1‑150 (MQ=255)
   aTCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                                                                 <  1:652421/150‑1 (MQ=255)
       ttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGGCTTCCCGCCACCCCCCCTGCCGtt                                                                                                                                             <  2:1341334‑M1/51‑1 (MQ=255)
       agagGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGCCCTGCCCCCCACCCCCCccgccgct                                                                                                                                             >  1:749980/1‑148 (MQ=255)
        tgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                             >  1:123937‑M1/99‑149 (MQ=255)
            ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                         >  1:364253‑M1/95‑149 (MQ=255)
            ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                        >  1:586247‑M1/95‑150 (MQ=255)
            ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                        >  1:638127‑M1/95‑150 (MQ=255)
            ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                        >  2:652464‑M1/95‑150 (MQ=255)
              cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                     <  1:912865/150‑1 (MQ=255)
               ctcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                      >  1:647691‑M1/92‑148 (MQ=255)
                 cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTCGCCCCCCCCCCCCGCGCCGTTGCTTTGCCAc                                                                                                                                  >  2:1384251/1‑150 (MQ=255)
                  aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                   >  1:1458800‑M1/89‑148 (MQ=255)
                  aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                  >  1:342117‑M1/89‑149 (MQ=255)
                  aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                 >  2:1155973‑M1/89‑150 (MQ=255)
                   ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCCGGCCCCTGCTTCGCaaccc                                                                                                                                <  2:452767‑M1/63‑3 (MQ=255)
                   ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCCGGTCGTTGTTTCGCaaaaa                                                                                                                                >  1:557591/1‑147 (MQ=38)
                    cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGtt                                                                                                                              >  1:371759/1‑150 (MQ=255)
                     gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGCCCTCCCCCCCCCCCCCCCCGCCGTTGCTTCCCCACATTc                                                                                                                             >  2:662225/1‑150 (MQ=255)
                      cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                            <  1:1408618/150‑1 (MQ=255)
                        gcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCGGCCCGCCACCCCCCCGGGCGCTGCTTCCgcaacttcca                                                                                                                             <  2:1460566‑M1/67‑6 (MQ=255)
                          gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                        >  1:1006297‑M1/81‑150 (MQ=255)
                          ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGGCTCCCCCCCCCCCCCCCGGCCGTTGGTTTGCTACGTTCAAATc                                                                                                                        >  2:1255826/1‑150 (MQ=38)
                           cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGGTGCCTTGCCACCCTCAAAAcc                                                                                                                       >  2:620579/1‑150 (MQ=255)
                             gCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                     <  2:234010/150‑1 (MQ=255)
                              cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCCGTGCTTCGCAACGTTCAAAtaccct                                                                                                                    >  2:123937/1‑145 (MQ=38)
                               ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGTTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                   <  2:1243472/150‑1 (MQ=255)
                               ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCGGTGCTTCGCCACGTTCAAAAcccccc                                                                                                                   >  2:1268764/1‑146 (MQ=37)
                                 cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                 <  1:1243307/150‑1 (MQ=255)
                                 cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCCGGCCGTTGCCTTGCAAAGTTTCAAATCGCTccc                                                                                                                 >  1:182977/1‑150 (MQ=255)
                                           aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                       <  2:1431589/150‑1 (MQ=255)
                                             ataaaacgaaaggcccagtctttcggctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                     >  2:661800‑M1/62‑150 (MQ=255)
                                               ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCGGGCCGGTGCTTTGCAAAGTTCAAAACCCCTCCCCGCGGGttttttct                                                                                                   >  1:678153/1‑145 (MQ=25)
                                                aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                   >  1:486606‑M1/59‑149 (MQ=255)
                                                 aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                 >  2:1390650‑M1/58‑150 (MQ=255)
                                                   caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGTTGCTTCTCAACCATTAAAAACGCCCCCCGCGGGGTTTTCctcctc                                                                                               >  2:963211/1‑150 (MQ=25)
                                                     tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                             <  2:12792/150‑1 (MQ=255)
                                                      aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                            >  1:452767‑M1/53‑150 (MQ=255)
                                                       aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                           >  1:657387‑M1/52‑150 (MQ=255)
                                                         gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                         <  2:1385733/150‑1 (MQ=255)
                                                          cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                        >  2:305937‑M1/49‑150 (MQ=255)
                                                                aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGATTTTTTCTACTCAGGAGAGGGTTCa                                                                                  >  1:1021711/1‑150 (MQ=21)
                                                                aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCCTTCTTTGCAACGGTCAAATCCCCCCCCCGCGGGGTTGTCCTCCCCCCGAGAGCGTTCa                                                                                  >  2:593412/1‑150 (MQ=21)
                                                                   tcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                  >  1:348976‑M1/40‑147 (MQ=255)
                                                                     cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCTCCCCGCGGGTTTGTCCTACTCCGGGGAGCGTTTACCGa                                                                              >  1:652464/1‑149 (MQ=18)
                                                                      gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCCCCCGGCCGCTGCTTCGCAACCTTCAAAATCCCCCCCCGCGGGGTTTTTCTTCTCCGGAGAGGGtttcccccc                                                                             >  1:701435/1‑142 (MQ=16)
                                                                      actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                            >  1:1212943‑M1/37‑150 (MQ=255)
                                                                         gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCCCTGCTTTCCAACCTTCAAATCCCCCCCCCGCGGGGTTTTTCTTCTCCGGGGAGCGTTCACCCacaaa                                                                          >  1:1320314/1‑149 (MQ=18)
                                                                          agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                        >  1:1002719‑M1/33‑150 (MQ=255)
                                                                          agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                        >  1:1192025‑M1/33‑150 (MQ=255)
                                                                           gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                       >  1:1460566‑M1/32‑150 (MQ=255)
                                                                             ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                     >  1:593412‑M1/30‑150 (MQ=255)
                                                                             ctttcgttttatttgatgcctgggaaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                     >  1:564983‑M1/30‑150 (MQ=255)
                                                                               cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                   <  1:662225/150‑1 (MQ=255)
                                                                               cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCCCCCCCCCCGCCGCTGCTTTGCCACCTTCAAATCCCCCCCCCGCGGGTTTTTTCTTCTCCGGAGAGCGGTCACCCCCAAacaacaa                                                                    >  2:1002719/1‑148 (MQ=12)
                                                                                 cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                  >  1:1315145‑M1/26‑149 (MQ=255)
                                                                                 cgttttatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                  >  2:951800‑M1/26‑149 (MQ=255)
                                                                                    aaaCATGAGAATTCCGTTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                               <  2:93329/149‑1 (MQ=255)
                                                                                    aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGCTGCTTCGCAACGGTCAAAAACCCTCCCCGCGGGTTTGTTCTACTCCGGAGAGCGGTTCCCGACAAAAAAAAGAaaaa                                                               >  1:521223/1‑149 (MQ=17)
                                                                                      tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                             >  2:521223‑M1/21‑149 (MQ=255)
                                                                                      aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                            <  2:424310/150‑1 (MQ=255)
                                                                                        atcagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGTTGCTTTGCAAAGTTTAAATCCGCCCCCCGCGGGTTTTTCCTACTCAGGAGAGCGGTCACCGCCAAACAACAGATAAAAAGa                                                           >  1:678987/4‑149 (MQ=14)
                                                                                        aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCCCGCGGGTTTTTTCTACTCAGGAGAGCGTTCACCGACAAACAACAAATaaaaaca                                                           >  1:548176/1‑146 (MQ=14)
                                                                                        aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCTTGCTTTCCAACCTTCCAATCCCCCCCCCGCGGGGTTTTTCTTCTCCCGAGGGCGTTCTCCCACCAACAACAAATaaaaaca                                                           >  2:348976/1‑146 (MQ=9)
                                                                                           agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                        <  2:695369/149‑1 (MQ=255)
                                                                                            gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCGCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                      <  2:451056/150‑1 (MQ=255)
                                                                                                ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                  >  1:1221516‑M1/11‑150 (MQ=255)
                                                                                                     gaattaaTTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCGTCGCAAAGTTCAAAACCGCCCCCCGCGGGTTTTTTTTACTCACGAGAGCGGTTACCCACAAACAACAAATAAAAAAAAAAGCCCCATCttt                                             >  2:782327/6‑150 (MQ=255)
                                                                                                     gaattaaTTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCCCGCGGCTTTTTTCTTCTCCGGAGAGCGTTCACCGACAAACAACAAATAAAAAGAAAAACCCAGTCctt                                             >  2:1329502/6‑150 (MQ=11)
                                                                                                      ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCCGCCGTTGCTTTTCCACGGTCCAAACCGCCCCCCGCGGGTTTGTTCTTCTCAGGGGAGGGGTCCCCGACAAAAAAAAAAGAAAAAGAAAGGCCCCGCCTTTc                                            >  2:497885/1‑150 (MQ=2)
                                                                                                      ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCTTCAAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGCGTTCACCCACCAAAAACAAATAAAAAGAAAAGCCCCGACTTTc                                            >  1:951800/1‑150 (MQ=11)
                                                                                                      ccccAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCACCCCCCCGGCCCTTGCTTTGCCAAGTTCAAATCCCCCCCCCGCGGGTTTTTCCTTCTCAGGAGAGCGTTCCCCGCCAAACAACAGAAAAAAAGAAAAACCCAGTCtttt                                            >  1:874230/1‑149 (MQ=2)
                                                                                                      aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                            >  1:455501‑M1/5‑150 (MQ=255)
                                                                                                         cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                         <  1:840040/150‑1 (MQ=33)
                                                                                                         cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                         <  2:804193/150‑1 (MQ=33)
                                                                                                         cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                         <  2:331995/150‑1 (MQ=33)
                                                                                                          aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                        <  2:875172/150‑1 (MQ=32)
                                                                                                           aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                       <  1:874018/150‑1 (MQ=32)
                                                                                                                                          cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat        >  1:305937/1‑150 (MQ=33)
                                                                                                                                           gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat        >  2:295605/1‑149 (MQ=33)
                                                                                                                                                ccagcGGGCCGTTGATTCGCAACGTTCAAATCCGTTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa  <  2:6209/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAA  >  CP009273/359998‑360291

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