Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,886,644 C→T 6.6% D423D (GAC→GAT rmuC → DNA recombination limiting inversions protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,886,6440CT 6.6% 267.4 / 6.4 92D423D (GAC→GATrmuCDNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base T (3/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCTGCT  >  NC_002947/1886496‑1886793
                                                                                                                                                    |                                                                                                                                                      
tGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCCGCTCGTTCGCGAAAGGGCCGAGCGGGCCTGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACg                                                                                                                                                       <  1:2100113/150‑1 (MQ=255)
  aTTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGGTATAGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCCCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGTCCTGGATGaa                                                                                                                                                     <  2:2117319/150‑1 (MQ=255)
    tGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACt                                                                                                                                                   <  2:851396/150‑1 (MQ=255)
      tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGCCTAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                                 <  1:194315/150‑1 (MQ=255)
      ggtgCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                                 <  2:869802/147‑1 (MQ=255)
       cAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAcagca                                                                                                                                                <  2:105905/150‑5 (MQ=255)
        aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                                              <  2:1148698/150‑1 (MQ=255)
         gCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                             <  2:1060580/150‑1 (MQ=255)
         gCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                             >  1:1060580/1‑150 (MQ=255)
           ccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGc                                                                                                                                           >  2:1462470/1‑150 (MQ=255)
                cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCtt                                                                                                                                      <  1:1734869/150‑1 (MQ=255)
                   gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgca                                                                                                                                   <  2:2363799/150‑1 (MQ=255)
                       cTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                               <  1:1697311/150‑1 (MQ=255)
                          gCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTTGGCAGCCGCTTgcagca                                                                                                                                <  2:486472/146‑1 (MQ=255)
                          gCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTTGGCAGCCGCTTgcagca                                                                                                                                >  1:486472/1‑146 (MQ=255)
                            cACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACg                                                                                                                                                       >  2:1031384/1‑122 (MQ=255)
                            cACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACg                                                                                                                                                       <  1:1031384/122‑1 (MQ=255)
                              ccTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                             <  1:2089269/145‑1 (MQ=255)
                              acTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                             >  2:2089269/2‑145 (MQ=255)
                                      gTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTAt                                                                                                                >  1:470349/1‑150 (MQ=255)
                                        tATCGACAGCCTGTGGAAGTAGGAGCGTAAATGGCAGATCTTTTTCGAAATTGCCGAGCTGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCtttttg                                                                                                              <  1:1892272/150‑5 (MQ=255)
                                          tcggcAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                            <  1:1628534/130‑1 (MQ=255)
                                          tCGGCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                            >  2:1628534/1‑134 (MQ=255)
                                          tCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCgg                                                                                                            <  2:541318/150‑1 (MQ=255)
                                             aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                                 >  2:1241072/1‑111 (MQ=255)
                                             aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                                 <  1:1241072/111‑1 (MQ=255)
                                             aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAg                                                                                                         >  1:1732804/1‑150 (MQ=255)
                                             aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAg                                                                                                         >  2:2272269/1‑150 (MQ=255)
                                                gCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTCCAGCAGGTGGACAAGGCTTATGCGGCAgcgc                                                                                                      >  2:506007/1‑150 (MQ=255)
                                                gCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATCGACGAGCGGGTGGGCTGGCTCTACGACAAGATCATGCTGTTCATCCAGGAACATAGGAAATGGGTCATCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgc                                                                                                      <  1:1252476/150‑1 (MQ=255)
                                                   tgtgGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGca                                                                                                   >  2:788869/1‑150 (MQ=255)
                                                     tgGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaac                                                                                                 >  1:1655145/1‑150 (MQ=255)
                                                       gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaacaa                                                                                               >  1:1040786/1‑150 (MQ=255)
                                                         aGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGc                                                                                             >  2:1541311/1‑150 (MQ=255)
                                                            aGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgt                                                                                          <  2:1315722/150‑1 (MQ=255)
                                                            aGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAAGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgt                                                                                          >  2:325147/1‑150 (MQ=255)
                                                            aGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCAAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgt                                                                                          >  1:1350312/1‑150 (MQ=255)
                                                             ggAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgtg                                                                                         <  1:2326467/150‑1 (MQ=255)
                                                             ggAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgtg                                                                                         >  2:1176686/1‑150 (MQ=255)
                                                                  gTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAg                                                                                    >  2:447122/1‑150 (MQ=255)
                                                                      aGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                             <  1:2280863/89‑1 (MQ=255)
                                                                      aGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                             >  2:2280863/1‑89 (MQ=255)
                                                                       gggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTTCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAATGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                               >  2:1816516/1‑150 (MQ=255)
                                                                        ggCAGAACGCCCTCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGCACAAGGCTTATGTGGCTGCGCGCAACAAGCTGTGCGAAGGGCGTg                                                                              >  2:1469045/1‑150 (MQ=255)
                                                                        ggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTg                                                                              <  2:1700890/150‑1 (MQ=255)
                                                                               cGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCt                                                                       <  1:1585663/150‑1 (MQ=255)
                                                                                    gcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                                  >  1:2168903/1‑150 (MQ=255)
                                                                                      gAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGc                                                                <  1:1462470/150‑1 (MQ=255)
                                                                                          ttGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTa                                                            >  2:1189542/1‑150 (MQ=255)
                                                                                             ccGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCg                                                         >  1:1124219/1‑150 (MQ=255)
                                                                                               gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTTTGCGAATGGCGTTGCAACCTGGTTAACCATAGCGac                                                       >  2:1110225/1‑149 (MQ=255)
                                                                                                 gCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagca                                                     <  1:2272269/150‑1 (MQ=255)
                                                                                                              ctctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctc                                        <  2:1732804/150‑2 (MQ=255)
                                                                                                              ctctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGctgctg                                        <  2:565547/150‑1 (MQ=255)
                                                                                                               tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctcg                                       <  1:788869/150‑3 (MQ=255)
                                                                                                               tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                       >  1:1985334/1‑150 (MQ=255)
                                                                                                               tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                       >  1:2065127/1‑150 (MQ=255)
                                                                                                               tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCATAAGCTGCTgg                                       >  1:1911028/1‑150 (MQ=255)
                                                                                                               tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTCCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                       <  1:506007/150‑1 (MQ=255)
                                                                                                                   cGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                                   >  1:520695/1‑150 (MQ=255)
                                                                                                                     aCAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGccc                                 <  1:447122/150‑1 (MQ=255)
                                                                                                                      cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAgg                                                                                                                     <  2:1298705/65‑1 (MQ=255)
                                                                                                                      cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAgg                                                                                                                     >  1:1298705/1‑65 (MQ=255)
                                                                                                                      cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCg                                >  2:1599246/1‑150 (MQ=255)
                                                                                                                          ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGcc                            >  2:587495/1‑150 (MQ=255)
                                                                                                                              tGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCa                        >  1:1799487/1‑150 (MQ=255)
                                                                                                                               gCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa                       <  1:969871/150‑1 (MQ=255)
                                                                                                                               ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                      <  1:2052042/100‑1 (MQ=39)
                                                                                                                               ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                      >  2:2052042/4‑103 (MQ=39)
                                                                                                                                 tGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAaga                     <  1:1176686/150‑1 (MQ=255)
                                                                                                                                  gTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCGCCGTGAATCGCGGTGCAGTCTGGTTGGCCGTAGAGAGCACCTTACGCTGGAGGGCGCCCGTGCCAGCAagag                    >  1:1252117/1‑150 (MQ=255)
                                                                                                                                  gTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAagag                    <  1:1110225/150‑1 (MQ=255)
                                                                                                                                   ttCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGc                   >  1:379592/1‑150 (MQ=255)
                                                                                                                                   ttCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTACGCTGCTGGGCGCCCGTGCCAGCAAGAGc                   <  1:574433/150‑1 (MQ=255)
                                                                                                                                     cATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCtt                 >  2:2038273/1‑150 (MQ=255)
                                                                                                                                      aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                               <  2:2023409/87‑1 (MQ=255)
                                                                                                                                      aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                               >  1:2023409/1‑87 (MQ=255)
                                                                                                                                      aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCGCGTGTCGGGAAGAGCTTg                >  2:1870025/1‑150 (MQ=255)
                                                                                                                                      aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTg                >  1:1214969/1‑150 (MQ=255)
                                                                                                                                       tCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaac                                                                                                 >  2:1810012/1‑68 (MQ=255)
                                                                                                                                       tCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaac                                                                                                 <  1:1810012/68‑1 (MQ=255)
                                                                                                                                         cAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGcc              >  2:1510293/1‑149 (MQ=255)
                                                                                                                                         cAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGcc              <  1:1510293/149‑1 (MQ=255)
                                                                                                                                            gACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGct          <  2:1655145/150‑2 (MQ=255)
                                                                                                                                               cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGAc       >  2:294877/1‑150 (MQ=255)
                                                                                                                                                tGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGACAGCGCGCAACAAGCTGTGCGAATGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGAcc      <  1:1592419/150‑1 (MQ=255)
                                                                                                                                                 ggACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCtacg                                               <  2:41206/108‑4 (MQ=255)
                                                                                                                                                 ggACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCtacg                                               >  1:41206/1‑105 (MQ=255)
                                                                                                                                                 ggACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACct     >  2:935471/1‑150 (MQ=255)
                                                                                                                                                  gACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACctg    >  1:1128430/1‑150 (MQ=255)
                                                                                                                                                   aCGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctc                                        <  1:754444/113‑2 (MQ=255)
                                                                                                                                                   aCGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctc                                        >  2:754444/1‑112 (MQ=255)
                                                                                                                                                    cGAACTGGGCAGCAGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACctgct  <  2:1574828/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCTGCT  >  NC_002947/1886496‑1886793

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: