Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 9.6% intergenic (‑107/+488) 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,1730CA 9.6% 573.7 / 12.3 206intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (51/125);  new base A (20/9);  total (71/135)
Fisher's exact test for biased strand distribution p-value = 7.82e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360025‑360321
                                                                                                                                                    |                                                                                                                                                    
ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                                      <  1:2819458‑M1/70‑1 (MQ=255)
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                                      >  2:278288/1‑149 (MQ=255)
   cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                   <  2:4265048/149‑1 (MQ=255)
    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                  >  2:1785335/1‑149 (MQ=255)
         aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                             <  1:1723844‑M1/79‑1 (MQ=255)
         aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTAAAATACGCTCCcgg                                                                                                                                             >  1:305766/1‑149 (MQ=255)
           cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                           <  2:282220/149‑1 (MQ=255)
             tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGTTCAAATCCGCTCCCGGCGGa                                                                                                                                         >  1:2164782/1‑149 (MQ=255)
             tCCAGCGCCACCATCCAGGGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                         <  1:3886140/149‑1 (MQ=255)
             cccagataaaccgacaggccccgtctttcgaatgagcatttctttttatttgatccctggggaataATTTCCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                                         >  2:2969237‑M1/67‑149 (MQ=255)
              acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGGTCAAATACGCTcccggcggct                                                                                                                                        <  2:2669947‑M1/84‑3 (MQ=255)
               cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                       <  2:2466099/149‑1 (MQ=255)
                 gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCAACCTTCAAATACGCTCCCGGCGGATTTg                                                                                                                                     >  1:4271232/1‑149 (MQ=25)
                   gcCACCATCCCGTGCAGGAGCTCCGATGATCAGCTGTCAAACTTGAGAATTCCGATCCACAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                   <  2:1028321/149‑1 (MQ=21)
                    aaaacgaaaggcccagtctttcgactgagcctgtcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCCCACCCGGCGGAGTTGTcc                                                                                                                                  <  1:3261061‑M1/90‑1 (MQ=255)
                     caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATGCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                 <  1:3657828/149‑1 (MQ=255)
                       ccgaacggcccagtcgttcgactgagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               >  2:1695901‑M1/57‑149 (MQ=255)
                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAACCCTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               <  2:3026828/149‑1 (MQ=255)
                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCTTCAAAACCGCCCCCCGCGGATTTGTACTaa                                                                                                                               >  1:1254259/1‑148 (MQ=25)
                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               <  2:612541/149‑1 (MQ=255)
                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               <  2:1592934/149‑1 (MQ=255)
                       acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCCACGGTCAAATACGCGCCCGGCGGctttgtgctga                                                                                                                               <  1:3183388‑M1/93‑6 (MQ=255)
                         gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAAAACGCTCCCCGCGGCTTGGACCTACTc                                                                                                                             <  2:2210540‑M1/95‑1 (MQ=255)
                          aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                            <  2:1821947‑M1/96‑1 (MQ=255)
                           aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                           >  2:1724247‑M1/53‑149 (MQ=255)
                             aGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                         <  1:3882332/149‑1 (MQ=21)
                               tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGATCAAATAACCTCCCGGCGGATATGTAATACTCAGgaga                                                                                                                       >  1:257990/1‑149 (MQ=21)
                               cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                       <  2:1017812‑M1/101‑1 (MQ=255)
                                gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGAGCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                      <  2:2097912/149‑1 (MQ=255)
                                  agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                    <  1:634229‑M1/104‑1 (MQ=255)
                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                   <  1:306066/149‑1 (MQ=255)
                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACCGTCAAATCCCCGCCCGGCGGATTTGTCCTACTCAGGAGAGccg                                                                                                                   >  2:1268364/1‑147 (MQ=21)
                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                   >  2:645097/1‑149 (MQ=255)
                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                  >  2:128899/1‑149 (MQ=255)
                                     aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                 <  1:2129820/149‑1 (MQ=255)
                                     aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                 <  1:1415446/149‑1 (MQ=255)
                                      tttcgactgagcctttcgttttatgtgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCa                                                                                                                <  1:2483490‑M1/108‑1 (MQ=255)
                                       ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                               >  1:1821947‑M1/41‑149 (MQ=255)
                                       cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                               <  2:3233116/149‑1 (MQ=255)
                                         cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                             <  2:3723694‑M1/111‑1 (MQ=255)
                                         ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                             <  2:404807/149‑1 (MQ=255)
                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                                            <  1:4315509‑M1/112‑1 (MQ=255)
                                           actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                           <  2:3883477‑M1/113‑1 (MQ=255)
                                             tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                         <  1:681318/149‑1 (MQ=255)
                                             tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                         <  2:4007011‑M1/115‑1 (MQ=255)
                                              gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAAACGCTCCCCGCGGAGTTGTCCGACTCAGGAGAGCGTTCACCCacaaa                                                                                                        >  2:817732/1‑149 (MQ=18)
                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGacaaac                                                                                                       <  2:295492‑M1/117‑1 (MQ=255)
                                                  ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCAACCCTCCGGGCCGTTGCTTCGCAAAGTTAAAATACGCTCCCGGCGGATTTTTACTACTCAGGAGTGCGTTAACCGACAaaaaac                                                                                                    <  1:1428283‑M1/120‑1 (MQ=255)
                                                    ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                  <  1:4007224‑M1/122‑1 (MQ=255)
                                                     tcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGACCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGATCGTTCACCGACAAACAACAGa                                                                                                 >  2:818787‑M1/27‑149 (MQ=255)
                                                       gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                               <  2:3201495‑M1/125‑1 (MQ=255)
                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                             <  2:4378333/149‑1 (MQ=255)
                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACCTTCAAATACGCGCCCGCCGGAGTTGTACAAATCACGAGAGCGTTCACCCACAaaaaacccaaaaa                                                                                             >  1:3896391/1‑142 (MQ=17)
                                                            cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTCAAAAAAGCGCCCGGCCGGTTTTTGCTACACAGGAGAGGGTTCACCGACAAACAACAGATaaaaag                                                                                          >  1:1004222/1‑147 (MQ=17)
                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCCGGCCGAGTTGGCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                          <  2:324589‑M1/130‑1 (MQ=255)
                                                             tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCCGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAgataaaaaga                                                                                         <  2:1949631‑M1/131‑4 (MQ=255)
                                                             tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                                         <  2:1950090‑M1/131‑1 (MQ=255)
                                                             aTGAGAATTCCGATCCCCAATTCCTGGCAGTGTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGa                                                                                         <  2:1300611/149‑1 (MQ=14)
                                                                gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                               >  1:1373184‑M1/16‑124 (MQ=255)
                                                                gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                               <  2:1373184‑M1/109‑1 (MQ=255)
                                                                 atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                     <  2:325217‑M1/135‑1 (MQ=255)
                                                                  tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                                    >  1:3069253‑M1/14‑149 (MQ=255)
                                                                  aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                                    <  2:3843094/149‑1 (MQ=255)
                                                                     tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                 <  1:1660775/149‑1 (MQ=40)
                                                                        gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAAACCGCGACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAAAGAaaggcccccta                                                                              <  1:3252538‑M1/142‑5 (MQ=255)
                                                                         aTCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             <  1:1394915/149‑1 (MQ=38)
                                                                          taattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            >  2:2598054‑M1/6‑149 (MQ=255)
                                                                          tCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            <  1:2777244/149‑1 (MQ=255)
                                                                          tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCAAACGTAAAAAAAGCACCCCGCGGATTTTTGCAAAACAGGGGAGCGTTGACCCGCAAAAAACAGATAAAAAGAAATGGCCaaggtg                                                                            >  1:1310773/1‑144 (MQ=2)
                                                                          gaactAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            >  2:3160060‑M1/6‑149 (MQ=255)
                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATACCCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAACAGATTAAAAGTAAGGCCAAGtagtt                                                                           >  2:2369273/1‑145 (MQ=11)
                                                                            attAATCCCTGGCAGTTTTGGGCGGGGGGCCGGCCCGCCACCCCCCGGGCCGTTGCTTCGCAACGTTCACATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          >  1:326174‑M1/4‑149 (MQ=255)
                                                                             ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                         <  2:4373831‑M1/147‑1 (MQ=255)
                                                                             ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                         >  2:3283950/1‑149 (MQ=34)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCAACCTTAAAAACCGCGCCCCGCGGATTTTTGCAAAACACGAGAGCGTTCACCCACAAACAAAAGATAAAAAGAAAAGCCCCGTCattaggc                                                                       >  1:2849147/1‑145 (MQ=2)
                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                      <  2:1182195/149‑1 (MQ=32)
                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                      <  1:1785335/149‑1 (MQ=32)
                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                      >  2:2878698/1‑149 (MQ=32)
                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                    <  1:4350363/149‑1 (MQ=11)
                                                                                  tCCTGGCAGTGTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                    <  2:3074470/149‑1 (MQ=11)
                                                                                       gctgTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                                               <  1:201700/146‑1 (MQ=1)
                                                                                                        tcgtGCCCCCAACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAAAAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                                              <  1:3452598/146‑1 (MQ=1)
                                                                                                             ccccccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                         <  1:3860175/149‑1 (MQ=14)
                                                                                                             cccGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCACCAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                         <  1:2064063/149‑1 (MQ=9)
                                                                                                              ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                        <  2:2452455/149‑1 (MQ=32)
                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  1:2758710/149‑1 (MQ=33)
                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  2:908274/149‑1 (MQ=33)
                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  1:128899/149‑1 (MQ=33)
                                                                                                                  cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                    <  2:540351/149‑1 (MQ=35)
                                                                                                                    cccTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGACAAAAAACAGATAAAACGAAAGGCCCAGTATTTAGACTTAGCCATTAGTTTTATTTGTTGCCTGgggaattgtg                                  >  2:1317752/1‑147 (MQ=11)
                                                                                                                          gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACGTGGGaattaa                            >  2:1038082/1‑149 (MQ=255)
                                                                                                                             ccGGTGCGTCGCAACGAACAAATCCGCGCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc                         <  1:3222029/149‑1 (MQ=14)
                                                                                                                              cggtGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCCGGAGTGCGGTCACCGACAAACACCCGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                        <  1:3003670/146‑1 (MQ=14)
                                                                                                                               ggtgatgcGCACCGTTCAAAGCCGCTCCCGGCGGTATTTTCCTACTCAGCAGATCGTTCACAGAAAAACAACAGATCAAAGGTAAGGCCCAGTCTTTCTACTGAGCATTTCGTTTTATTTGATGCCTGGGTAATTGGGAATTAATTCtt                       <  1:2210540/142‑1 (MQ=14)
                                                                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGTTGCCTGGGGAATTGGGAATTAATTCtt                       >  1:739994/1‑149 (MQ=255)
                                                                                                                                       gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGa               >  1:2306599/1‑149 (MQ=255)
                                                                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa              >  2:3681640/1‑149 (MQ=255)
                                                                                                                                          aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg            >  2:3180021/1‑149 (MQ=255)
                                                                                                                                               cAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt       <  2:1428283/149‑1 (MQ=255)
                                                                                                                                                 aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg     >  1:841583/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360025‑360321

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