Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,064:1 +C 11.4% coding (103/822 nt) PP_3760 ← methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,289,0641.C11.4% 312.8 / 36.1 88coding (103/822 nt)PP_3760methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (38/40);  new base C (5/5);  total (43/45)
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.22e-01

TCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGC  >  NC_002947/4288919‑4289212
                                                                                                                                                   |                                                                                                                                                    
tCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATc                                                                                                                                                  <  2:1375628/150‑1 (MQ=255)
   gggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGaa                                                                                                                                               <  2:597058/150‑1 (MQ=255)
        aCGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                          <  2:1037098/150‑1 (MQ=255)
        aCGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGcc                                                                                                                                          <  2:2214582/150‑1 (MQ=255)
        aCGGAACATCTCGCTGACCGGGATCGCTTCCTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGcc                                                                                                                                          >  1:1871190/1‑150 (MQ=255)
          ggAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCgg                                                                                                                                        >  2:1811533/1‑150 (MQ=255)
           gAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                       <  2:1253382/150‑1 (MQ=255)
                tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                          >  1:795371/1‑142 (MQ=255)
                tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                          <  2:795371/142‑1 (MQ=255)
                 ctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGTCGGAATAAtt                                                                                                                                 >  2:1454863/1‑150 (MQ=255)
                  tcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTg                                                                                                                                >  2:668988/1‑150 (MQ=255)
                   cGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGc                                                                                                                               <  1:1888555/150‑1 (MQ=255)
                    gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGGGCATCTACCCTAGCACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                              <  2:1871190/150‑1 (MQ=255)
                    gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                              <  2:1409268/150‑1 (MQ=255)
                    gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                              <  1:1302849/150‑1 (MQ=255)
                        aCCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCCAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGaaa                                                                                                                          >  2:37313/1‑150 (MQ=255)
                             gATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                                     <  1:1811533/150‑1 (MQ=255)
                               tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGc                                                                                                                   >  1:258265/1‑150 (MQ=255)
                                 gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                       <  1:2408391/128‑1 (MQ=255)
                                 gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                       >  2:2408391/1‑128 (MQ=255)
                                  tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                <  2:1539787/150‑1 (MQ=255)
                                  tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                >  1:863004/1‑150 (MQ=255)
                                  acAGATACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                >  2:705798/2‑150 (MQ=255)
                                  aaaGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                >  1:84818/3‑150 (MQ=255)
                                        cgtggagcgACTGCATGAACATGCCGGGGTCGTGCAGGACGCGC‑TCCTGCAACGCGGACCCCGTCAGTCAATCGAACTGACGCAGCGCGTGGAGAGTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCGGTACTTGa                                                                                                          <  2:47575/145‑1 (MQ=255)
                                        aCTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                          <  1:432363/150‑1 (MQ=255)
                                        aCTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                          >  2:462556/1‑150 (MQ=255)
                                         cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGA‑TCGTGTAAATCGCACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                         <  2:64372/150‑1 (MQ=255)
                                          tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                          >  2:73828/1‑148 (MQ=255)
                                          tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                          <  1:73828/148‑1 (MQ=255)
                                          tGCAGAAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAgg                                                                                                        >  2:1722686/1‑150 (MQ=255)
                                            cagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTa                                                                                                      >  2:2086801/1‑150 (MQ=255)
                                                  cTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                          <  1:869277/108‑1 (MQ=255)
                                                  cTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                          >  2:869277/1‑108 (MQ=255)
                                                     cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                              <  2:2086157/149‑1 (MQ=255)
                                                     cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                              >  1:2086157/1‑149 (MQ=255)
                                                     cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                             >  2:1311578/1‑150 (MQ=255)
                                                     cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                             >  1:1086988/1‑150 (MQ=255)
                                                      gtGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGCTCGCCTCg                                                                                            >  1:348916/2‑150 (MQ=255)
                                                       tGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GGTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                           >  1:836908/1‑150 (MQ=255)
                                                       tGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                           <  2:912961/150‑1 (MQ=255)
                                                        gAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGAt                                                                                          >  1:1606948/1‑150 (MQ=255)
                                                         aaCATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATc                                                                                         <  2:258265/150‑1 (MQ=255)
                                                                ccGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc                                                                                  <  2:1606948/150‑1 (MQ=255)
                                                                 cGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                   <  1:2268620/100‑1 (MQ=255)
                                                                 cGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                   >  2:2268620/1‑100 (MQ=255)
                                                                      cGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATGGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                            <  2:84818/150‑1 (MQ=255)
                                                                        tGCAGGA‑GCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCg                                                                          >  2:607726/1‑150 (MQ=255)
                                                                         gCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGa                                                                         >  1:1901974/1‑150 (MQ=255)
                                                                              gagcgcgtCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCCCGATGTc                                                                     >  1:1243419/8‑150 (MQ=255)
                                                                              aCGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GGTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCg                                                                    <  2:836908/150‑1 (MQ=255)
                                                                                   c‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                               <  2:1086988/150‑1 (MQ=255)
                                                                                        tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACtt                                                                                                            >  1:83381/1‑101 (MQ=255)
                                                                                        tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACtt                                                                                                            <  2:83381/101‑1 (MQ=255)
                                                                                        tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                             >  2:138981/1‑116 (MQ=255)
                                                                                        tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                             <  1:138981/116‑1 (MQ=255)
                                                                                              gCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                                     <  1:1311578/150‑1 (MQ=255)
                                                                                              gCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                                     <  1:462556/150‑1 (MQ=255)
                                                                                                aGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTAGCTAGt                                                   >  2:1116579/1‑150 (MQ=255)
                                                                                                     ccGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                              <  1:2086801/150‑1 (MQ=255)
                                                                                                          aGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCt                                         >  1:2012687/1‑150 (MQ=255)
                                                                                                           gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGa                                                                               <  1:1122601/111‑1 (MQ=255)
                                                                                                           gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                        >  1:2361938/1‑150 (MQ=255)
                                                                                                           gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                        <  1:579632/150‑1 (MQ=255)
                                                                                                           gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTAGGa                                                                               >  2:1122601/1‑111 (MQ=255)
                                                                                                                tCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                    >  1:1133041/1‑69 (MQ=255)
                                                                                                                tCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                    <  2:1133041/69‑1 (MQ=255)
                                                                                                                 cGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTg                                  >  1:1693592/1‑150 (MQ=255)
                                                                                                                   aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCg                                >  1:1583187/1‑150 (MQ=255)
                                                                                                                   aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCg                                >  1:2072852/1‑150 (MQ=255)
                                                                                                                    aCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCAAGTTGCgg                               <  2:1901974/150‑1 (MQ=255)
                                                                                                                         cGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGcc                          <  2:348916/150‑1 (MQ=255)
                                                                                                                          gCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGTGGCAGCTAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCCCCAGTTGCGGCAGCCa                         >  1:150731/1‑150 (MQ=255)
                                                                                                                          gCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                         <  2:2012687/150‑1 (MQ=255)
                                                                                                                             gcgcGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGac                      >  2:1416789/1‑150 (MQ=255)
                                                                                                                               gcgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGAt                                                                        <  1:19076/98‑1 (MQ=255)
                                                                                                                               gcgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGAt                                                                        >  2:19076/1‑98 (MQ=255)
                                                                                                                                cgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACg                   <  1:1688057/150‑1 (MQ=255)
                                                                                                                                  tGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGTTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTACAGTTACGGCCACCAGACACGGa                 >  2:357723/1‑150 (MQ=255)
                                                                                                                                   ggAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGAt                <  2:1693592/150‑1 (MQ=255)
                                                                                                                                    gagaATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATc               >  1:859968/1‑150 (MQ=255)
                                                                                                                                         tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACga          <  1:1416789/150‑1 (MQ=255)
                                                                                                                                         tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACga          <  1:1559345/150‑1 (MQ=255)
                                                                                                                                           cGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCCCCAGTTGCGGCAGCCAGACACGGATCAACgaga        >  2:479256/1‑150 (MQ=255)
                                                                                                                                            ggCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGaa       >  1:713170/1‑150 (MQ=255)
                                                                                                                                            ggCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGaa       <  1:705798/150‑1 (MQ=255)
                                                                                                                                                cTT‑GATCGAAGCGTCGGAATAATTGCGGAAATCGTATCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAgg   <  1:1454863/150‑1 (MQ=255)
                                                                                                                                                cTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAgg   <  2:2072852/150‑1 (MQ=255)
                                                                                                                                                 tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                            <  1:938243/76‑1 (MQ=255)
                                                                                                                                                 tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                            >  2:938243/1‑76 (MQ=255)
                                                                                                                                                 tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc  >  2:492447/1‑150 (MQ=255)
                                                                                                                                                 tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc  >  2:993092/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGC  >  NC_002947/4288919‑4289212

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: ‑