Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,886,644 C→T 5.4% 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 5.4% 330.7 / 5.9 111D423D (GAC→GATrmuCDNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base C (54/50);  new base T (4/2);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 114 total aligned reads displayed.

GTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCT  >  NC_002947/1886495‑1886790
                                                                                                                                                     |                                                                                                                                                  
gTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                    >  1:473247/1‑150 (MQ=255)
 tGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACg                                                                                                                                                   >  1:1935592/1‑150 (MQ=255)
      gTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTg                                                                                                                                              <  2:2152590/150‑1 (MQ=255)
      gTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTGATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTg                                                                                                                                              <  2:1688720/150‑1 (MQ=255)
       tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCCGCAAGGCCACGAAATTGCCAAGCGGCGTGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                             <  1:1486268/150‑1 (MQ=255)
       tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                             >  1:2265001/1‑150 (MQ=255)
       tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTgg                                                                                                                                             <  2:1056681/150‑1 (MQ=255)
          gCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCa                                                                                                                                          <  2:1935592/150‑1 (MQ=255)
                 cACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCtt                                                                                                                                   >  2:1632014/1‑150 (MQ=255)
                    gTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgca                                                                                                                                >  1:1149260/1‑150 (MQ=255)
                     ttGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgcag                                                                                                                               >  1:218085/1‑150 (MQ=255)
                        cTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            <  1:1173591/150‑1 (MQ=255)
                         tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAgg                                                                                                                           >  1:515933/1‑150 (MQ=255)
                         tGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAACAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCAGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAgg                                                                                                                           >  1:288695/1‑150 (MQ=255)
                         tGGCCACCTGGCGGGTTATCGACAGCCTGTGGAAGAAGGGGAGTCAAGAGTAGAACGCCCGCGAAATTGCAGCGATTGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAgg                                                                                                                           <  1:2378292/150‑1 (MQ=255)
                          ggCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          >  1:1926256/1‑150 (MQ=255)
                           gCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTg                                                                                                                         >  1:254549/1‑150 (MQ=255)
                             cACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGa                                                                                                                       <  2:754347/150‑1 (MQ=255)
                                 ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                    >  1:771416/1‑117 (MQ=255)
                                 ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                    <  2:771416/117‑1 (MQ=255)
                                  tGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAgg                                                                                                                  <  2:467106/150‑1 (MQ=255)
                                      ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTa                                                                                                              >  2:1614511/1‑150 (MQ=255)
                                       gTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTAt                                                                                                             >  1:1213023/1‑150 (MQ=255)
                                       gTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGACAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACACGGCTGAt                                                                                                             >  1:1043615/1‑150 (MQ=255)
                                         tATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGc                                                                                                           <  2:1149260/150‑1 (MQ=255)
                                          agcGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCg                                                                                                          <  2:207995/148‑1 (MQ=255)
                                           tCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCgg                                                                                                         >  2:883960/1‑150 (MQ=255)
                                           tCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATGGCCGAGCGGGACGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCCCCCCCTGTCATCAGGTGGACAAGGCTTATGCgg                                                                                                         <  1:663349/150‑1 (MQ=255)
                                             gACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGACCTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGTCa                                                                                                       <  2:157373/150‑1 (MQ=255)
                                             gACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCa                                                                                                       <  2:804934/150‑1 (MQ=255)
                                                  ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          <  1:1151473/126‑1 (MQ=255)
                                                  ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          >  2:1151473/1‑126 (MQ=255)
                                                  ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          >  2:2423090/1‑126 (MQ=255)
                                                  ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          <  1:2423090/126‑1 (MQ=255)
                                                  ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgcg                                                                                                  >  2:589886/1‑150 (MQ=255)
                                                    tgtgGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGca                                                                                                >  2:205747/1‑150 (MQ=255)
                                                        gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaacaa                                                                                            >  2:265930/1‑150 (MQ=255)
                                                        gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGATGGCGCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGcaacaa                                                                                            <  2:2157564/150‑1 (MQ=255)
                                                         aaGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAg                                                                                           <  1:1320358/150‑1 (MQ=255)
                                                          aGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGc                                                                                          >  2:2279852/1‑150 (MQ=255)
                                                          aGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGc                                                                                          <  2:2265001/150‑1 (MQ=255)
                                                           gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCt                                                                                         >  2:1764805/1‑150 (MQ=255)
                                                           gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCt                                                                                         >  1:1140294/1‑150 (MQ=255)
                                                             aGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgt                                                                                       <  2:1568040/150‑1 (MQ=255)
                                                               gAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            >  2:2151430/1‑111 (MQ=255)
                                                               gAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            <  1:2151430/111‑1 (MQ=255)
                                                                  cGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGaa                                                                                  <  2:254549/150‑1 (MQ=255)
                                                                  cGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGaa                                                                                  <  1:1764805/150‑1 (MQ=255)
                                                                   gTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAg                                                                                 <  2:515933/150‑1 (MQ=255)
                                                                    tCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAATGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAgg                                                                                >  1:99176/1‑150 (MQ=255)
                                                                      aaGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGc                                                                              >  2:245964/1‑150 (MQ=255)
                                                                        gggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                            <  1:1448902/150‑1 (MQ=255)
                                                                        gggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                            >  2:433573/1‑150 (MQ=255)
                                                                        gggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                            <  1:245964/150‑1 (MQ=255)
                                                                           cAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGc                                                                         >  1:2315793/1‑150 (MQ=255)
                                                                                 gCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   >  1:2424007/1‑150 (MQ=255)
                                                                                 gCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   >  1:269007/1‑150 (MQ=255)
                                                                                 gCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAACCAGTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   <  1:1614511/150‑1 (MQ=255)
                                                                                   ccGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                 <  1:1632014/150‑1 (MQ=255)
                                                                                       gAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGc                                                             <  2:1926256/150‑1 (MQ=255)
                                                                                        aaaTTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGcc                                                            <  2:218085/150‑1 (MQ=255)
                                                                                               cGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTg                                                                   <  2:1367487/136‑1 (MQ=255)
                                                                                                gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAAGCGGGGTAGCCGTAGCGag                                                    >  2:1131934/1‑150 (MQ=255)
                                                                                                gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGag                                                    <  2:1588468/150‑1 (MQ=255)
                                                                                                gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGag                                                    <  1:97938/150‑1 (MQ=255)
                                                                                                 aGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagc                                                   <  1:205747/150‑1 (MQ=255)
                                                                                                        cGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGcagcag                                                                                                                            <  2:407064/70‑1 (MQ=255)
                                                                                                         ggctggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAg                                           <  2:288695/150‑1 (MQ=255)
                                                                                                              gctCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgct                                      >  1:1965787/1‑150 (MQ=255)
                                                                                                               ctctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctc                                     <  1:608047/150‑2 (MQ=255)
                                                                                                               ctctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGctgctg                                     >  2:613998/1‑150 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATGGCGGCCCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                    >  1:663660/1‑150 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          >  2:2128403/1‑64 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGt                                                                                                                          <  1:2128403/64‑1 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGc                                                                                     >  2:399100/1‑101 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGc                                                                                     <  1:399100/101‑1 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                 <  2:494044/121‑1 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                 >  1:494044/1‑121 (MQ=255)
                                                                                                                tctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTgg                                    <  1:1580023/150‑1 (MQ=255)
                                                                                                                  tACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCGTAAGCTGCTGGgc                                  >  2:45453/1‑150 (MQ=255)
                                                                                                                    cGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGCGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgcgc                                <  2:1140294/150‑1 (MQ=255)
                                                                                                                       cAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGGGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCg                             >  2:2028024/1‑150 (MQ=255)
                                                                                                                          gTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                          <  2:2315793/150‑1 (MQ=255)
                                                                                                                           ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGt                            <  2:37223/147‑1 (MQ=255)
                                                                                                                           ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGt                            >  1:37223/1‑147 (MQ=255)
                                                                                                                             tgcccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAg                       >  2:148972/7‑150 (MQ=255)
                                                                                                                                ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCCGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCaa                    <  1:148972/147‑1 (MQ=255)
                                                                                                                                ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtg                                                                                        >  2:290039/4‑82 (MQ=37)
                                                                                                                                ccagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtg                                                                                        <  1:290039/79‑1 (MQ=38)
                                                                                                                                  agTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAaga                  >  2:2054073/2‑150 (MQ=255)
                                                                                                                                   gTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAACCTATGCGAAGGGCGTGGCACCCTGGTTAGCCGTCGCGAGCGGCTTAAGCTGCTGGGCGCCCGTGCCAGCAagag                 <  1:1994878/150‑1 (MQ=255)
                                                                                                                                    ttCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGc                >  2:1385701/1‑150 (MQ=255)
                                                                                                                                       aTCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTg             <  1:589886/150‑1 (MQ=255)
                                                                                                                                        tCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGc            >  2:314610/1‑150 (MQ=255)
                                                                                                                                        tCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGc            >  2:1999363/1‑150 (MQ=255)
                                                                                                                                                cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACaa                                                                                                                    <  1:1074907/38‑1 (MQ=255)
                                                                                                                                                cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACaa                                                                                                                    >  2:1074907/1‑38 (MQ=255)
                                                                                                                                                cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGaa                                                                                  <  1:2115627/72‑1 (MQ=255)
                                                                                                                                                cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGaa                                                                                  >  2:2115627/1‑72 (MQ=255)
                                                                                                                                                  ggACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACct  >  2:101595/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCT  >  NC_002947/1886495‑1886790

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