Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,691,283 A→C 5.6% D193E (GAT→GAG PP_3257 ← putative Methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,691,2830AC 5.6% 326.4 / 8.7 108D193E (GAT→GAGPP_3257putative Methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (54/48);  new base C (3/3);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 118 total aligned reads displayed.

GCGGTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGCTTCCAGTG  >  NC_002947/3691134‑3691432
                                                                                                                                                     |                                                                                                                                                     
   gTCCAGTACCACCAGCACTTCCGCGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATcg                                                                                                                                                    >  1:2231399/1‑150 (MQ=255)
        gTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCgacgac                                                                                                                                               >  2:1511093/1‑150 (MQ=255)
          accaccAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGATCTcgacgagac                                                                                                                                             <  2:482843/150‑4 (MQ=255)
            caccaGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCtt                                                                                                                                           <  2:38316/150‑1 (MQ=255)
               caGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGCCGCTtggc                                                                                                                                        >  1:1657570/1‑150 (MQ=255)
                ggCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCTCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtggct                                                                                                                                       >  1:2357855/2‑150 (MQ=255)
                   aCTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtggctggc                                                                                                                                    >  1:1086132/1‑150 (MQ=255)
                     ttCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCAt                                                                                                                                  >  1:200812/1‑150 (MQ=255)
                     ttCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCAt                                                                                                                                  >  2:767890/1‑150 (MQ=255)
                       ccACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGa                                                                                                                                <  2:146993/150‑1 (MQ=255)
                        cACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGCCAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGAc                                                                                                                               <  1:783506/150‑1 (MQ=255)
                        cACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGAc                                                                                                                               <  2:209438/150‑1 (MQ=255)
                          cgttgaGTCCCACTTGCTCCGCCTCGCGCGCATCCACCAGAACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATTTCGACGGCAATGATACGCCCGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCa                                                                                                                             <  2:47876/145‑1 (MQ=255)
                                gTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTc                                                                                                                       <  1:1511093/150‑1 (MQ=255)
                                gTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTc                                                                                                                       >  2:2042416/1‑150 (MQ=255)
                                 tCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTcc                                                                                                                      >  1:2024253/1‑150 (MQ=255)
                                  cccGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCa                                                                                                                     <  2:1509291/150‑1 (MQ=255)
                                   ccGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAg                                                                                                                    >  2:825409/1‑150 (MQ=255)
                                    cGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGc                                                                                                                   <  1:1790643/150‑1 (MQ=255)
                                     gCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCa                                                                                                                  >  2:1521342/1‑150 (MQ=255)
                                     gCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCa                                                                                                                  >  1:2354515/1‑150 (MQ=255)
                                         gcggcgcCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCaaaaa                                                                                                              >  2:2209945/4‑150 (MQ=255)
                                           tgcGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGt                                                                                                            >  1:70562/1‑150 (MQ=255)
                                           tgcGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGt                                                                                                            >  1:1057921/1‑150 (MQ=255)
                                             cgcCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTcc                                                                                                          >  1:1916247/1‑150 (MQ=255)
                                              gcCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCg                                                                                                         >  2:2381350/1‑150 (MQ=255)
                                                 tGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGTCATAACCAGTATTCCAGCAAAAAGTACGcct                                                                                                      >  2:396728/1‑150 (MQ=255)
                                                     gggCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATcga                                                                                                                                                   <  1:1249525/101‑1 (MQ=255)
                                                     gggCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATcga                                                                                                                                                   >  2:1249525/1‑101 (MQ=255)
                                                         aGCCACCAGCACCGCCCGCAGCCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGc                                                                                              <  1:2152524/150‑1 (MQ=255)
                                                           ccaccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCg                                                                                            >  2:2015709/1‑150 (MQ=255)
                                                            caccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGATCTcgacgagac                                                                                                                                             <  2:1648802/100‑4 (MQ=255)
                                                            caccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGATCTcgacgagac                                                                                                                                             >  1:1648802/1‑97 (MQ=255)
                                                             accaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTtg                                                                                                                                          >  2:1264376/1‑102 (MQ=255)
                                                             accaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGaa                                                                                          >  2:159098/1‑150 (MQ=255)
                                                             accaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATCACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGaa                                                                                          >  2:1093461/1‑150 (MQ=255)
                                                              ccaGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTGATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGaaa                                                                                         >  2:2269516/1‑150 (MQ=255)
                                                                   aCCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                                    >  1:991343/1‑150 (MQ=255)
                                                                     cgcccgcAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGa                                                                                  <  1:619354/150‑1 (MQ=255)
                                                                     cgcccgcAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGa                                                                                  <  1:767890/150‑1 (MQ=255)
                                                                         cgcAGGCATCGCCCAAGGCAGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTACGCCTCCTCATGCCGAAACTGCTGATGCa                                                                              >  2:1100802/1‑150 (MQ=255)
                                                                            aGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACg                                                                           >  1:1128294/1‑150 (MQ=255)
                                                                              gCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGc                                                                                                        <  2:1389629/119‑1 (MQ=255)
                                                                              gCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGc                                                                                                        >  1:1389629/1‑119 (MQ=255)
                                                                              gCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGc                                                                         <  2:1086132/150‑1 (MQ=255)
                                                                                 tCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCaaa                                                                      <  2:336503/150‑1 (MQ=255)
                                                                                     ccAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCgaga                                                                                                                                                       <  1:631809/65‑1 (MQ=255)
                                                                                     ccAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTc                                                                  >  1:394250/1‑150 (MQ=255)
                                                                                     ccAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTc                                                                  <  2:1057921/150‑1 (MQ=255)
                                                                                        aGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCggg                                                               <  2:2231399/150‑1 (MQ=255)
                                                                                        aGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCggg                                                               <  1:825409/150‑1 (MQ=255)
                                                                                            ggcggGCTTCTTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTTGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCaa                                                           <  1:1072858/150‑1 (MQ=255)
                                                                                            ggcggGCTTCGTTATGTGCAGGGACGCGGCCATCAATCATGGGCACGTCTCGTCGAGATCGACGAAGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCaa                                                           >  1:1269271/1‑150 (MQ=255)
                                                                                             gcggGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGcctcct                                                                                                   <  2:1796821/109‑1 (MQ=255)
                                                                                             gcggGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGcctcct                                                                                                   >  1:1796821/1‑109 (MQ=255)
                                                                                                ggCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAgg                                                       <  1:396728/150‑1 (MQ=255)
                                                                                                 gCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGc                                                      >  1:57789/1‑150 (MQ=255)
                                                                                                 gCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGc                                                      <  1:2381350/150‑1 (MQ=255)
                                                                                                     cGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTg                                                  <  1:490576/150‑1 (MQ=255)
                                                                                                       ttATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACGAGTATTACAGGAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCa                                                >  1:1748375/1‑150 (MQ=255)
                                                                                                           gtgCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTc                                            <  2:1916247/150‑1 (MQ=255)
                                                                                                            tgCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcg                                           >  1:417220/1‑150 (MQ=255)
                                                                                                            tgCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcg                                           <  1:159098/150‑1 (MQ=255)
                                                                                                             gCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAg                                                        <  2:843192/136‑1 (MQ=255)
                                                                                                             gCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAg                                                        >  1:843192/1‑136 (MQ=255)
                                                                                                              cAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTcgcg                                         >  1:1213905/1‑150 (MQ=255)
                                                                                                                 ggATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGAcc                                      >  1:909265/1‑150 (MQ=255)
                                                                                                                    tGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTc                                   <  1:2015709/150‑1 (MQ=255)
                                                                                                                      aCAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCt                                 >  2:1990159/1‑150 (MQ=255)
                                                                                                                        aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGa                                                                                  <  2:824483/99‑1 (MQ=255)
                                                                                                                        aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGa                                                                                  >  1:824483/1‑99 (MQ=255)
                                                                                                                        aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAg                                                        >  2:1515271/1‑125 (MQ=255)
                                                                                                                        aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAg                                                        <  1:1515271/125‑1 (MQ=255)
                                                                                                                        aaCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGc                               >  2:1836669/1‑150 (MQ=255)
                                                                                                                         aCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGcc                                                                                                       >  1:284773/1‑77 (MQ=255)
                                                                                                                         aCGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGcc                                                                                                       <  2:284773/77‑1 (MQ=255)
                                                                                                                          cGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCg                                                     >  1:145304/1‑126 (MQ=255)
                                                                                                                          cGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCgg                             <  1:289344/150‑1 (MQ=255)
                                                                                                                                tCATTGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCaa                                                                             >  2:705757/1‑96 (MQ=255)
                                                                                                                                tCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCaa                                                                             <  1:705757/96‑1 (MQ=255)
                                                                                                                                tCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCgg                             >  1:463670/1‑144 (MQ=255)
                                                                                                                                tCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCgg                             <  2:463670/144‑1 (MQ=255)
                                                                                                                                  aTGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCcg                     >  1:2395132/1‑150 (MQ=255)
                                                                                                                                     ggCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGa                  <  1:1484528/150‑1 (MQ=255)
                                                                                                                                      gCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGAt                 <  2:2354515/150‑1 (MQ=255)
                                                                                                                                       cACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATa                <  1:1836669/150‑1 (MQ=255)
                                                                                                                                         cGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAgg              <  2:2357855/150‑1 (MQ=255)
                                                                                                                                          gTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCATCTTCGCCGCCCCGCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAggg             <  2:1269271/150‑1 (MQ=255)
                                                                                                                                           tctcGTCCAGATCGACGACGCTTGGCTGGCATGACCATTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATATGCGGTGCACTTCGCGACCGTCCTGCGGGTAACCCGCGATAGGGt            <  1:626190/150‑1 (MQ=255)
                                                                                                                                            gtcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                                    <  1:2157145/76‑1 (MQ=255)
                                                                                                                                            gtcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAActgct                                                                                    >  2:2157145/2‑77 (MQ=255)
                                                                                                                                             tcgtcgATCTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCTCTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGc          >  2:2106745/1‑150 (MQ=255)
                                                                                                                                               gtcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGCtt        >  2:1450017/1‑150 (MQ=255)
                                                                                                                                               gtcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCAGGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCg                                                     <  2:1750561/105‑1 (MQ=255)
                                                                                                                                                tcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGc                                                      >  2:1910009/1‑103 (MQ=255)
                                                                                                                                                tcgAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGc                                                      <  1:1910009/103‑1 (MQ=255)
                                                                                                                                                  gagaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGc                                                 >  2:201500/1‑106 (MQ=255)
                                                                                                                                                    gaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCgg                                                                <  2:1741170/89‑1 (MQ=255)
                                                                                                                                                    gaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCgg                                                                >  1:1741170/1‑89 (MQ=255)
                                                                                                                                                    gaTCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGCTTCCAGt   >  2:1403406/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGTCCAGTACCACCAGCACTTCCACGTGCAGTCCCGCTTGCTGCGCCTGCCGGGCAGCCACCAGCACCGCCCGCAGGCATCGCCCAAGGCGGCGGGCTTCGTTATGTGCAGGGATGACAACGGCAATCATGGGCACGTCTCGTCGAGATCGACGACGCTTGGCTGGCATGACCAGTATTCCAGCAAAAAGTCCGCCTCCTCATGCCGAAACTGCTGATGCAACGGCAAATGTCGGGTCAATAGGCGGTGCACTTCGCGACCGTCCTGCGGGCAACCCGCGATAGGGTGCTTCCAGTG  >  NC_002947/3691134‑3691432

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