Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,165 G→C 6.1% intergenic (‑99/+496) 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,1650GC 6.1% 1081.6 / 7.1 350intergenic (‑99/+496)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (104/219);  new base C (18/8);  total (122/228)
Fisher's exact test for biased strand distribution p-value = 3.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGC  >  CP009273/360016‑360314
                                                                                                                                                     |                                                                                                                                                     
 ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattTATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCAGGGCCGTTGGTTAGCACCGt                                                                                                                                                      <  2:720896‑M1/63‑1 (MQ=255)
  gagagcgttcaccgacacacaacagataaacctaaaggcccagtcgttcgcctgagcctttcgcattatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                                     >  2:5181914‑M1/87‑150 (MQ=255)
  cagagcggtcaccgacaaacaacagatcaacccacaggcccagtctttcgccggatcctttcgttttaatttatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                                     >  2:4487018‑M1/87‑150 (MQ=255)
  cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGtt                                                                                                                                                     >  2:4050865/1‑150 (MQ=255)
   agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                    >  2:1617487‑M1/86‑150 (MQ=255)
     ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTCGCAACCATCaa                                                                                                                                                  >  1:2574726/1‑150 (MQ=38)
     agcgttcaccgacacccaccagcgaaaccgacaggcccagtctttcgactgagcctgtcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                                  >  2:781449‑M1/84‑150 (MQ=255)
         cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCCACGTTCAAATcc                                                                                                                                              >  2:3353197/1‑150 (MQ=255)
          aGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGACTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCg                                                                                                                                             >  2:138514/1‑150 (MQ=255)
              gCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                         <  2:4267856/150‑1 (MQ=255)
                ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                       <  2:2108213/150‑1 (MQ=255)
                ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCCACGTTCAAATCCGCTCCcg                                                                                                                                       >  2:4081579/1‑150 (MQ=255)
                   ccATCCCGCGCCACCATCCCGTGCTGGAGCTCCGATTATAAGCTGTAAAACTTGAGAATTCCGATCCCCAATGCCGGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                    <  1:2033318/150‑1 (MQ=21)
                     caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCGACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                  <  2:3026391‑M1/83‑1 (MQ=255)
                      tCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                 <  2:2413695/150‑1 (MQ=255)
                        cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                               <  2:2610752/150‑1 (MQ=255)
                           cgcCACCATCCAGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                            <  2:1832871/150‑1 (MQ=255)
                            taaaacgaaaggcccagtctttcgactgagcctttcgtgttctttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGCTGTGTcc                                                                                                                           <  2:5010191‑M1/90‑1 (MQ=255)
                              cacaaTCCTGTCCTGTAGCTCCTATGATACGCTGAAAAACATCAGAATCCCGACCCCCAATCCCTGGCTGTTTTTGGCGGGCGCCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                         <  1:3141222/150‑1 (MQ=17)
                               aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCGAATACGCGACCGGCGGctttgtgctca                                                                                                                        <  1:224264‑M1/93‑6 (MQ=255)
                                  aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                     <  2:4328055/150‑1 (MQ=255)
                                   aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTGAAATCCGCTCCCGGCGGGGTTGTGCTActccg                                                                                                                    <  2:3642261‑M1/97‑3 (MQ=255)
                                     cggcccaggatttcgactgagcatttctttttattttatgcctgggtaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                  >  1:2027435‑M1/52‑150 (MQ=255)
                                     cgcccccggctttcgactgagccgttcgttttatttggtgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                  >  2:3534691‑M1/52‑150 (MQ=255)
                                     aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATACGCACCCCGCGGAGTTGTACAACTCAGga                                                                                                                  <  2:49389‑M1/99‑1 (MQ=255)
                                        tGCTGGAGCTCCGTTTATAAGCTTTCAAACATGATAAGTCCGATCCCCATTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                               <  1:4980403/150‑1 (MQ=17)
                                         gCAGGAGCTCCGATGATAAGCGGGCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                              <  2:4742226/150‑1 (MQ=255)
                                           aGGAGCTCCGATGATACGCTGTCACACAGGAGAATTCCGATCCCCACTTCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGt                                                                                                            <  2:2021302/150‑1 (MQ=18)
                                            ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                           <  1:1732570/150‑1 (MQ=255)
                                            ggAGCGCCGTTGATACGCTGACAAACATGAGAATGCCGAGCCCCACTTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                           <  1:5014034/150‑1 (MQ=17)
                                               tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                             >  1:391486‑M1/42‑129 (MQ=255)
                                               tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATCCCCACCCCGCGGATTTGGGCTACTACGGAGAGCGTTCAc                                                                                                        <  2:1304576‑M1/109‑1 (MQ=255)
                                                ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                       <  2:4171100‑M1/110‑1 (MQ=255)
                                                cTCCGATGATAAGCTGTCACACATGAGAAGTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                       <  2:2909324/150‑1 (MQ=255)
                                                   gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                    >  2:1553132‑M1/38‑150 (MQ=255)
                                                   gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                    >  2:2616146‑M1/38‑150 (MQ=255)
                                                    gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAAAACGCACCCGGCGGAATTGTCCAAAAAAGGAGAGCGTTAACCcaca                                                                                                   >  1:5239167/1‑150 (MQ=18)
                                                      tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTTTCCTACTCAGGAGAGCGGTAACCGacaaa                                                                                                 <  2:283956‑M1/116‑1 (MQ=255)
                                                       gagcctttcgtattatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCATCGCAACGTTCAAATACGCTACCGGCGGATTTGTCCTACTCAGGATAGCGTTAACCGacaaac                                                                                                <  2:1958555‑M1/117‑1 (MQ=255)
                                                        atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCAACGGTAAAAAACGCACCCCGCGGCTGTTTGCGACTCAGGAGTGCGTTCACCGcaaaaca                                                                                               >  1:122182/1‑150 (MQ=18)
                                                        atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                               >  2:2395193/1‑150 (MQ=255)
                                                        agcatttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                               >  2:122182‑M1/33‑150 (MQ=255)
                                                          cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGATTTGTCCTACACAGGAGAGCGTTCACCCACAaacaac                                                                                             <  2:618152‑M1/120‑1 (MQ=255)
                                                          cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGACGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTGTCCTACTCAGGATAGCGTTCACCGTCAaacaac                                                                                             <  2:4957714‑M1/120‑1 (MQ=255)
                                                          aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCGCCCGGCGGAGTTTTCCTACTAAGGAGGGGGGTCACCGccaaaaaaa                                                                                             >  1:1607351/1‑149 (MQ=17)
                                                           aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                            <  2:1496715/150‑1 (MQ=255)
                                                           aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                            <  2:1572099/150‑1 (MQ=255)
                                                             cTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                          <  2:5106407/150‑1 (MQ=255)
                                                              tGTCAAACATGAGATTTCCGATCCCCATTTCCTGGAAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                         <  1:2504017/150‑1 (MQ=255)
                                                              tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTTCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                         <  2:3996065/150‑1 (MQ=255)
                                                               gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCACCGTTCAAAGACGCTCCCCGCGGATTTGTCCGACGCAGGTGAGCGTTCACCCACAAACAACAGaaa                                                                                        >  1:766135/1‑148 (MQ=17)
                                                                tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGGTAAAATCCGCTCCCGGCGGATTTGTACAACAAAGGAGAGCGTTAACCGACAAAAAACAGATaa                                                                                       >  1:3452758/1‑150 (MQ=17)
                                                                tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTCGCCACCGGCAAATACGCTACCCGGGGCGTTGTGCGTCTCGGGGGGGCGGTGACCGCCAACCAACAGATaa                                                                                       >  1:2710280/1‑150 (MQ=12)
                                                                 ttttattttaggcctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                      >  2:3131710‑M1/24‑150 (MQ=255)
                                                                 ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAACATCAAAACCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATaaa                                                                                      <  2:911130‑M1/127‑1 (MQ=255)
                                                                 cAAACATGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                      <  2:1892567/150‑1 (MQ=255)
                                                                 cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                      <  2:227479/150‑1 (MQ=255)
                                                                 cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                      <  2:1486456/150‑1 (MQ=255)
                                                                  aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATACGCTCCCGGCGGATTTGTGATAAAAAGGATAGCGGTAACCGACAAACAAAAGATaaaa                                                                                     >  1:1236115/1‑150 (MQ=17)
                                                                         gttgcgggggtatttAATTCGTGGCTGTTTAGGGCGGGGGTCCTCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                              >  1:3662591‑M1/16‑150 (MQ=255)
                                                                           aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                            <  2:647937/150‑1 (MQ=255)
                                                                            aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                           <  2:128226/150‑1 (MQ=255)
                                                                             cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                          <  2:3449936‑M1/139‑1 (MQ=255)
                                                                              tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACTTAAAATACGCTCACCGCGGATTTGTGCGACTCACGATAGCGGTGACCGACAAAAAAAACATAAATAGGAAAGCCCac                                                                         >  1:1524688/1‑149 (MQ=9)
                                                                               ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTAGCAAAGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCACGAGAGCGTTCACCGACCAACCACAGATAAAACGAAAGGCCCCGt                                                                        >  2:4909713/1‑150 (MQ=14)
                                                                                ggctcattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                       >  2:3715924‑M1/9‑150 (MQ=255)
                                                                                cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                       <  2:4826264/150‑1 (MQ=39)
                                                                                 gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                      <  2:5240823‑M1/143‑1 (MQ=255)
                                                                                 gctcacaaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                      <  2:1741991/143‑1 (MQ=12)
                                                                                  aTCCCCTATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                     <  2:5073106/150‑1 (MQ=25)
                                                                                  aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                     <  2:2446563/150‑1 (MQ=38)
                                                                                  aTCCCCAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                     <  1:1270297/150‑1 (MQ=38)
                                                                                   tCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                    <  1:1172927/150‑1 (MQ=37)
                                                                                     attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                  >  1:842087‑M1/4‑150 (MQ=255)
                                                                                      ttaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCGTCAAATCCGCGCCCCGCGGATTTGTACGACTCAGGAGAGCGGTACCCCACAAAAAACAAAAAAAAAGAAAGGCCCCGTCTTTAGa                                                                 >  1:4679382/3‑150 (MQ=11)
                                                                                      ccAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                 <  1:2934422/150‑1 (MQ=21)
                                                                                       tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                <  2:3497273‑M1/149‑1 (MQ=255)
                                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACGCAGGAGTGCGTTCACCGACAAACACCAGATAATACGAAAGGCCCCGTCGTTCGACt                                                               >  2:3249566/1‑150 (MQ=11)
                                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCTCCCCGCGGATTTTTCCTACTCAGGAGAGCGTTCACCCACAAACAAAACATAAAAAGAAAGGCCCACGCTTTCTAca                                                               >  1:1665501/1‑149 (MQ=11)
                                                                                        aaTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                               <  1:4280916/150‑1 (MQ=17)
                                                                                          ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                             <  2:3915092/150‑1 (MQ=31)
                                                                                          ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAAAACGCTCCCGGCGGCTTTGTACTACTCAAGAGAGCGTTCACCGACAAACAACAGATAATAAGAAAGGCCCAGTGtttataaagg                                                             >  2:2207769/1‑143 (MQ=11)
                                                                                           tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                          <  1:2100425/104‑1 (MQ=11)
                                                                                                     tttgggggggctTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGttt                                                  <  1:1031007/143‑1 (MQ=0)
                                                                                                                    tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                   <  1:3635476/150‑1 (MQ=11)
                                                                                                                    tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCATTCGTTTTTTTTGATGCCTggg                                   >  2:2719508/1‑150 (MQ=2)
                                                                                                                    tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                   <  1:3373703/150‑1 (MQ=11)
                                                                                                                      cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAACCACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                 <  2:1050448/150‑1 (MQ=12)
                                                                                                                       ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGACTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                <  2:2979311/150‑1 (MQ=17)
                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                               <  1:2395193/150‑1 (MQ=33)
                                                                                                                        cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                               <  1:1995109/150‑1 (MQ=33)
                                                                                                                         gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAAAAGATAAAACGAAATGCCCAGTCTTTAGACTGAGCCTTTCGTTTTATTTGATTCCTGgggactt                              >  2:4773578/1‑150 (MQ=11)
                                                                                                                              ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCCGGGAGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa                         <  1:403721/150‑1 (MQ=11)
                                                                                                                                       cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                               >  1:5117320/1‑135 (MQ=33)
                                                                                                                                       cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                <  1:4430719/150‑1 (MQ=255)
                                                                                                                                       cGTTGCTTCGCAACGGTCTTATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                <  1:3175754/150‑1 (MQ=255)
                                                                                                                                        gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACTTGGGAATTAATTCTTg               >  1:4222299/1‑150 (MQ=255)
                                                                                                                                          tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa             >  1:1809320/1‑150 (MQ=255)
                                                                                                                                             ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGACTTTCGACTGAGCCTTTCGTGTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc          >  2:4615836/1‑150 (MQ=255)
                                                                                                                                                gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa       <  1:4773578/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGC  >  CP009273/360016‑360314

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: