Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,156 C→G 7.5% intergenic (‑90/+505) 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,1560CG 7.5% 675.2 / 14.3 226intergenic (‑90/+505)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (61/146);  new base G (14/5);  total (75/151)
Fisher's exact test for biased strand distribution p-value = 1.96e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAA  >  CP009273/360008‑360303
                                                                                                                                                    |                                                                                                                                                   
 gCACGTCCCCGCTTGCCAGCGGCTTACCCGCCAGCGCCACCCGCCAGTGCAGGAGCACCGATGAAACGCTGTCAAACAGGACAATTCCGATCCCCACGTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                    <  2:2620494/149‑1 (MQ=25)
  ctactcaggagtgcgttcaccgacaaacaacagatacaccgaaaggccctgtgtttctactgaccctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   >  2:4148552‑M1/95‑149 (MQ=255)
  cACTTCACCGCTTGCCAGCGGCTTCCCATCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTTTCAAACATGAGAATTCCGATCCCAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   <  1:1419013/149‑1 (MQ=255)
  cACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   <  2:3242613/149‑1 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGc                                                                                                                                                >  1:1695901/1‑149 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                                >  2:2758710/1‑149 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGc                                                                                                                                                >  1:3844638/1‑149 (MQ=255)
      tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                               <  2:635492‑M1/59‑1 (MQ=255)
           agagcgttccccgacacacaacagctaaaccgacaggcccagtctttcgcctgcgcctgtcgttttatttgatgcctggggaattATTTCCGGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                          >  2:863417‑M1/86‑149 (MQ=255)
             ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTAGCAAccataa                                                                                                                                        >  1:1480688/1‑144 (MQ=255)
             agcgttcacccacaaccaccagctaaaacgacaggcccagtctttcgcctgagcctttctttttattgtatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                        >  2:3510275‑M1/84‑149 (MQ=255)
               gCCAGCGGCTTACCATCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                      <  2:2081826/149‑1 (MQ=255)
                gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTATATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGGTGAAAt                                                                                                                                     <  1:1215769‑M1/69‑1 (MQ=255)
                 cAGCGGCTTACCAACCCGCGCCACCATCCCGTGCAGGAGCTCCTTTTATATGGTGAAAAACATGCTAATTCCGACCCCAAATCCCTGGCATTTTTGGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                    <  1:1296124/149‑1 (MQ=16)
                   caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                  <  2:4375161‑M1/72‑1 (MQ=255)
                     ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                >  2:533868‑M1/76‑149 (MQ=255)
                        acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACTTTCAAATCCGCGccc                                                                                                                             <  1:1736170‑M1/77‑1 (MQ=255)
                          aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCAACcgg                                                                                                                           >  1:305775/1‑149 (MQ=255)
                           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                          <  2:3069253‑M1/80‑1 (MQ=255)
                            cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                         <  2:282220/149‑1 (MQ=255)
                              tCCAGCGCCACCATCCAGGGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                       <  1:3886140/149‑1 (MQ=255)
                              tCCAGCGCCACCAGCCAGGGCCGGAGCTCCTATGATAAGCTGTCACCCATGAGATTTCCGCGCCCCAATGCCTGGCTGTTTATGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                       <  2:2969132/149‑1 (MQ=18)
                                cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                     <  2:2466099/149‑1 (MQ=255)
                                cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                     <  2:1057398/149‑1 (MQ=255)
                                  gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAAAACGCTCCCGGCGGCTTTg                                                                                                                   >  1:4271106/1‑149 (MQ=25)
                                   cgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                  <  2:2363820/149‑1 (MQ=255)
                                    gcCACCATCCCGTGCAGGAGCTCCGATGATCAGCTGTCAAACTTGAGAATTCCGATCCACAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                 <  2:1028321/149‑1 (MQ=21)
                                     ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTGAAATCCGCTCCCGGCGGATTTGtgc                                                                                                                >  1:2100904/1‑147 (MQ=255)
                                     aaaacgaaaggcccagtctttcgactgagcctgtcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCCCACCCGGCGGAGTTGTcc                                                                                                                <  1:3261061‑M1/90‑1 (MQ=255)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCTCCCGGCGGTATTGTCCt                                                                                                               >  1:4348999/1‑149 (MQ=255)
                                        ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                             <  2:1592934/149‑1 (MQ=255)
                                        ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                             <  2:612541/149‑1 (MQ=255)
                                        ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACCTTCAAATACGCTCCCCGCGGAGTTGTCCTAc                                                                                                             >  1:1382203/1‑149 (MQ=25)
                                        acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCCACGGTCAAATACGCGCCCGGCGGctttgtgctga                                                                                                             <  1:3183388‑M1/93‑6 (MQ=255)
                                          aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                           <  2:436872/149‑1 (MQ=255)
                                          aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCAGCTCCCGGCGGATTTGTCCTACTc                                                                                                           <  2:1782705/149‑1 (MQ=255)
                                           aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                          <  2:1821947‑M1/96‑1 (MQ=255)
                                                tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGATCAAATAACCTCCCGGCGGATATGTAATACTCAGgaga                                                                                                     >  1:257990/1‑149 (MQ=21)
                                                 gccggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACCTTCAAAAACGCACCCCGCGGATTTGTCCGACGCAGgaggg                                                                                                    >  1:3866592/4‑147 (MQ=21)
                                                      aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                               <  1:1415476/149‑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)
                                                            actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                         <  2:3883477‑M1/113‑1 (MQ=255)
                                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                       <  1:681318/149‑1 (MQ=255)
                                                              tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                       >  2:1723844‑M1/35‑149 (MQ=255)
                                                               gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCGCCCGGCGGGTTTGTCCTACTCAGGAGAGCGTTGACCGccaaa                                                                                      >  2:3023359/1‑149 (MQ=18)
                                                                     ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                <  1:4007224‑M1/122‑1 (MQ=255)
                                                                      tcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGACCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGATCGTTCACCGACAAACAACAGa                                                                               >  2:818787‑M1/27‑149 (MQ=255)
                                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                              <  2:3241591/149‑1 (MQ=255)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                             <  2:3201495‑M1/125‑1 (MQ=255)
                                                                         cAAACTTGAGATTTCCGATCCCCATTTCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                            <  1:2138529/149‑1 (MQ=14)
                                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACCTTCAAATACGCGCCCGCCGGAGTTGTACAAATCACGAGAGCGTTCACCCACAaaaaacccaaaaa                                                                           >  1:3896391/1‑142 (MQ=17)
                                                                           aaCATGAGATTTCCGATCCCCACTTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                          <  1:2406801/149‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAAACTTAAAATACCCACCCGGCGGTGTTTTGCGAAAACAGTGAGCGTTAACCCACAAAAAACAGATaaaa                                                                          >  2:752226/1‑149 (MQ=12)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTCAAAAAAGCGCCCGGCCGGTTTTTGCTACACAGGAGAGGGTTCACCGACAAACAACAGATaaaaag                                                                        >  1:1004222/1‑147 (MQ=17)
                                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACGCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                        <  2:324616‑M1/130‑1 (MQ=255)
                                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATCCCCTCCCCGCCGATTTGTCCGACACACGAGAGCGTTCACCCACAAACAACACAAAAAacc                                                                        <  1:2912872‑M1/130‑2 (MQ=255)
                                                                              tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCCGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAgataaaaaga                                                                       <  2:1949631‑M1/131‑4 (MQ=255)
                                                                              aggagaAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                       <  2:1300609/147‑1 (MQ=255)
                                                                              aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                       <  2:2719994/149‑1 (MQ=255)
                                                                                tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACACAAAAAACGaaa                                                                     <  2:1277382‑M1/133‑1 (MQ=255)
                                                                                gttatttCCGATACAAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                     <  2:2100904/144‑1 (MQ=255)
                                                                                  gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGGTGGCCACGGTGAAAAAACCTACCCGCGGGTTTGTGCGAATGAGGCGAGCGTTGACCGAGCAAAAACACAAAAAACGgaatg                                                                   >  1:3739242/1‑147 (MQ=9)
                                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGGTTTGTCCTACGCAGGGGCGCGGTCACCGACAACCAACAGATAAAACGAAAgg                                                                   <  2:326174‑M1/135‑1 (MQ=255)
                                                                                   tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                  >  1:3069253‑M1/14‑149 (MQ=255)
                                                                                    aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                 <  2:2786968/149‑1 (MQ=255)
                                                                                     gatgtgtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                >  2:1310773‑M1/12‑149 (MQ=255)
                                                                                      tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                               <  2:1325889/149‑1 (MQ=40)
                                                                                      tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                               <  1:1660775/149‑1 (MQ=40)
                                                                                         gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAAACCGCGACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAAAGAaaggcccccta                                                            <  1:3252538‑M1/142‑5 (MQ=255)
                                                                                           taattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                          >  2:2598054‑M1/6‑149 (MQ=255)
                                                                                           gaactAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                          >  2:3160060‑M1/6‑149 (MQ=255)
                                                                                           cccccAATTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                          <  1:913294/148‑1 (MQ=11)
                                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTGCACCGCCACACAACAGATAAAACGAAAGGCCCCGTCttt                                                         <  2:3953815‑M1/145‑1 (MQ=255)
                                                                                             cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                        <  2:343007/149‑1 (MQ=35)
                                                                                             attAATCCCTGGCAGTTTTGGGCGGGGGGCCGGCCCGCCACCCCCCGGGCCGTTGCTTCGCAACGTTCACATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                        >  1:326174‑M1/4‑149 (MQ=255)
                                                                                              ttaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGATTTTTACTACAAAGGAGAGCGTTAACCGACAAAAAAAAAATAAAACTAAAGGCCCCGTCTTTCg                                                       >  1:2178522/3‑149 (MQ=11)
                                                                                              ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAAAACGCTCCCGGCGGGTTTGTCCAAAAAAGGAGAGCGTTGAACGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCg                                                       >  1:305100/1‑149 (MQ=11)
                                                                                               tgaggcgtGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                      <  1:2706809/142‑1 (MQ=11)
                                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCCGTCTTTCGa                                                      <  2:2937311‑M1/148‑1 (MQ=255)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTACTACTCAGGAGAGCGTTAACCCACAAACAACAAAAAAAACGAAAGGCCCAGTCTTTAGAc                                                     >  1:2849131/1‑149 (MQ=11)
                                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTAGACTg                                                   >  2:1407629/1‑149 (MQ=14)
                                                                                                   tCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                  <  2:1185296/149‑1 (MQ=11)
                                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                     >  2:3191435/1‑130 (MQ=11)
                                                                                                        gctgTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                             <  1:201700/146‑1 (MQ=1)
                                                                                                                gggggggcgTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                     <  1:1268364/145‑1 (MQ=0)
                                                                                                                         tcgtGCCCCCAACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTACCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAAAAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                            <  1:3452598/146‑1 (MQ=1)
                                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCAACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                       <  1:4284605/149‑1 (MQ=14)
                                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                      <  2:2452455/149‑1 (MQ=32)
                                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                    <  2:908274/149‑1 (MQ=33)
                                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                    <  1:2758972/149‑1 (MQ=33)
                                                                                                                                     cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg                >  2:2047182/1‑149 (MQ=37)
                                                                                                                                           gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACGTGGGaattaa          >  2:1038082/1‑149 (MQ=255)
                                                                                                                                            ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat         <  1:1407629/149‑1 (MQ=255)
                                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt        <  1:4373831/149‑1 (MQ=255)
                                                                                                                                              ccGGTGCGTCGCAACGAACAAATCCGCGCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc       <  1:3222029/149‑1 (MQ=14)
                                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGTTGCCTGGGGAATTGGGAATTAATTCtt     >  1:739994/1‑149 (MQ=255)
                                                                                                                                                   gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGTGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa  >  2:2835295/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAA  >  CP009273/360008‑360303

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: