Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,132 T→G 5.4% intergenic (‑66/+529) 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,1320TG 5.4% 502.2 / 4.3 210intergenic (‑66/+529)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (68/125);  new base G (2/15);  total (70/140)
Fisher's exact test for biased strand distribution p-value = 6.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

ACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGAT  >  CP009273/359986‑360274
                                                                                                                                                  |                                                                                                                                              
aCCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTg                                                                                                                                               <  1:3146926/148‑1 (MQ=255)
 ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGc                                                                                                                                              >  2:1613063/1‑148 (MQ=255)
 ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGc                                                                                                                                              <  2:2825694/148‑1 (MQ=255)
  atccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcacagtctttcgactgagactttcgttttatttgatgcctgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGcc                                                                                                                                             >  1:2380665‑M1/117‑148 (MQ=255)
   ttGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                                            <  1:3099907/148‑1 (MQ=255)
   tccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGccc                                                                                                                                            <  1:2239453‑M1/33‑1 (MQ=255)
   tccgctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctgtcgactgagcctttcggggtagttgatgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGccc                                                                                                                                            <  1:183749‑M1/33‑1 (MQ=255)
       ctcccggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGATGTTTATGGCGGGAGTCCTGcccgaca                                                                                                                                        <  2:1302310‑M1/37‑4 (MQ=255)
        gAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAc                                                                                                                                       >  2:64156/1‑148 (MQ=255)
         cgcggcggatttttcctactcacgagagcggtcaccgacaaacaacagatcaaccgaaaggcccagtctgtcgcctgagcctttcgttttatttgatgcgtggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      >  2:501796‑M1/110‑148 (MQ=255)
         cccggcggatttgtcctactcaggagagcgttcaccgacaaacaccagatacaccgaacggcccagtctttcgactgagcctttcgttttatttgctgcctggggaattAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      >  2:2266639‑M1/110‑148 (MQ=255)
         aGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      <  2:2294918/148‑1 (MQ=255)
           cGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                                    >  1:3240495/1‑148 (MQ=255)
           cGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                                    <  2:167004/148‑1 (MQ=255)
                gatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                                                               <  1:537020‑M1/46‑1 (MQ=255)
                gatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                                                               <  1:1250293‑M1/46‑1 (MQ=255)
                 ccAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCCGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                                              >  1:552918/1‑148 (MQ=255)
                 atttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                                              <  1:2434369‑M1/47‑1 (MQ=255)
                   ttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                                                                                                            <  1:1504209‑M1/49‑1 (MQ=255)
                   agagGCACTTCACCGGTTGCCAGCGGCTTACCAGCCAGCCCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACTTGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCg                                                                                                                            <  2:1498155/148‑1 (MQ=39)
                   agagGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                                                                                                            <  2:3197801/148‑1 (MQ=255)
                   agagGCACTTCACCGCTTGCCAGCGGATTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATTATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCg                                                                                                                            <  1:2584653/148‑1 (MQ=39)
                    gagGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                                                                                                           >  2:170472/1‑148 (MQ=255)
                      ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTg                                                                                                                         >  2:2697590/1‑148 (MQ=255)
                      ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTg                                                                                                                         >  2:943569/1‑148 (MQ=255)
                           ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                    <  2:3202260/148‑1 (MQ=255)
                             cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                  <  2:79987/148‑1 (MQ=255)
                               ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                      >  2:1081265‑M1/88‑142 (MQ=255)
                               ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                      <  1:1081265‑M1/55‑1 (MQ=255)
                                     cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                          <  2:3064684‑M1/67‑1 (MQ=255)
                                      gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      >  1:105353‑M1/81‑119 (MQ=255)
                                      gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      <  2:3444330‑M1/39‑1 (MQ=255)
                                      gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAcc                                                                                                                                      >  1:3444330‑M1/81‑119 (MQ=255)
                                       ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                        <  1:2990167‑M1/69‑1 (MQ=255)
                                         caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                      <  2:286772‑M1/71‑1 (MQ=255)
                                          cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTGGCCACGGTCaattgcg                                                                                                     >  1:2122306/1‑143 (MQ=38)
                                          cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                     <  2:1154182/148‑1 (MQ=255)
                                          cGGCTTACCAGCCAGCGCCACCATCCAGTGCAGGAGCTCCGAGGATAAGCTGTCAAACCTGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCGCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                     <  2:2140134/148‑1 (MQ=37)
                                          accgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCGGCCACCCTCCGGGGCGTTGCTTCGCAACGTTCAAATCCg                                                                                                     <  2:1988251‑M1/72‑1 (MQ=255)
                                           ggataaCCATCCCGCGCCACCATCCCGTGCTGGAGCTCCTATTATACGCTTTAAAACATCAGAATTCCGATCCCAATTTCGTGGCTGTTTTGGGGGGGGGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                    <  1:1257761/143‑1 (MQ=16)
                                           ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatacctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                    <  2:1061876‑M1/73‑1 (MQ=255)
                                            cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                   <  2:868223‑M1/74‑1 (MQ=255)
                                             cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                  >  2:3181816/1‑148 (MQ=255)
                                               caaacaacagataaatcgaaaggcccagtctttcgactgagcgttgcgttttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCCAATCCGCTccc                                                                                                <  2:2623394‑M1/77‑1 (MQ=255)
                                               caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTccc                                                                                                <  2:2099410‑M1/77‑1 (MQ=255)
                                                 aacaacagataaaacgaaaggcccagtctttagactgagcctttcgttttatttgatgcctggggaatgAATTCGGGGGAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTGGCGTCGCCACGGTGAAATACCctccccg                                                                                              <  2:1494476‑M1/79‑3 (MQ=255)
                                                   caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCcggcg                                                                                            <  2:739713‑M1/81‑1 (MQ=255)
                                                    tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                           <  1:1455529/148‑1 (MQ=255)
                                                     ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCTAATCCGCTCCCGGCGGa                                                                                          <  2:159160/148‑1 (MQ=255)
                                                     ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                          <  2:686531/148‑1 (MQ=255)
                                                     ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                          <  2:159122/148‑1 (MQ=255)
                                                     acagataccacgcacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATGCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                          >  2:183749‑M1/66‑148 (MQ=255)
                                                         cgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                      >  1:2893701/1‑148 (MQ=255)
                                                         ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAAAAAGCTCCCCGCGGatttt                                                                                      <  1:2813139‑M1/87‑2 (MQ=255)
                                                           ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACCTGAGAATTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                    <  1:203376/148‑1 (MQ=255)
                                                            caccaTCCCGTGCAGGAGCTCCGATGATAAGCTGACAAACATGATAATTCCGAACCCCAATTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                   <  1:1698517/148‑1 (MQ=21)
                                                               cgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGGCGTGGGGTCGCAACGTTCCAATCCGCTCCCGGCGGCTTtgtccgtc                                                                                <  2:2952171‑M1/93‑4 (MQ=255)
                                                                aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTAGCAACGTTCAAATACGCTACCGGCGGATTTGTCCTaaa                                                                               >  1:807661/1‑146 (MQ=255)
                                                                 ccccgTGCAGGAGCTCCGATGATACGCTGACAACCATGATAATTCCGATCCCCATTTCCTGGCTGTTTATGGGGGGGGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                              <  1:3347169/144‑1 (MQ=18)
                                                                   aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTGCTACTCAg                                                                            <  2:364413‑M1/97‑1 (MQ=255)
                                                                    cgTGCAGGAGCTCCGATGATAAGCTGTCAATCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCGGCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCCgg                                                                           >  2:3489326/2‑148 (MQ=21)
                                                                        ccgtctctcgactgagcatttcgttttttttgttgcctggttaaatATTTCCTGGCTGTTGATGGCGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                       >  1:1603046‑M1/47‑148 (MQ=255)
                                                                        cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                                                  <  1:3453152‑M1/43‑1 (MQ=255)
                                                                        cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                         <  2:546643‑M1/100‑1 (MQ=255)
                                                                        cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGagag                                                                       <  2:1604422‑M1/102‑1 (MQ=255)
                                                                          gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                                                                                                   <  2:2025929‑M1/42‑1 (MQ=255)
                                                                          gtctttcgacggagcctttcgttttatttgaggcctgggtaattACTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                     >  2:892451‑M1/45‑148 (MQ=255)
                                                                          ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                     <  1:2517913/148‑1 (MQ=255)
                                                                           gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                    >  2:446863/1‑148 (MQ=255)
                                                                            ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                   >  2:1504216‑M1/43‑148 (MQ=255)
                                                                            aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                   <  2:325920/148‑1 (MQ=255)
                                                                             gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACGTTACAATCGGCTCCCGGCGGCTTTTTCCTGCTAACGATAGCGtta                                                                  >  2:2380665/1‑147 (MQ=18)
                                                                                 ctatgatAAGCTGTCAAACATGAGAATTCCTATCCCCAAGGCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                              <  1:469638/146‑1 (MQ=255)
                                                                                    tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATCCGCTCCCGGCGGATATGTCCTACAAAGGAGAGAGTTCACCGaca                                                           >  2:2250496/1‑148 (MQ=18)
                                                                                     gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                              <  1:1752911/144‑1 (MQ=255)
                                                                                          gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCACCGGTCAAATACCCTCCCCGCGGCTTTTTGCTACTCAGGGGAGCGGTGGCCGccaaacaaa                                                     >  1:680384/1‑147 (MQ=17)
                                                                                           ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACTTGAAAACCGCGACCGGCGGGTTTTTCCGAAACAGGAGCGCGTTCACCGACAAacaaca                                                    <  1:2753802‑M1/121‑1 (MQ=255)
                                                                                            tcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaacac                                                   <  1:2880331‑M1/122‑2 (MQ=255)
                                                                                             gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCGTCAAAACCGCGCCCCGCGGCGGTGGGCGACAAAGGAGCGCGTTCACCCACCAaccacaca                                                  >  1:3326005/1‑146 (MQ=12)
                                                                                             gTCAAACAGGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                  <  1:870379/148‑1 (MQ=255)
                                                                                                tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                               >  1:2011294‑M1/23‑148 (MQ=255)
                                                                                                    aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                           <  2:3327528/148‑1 (MQ=255)
                                                                                                       agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                        <  2:1065844/148‑1 (MQ=255)
                                                                                                         aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                      <  1:3501685/148‑1 (MQ=255)
                                                                                                             ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                  <  2:1695788/148‑1 (MQ=40)
                                                                                                             ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                  <  2:2475868/148‑1 (MQ=40)
                                                                                                               gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                                                                                                   <  1:2718533‑M1/42‑1 (MQ=255)
                                                                                                               gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                                                                                                   >  2:2718533‑M1/8‑49 (MQ=255)
                                                                                                                aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCATCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                               <  1:1665943/148‑1 (MQ=38)
                                                                                                                aTCCCCAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                               <  2:299380/148‑1 (MQ=38)
                                                                                                                 taattAATGCCTGGCAGTTTTTGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                              >  1:996623‑M1/6‑148 (MQ=255)
                                                                                                                 tCCCCATTTCCTGGCTGTTTTGGGCGGGCGCCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                              <  1:1644467/148‑1 (MQ=11)
                                                                                                                 tCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                              <  2:2570812/148‑1 (MQ=21)
                                                                                                                 tCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                              <  2:3258766/148‑1 (MQ=21)
                                                                                                                  ccccAATTCCTGGATGTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                             <  1:467983/148‑1 (MQ=11)
                                                                                                                   attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                            <  2:3061020‑M1/145‑1 (MQ=255)
                                                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTGCTACACAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCACTCTTTCGTCTg                      >  2:310672/1‑148 (MQ=11)
                                                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAACCGCTCCCCGCGGATTTGTGCTAATCAGGCGAGCGGTCACCGACAAACAAAACATAAAACGAAAGGCCCACTCTTTCGACTg                      >  1:2265686/1‑148 (MQ=11)
                                                                                                                                   ttttGGGGGGGGTCCTGCACGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGt            <  1:2797717/145‑1 (MQ=0)
                                                                                                                                    tagggcgggcgTCCGCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGtt           <  1:943569/145‑1 (MQ=0)
                                                                                                                                                  |                                                                                                                                              
ACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGAT  >  CP009273/359986‑360274

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: