Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,604,251 C→G 5.7% P12R (CCG→CGG) ‡ PP_4076 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,604,2510CG 5.7% 338.6 / 8.4 105P12R (CCG→CGG) ‡PP_4076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (48/50);  new base G (3/3);  total (51/54)
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.81e-01

Only 100 of 107 total aligned reads displayed.

TCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604104‑4604396
                                                                                                                                                   |                                                                                                                                                  
tcGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                  >  1:2124038/1‑150 (MQ=255)
tcGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                  >  1:1828072/1‑150 (MQ=255)
tcGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                  >  1:1208402/1‑150 (MQ=255)
  ggACGTGGCGGGTGATCCGGGCTGGCACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTcc                                                                                                                                                >  2:685362/1‑150 (MQ=255)
        ggCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATc                                                                                                                                          <  2:2124038/150‑1 (MQ=255)
           gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCgg                                                                                                                                       >  2:2322921/1‑150 (MQ=255)
           gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGACCNCGCCCTTCAACTACATAAGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCTGACCCGTGCGCCATCCgg                                                                                                                                       >  1:2773917/1‑150 (MQ=255)
             gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTGCGCCATCCGGGa                                                                                                                                     >  1:1602689/1‑150 (MQ=255)
             gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                     <  2:1453775/150‑1 (MQ=255)
              tgctCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCGTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAg                                                                                                                                    <  1:1150525/147‑1 (MQ=255)
               gATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAgg                                                                                                                                   <  1:11654/150‑1 (MQ=255)
                 tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGccctaaag                                                                                                                                       >  1:1667246/1‑138 (MQ=255)
                 tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaag                                                                                                                                       <  2:1667246/144‑4 (MQ=255)
                 tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                     <  1:1384454/146‑1 (MQ=255)
                   cGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtc                                                                                                                               >  1:1718786/1‑150 (MQ=255)
                   cGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtc                                                                                                                               >  2:1310540/1‑150 (MQ=255)
                       cTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGc                                                                                                                          <  2:1208402/150‑1 (MQ=255)
                       cTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGc                                                                                                                          <  2:1187105/150‑1 (MQ=255)
                          gTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGCGCCGGTCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                       <  1:2385415/150‑1 (MQ=255)
                            aCCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                     <  1:685362/150‑1 (MQ=255)
                             ccAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaa                                                                                                                    >  1:2790540/1‑150 (MQ=255)
                             ccAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaa                                                                                                                    >  1:2907436/1‑150 (MQ=255)
                              cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                   >  2:787163/1‑150 (MQ=255)
                               aGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAg                                                                                                                  <  2:1344335/150‑1 (MQ=255)
                                    tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCtt                                                                                                             <  2:1718786/150‑1 (MQ=255)
                                    tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCtt                                                                                                             >  1:2869682/1‑150 (MQ=255)
                                      gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaa                                                                                                                                        <  2:277720/122‑3 (MQ=255)
                                      gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGGGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                           <  2:1909989/150‑1 (MQ=255)
                                       ctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGc                                                                                                          >  2:928868/1‑150 (MQ=255)
                                        tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                       >  2:516582/1‑136 (MQ=255)
                                        tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                       <  1:516582/136‑1 (MQ=255)
                                        tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTCGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGca                                                                                                         >  2:326567/1‑150 (MQ=255)
                                               gcgcCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTc                                                                                                  <  1:1506342/150‑1 (MQ=255)
                                                  cccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                               <  1:928868/150‑1 (MQ=255)
                                                  cccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                               >  1:35043/1‑150 (MQ=255)
                                                      cACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCa                                                                                           <  1:787163/150‑1 (MQ=255)
                                                       aCTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCAt                                                                                          >  1:1110227/1‑150 (MQ=255)
                                                        cTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATc                                                                                         >  2:2169404/1‑150 (MQ=255)
                                                         ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCAAGGAAAGCGCCTCGCACAAGTTCGCGCCCATCa                                                                                        <  2:2773917/150‑1 (MQ=255)
                                                          tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTc                                                                                                  <  1:2595703/139‑1 (MQ=255)
                                                          tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTc                                                                                                  >  2:2595703/1‑139 (MQ=255)
                                                          tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCaa                                                                                       >  1:1958538/1‑150 (MQ=255)
                                                           cAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAg                                                                                      >  1:904353/1‑150 (MQ=255)
                                                            aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                     <  1:39389/150‑1 (MQ=255)
                                                                   tGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCGATCCGGGAGGCGTCTCGCCAGGAAAgcgc                                                                                                               <  1:2862023/117‑1 (MQ=255)
                                                                   tGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCACCATCCGGGAGGCGTCTCGCCAGGAAAgcgc                                                                                                               >  2:2862023/1‑117 (MQ=255)
                                                                    gcgcgcGAACTTCGATATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGCTCGCTTGATGGCGCCCCCGTCCGCCATCCGGGAGACGTCTCCCCAGGAAAGCGCGACGCTCGGGTTCGCGGCGATCAAGTGCCACCTc                                                                             <  2:811936/148‑1 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                   <  1:1310540/150‑1 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                   >  1:2402383/1‑150 (MQ=255)
                                                                               tCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCt                                                                  <  2:1958538/150‑1 (MQ=255)
                                                                                cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTg                                                                 <  2:2869682/150‑1 (MQ=255)
                                                                                 gAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGa                                                                >  2:544968/1‑150 (MQ=255)
                                                                                  aaaTTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGCATGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGaa                                                               >  1:2202054/1‑150 (MQ=255)
                                                                                       tCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTc                                                          <  1:2169404/150‑1 (MQ=255)
                                                                                            ggTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCaa                                                     <  1:454383/150‑1 (MQ=255)
                                                                                             gTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCAt                                                                                                                                           >  2:3070863/1‑64 (MQ=255)
                                                                                             gTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCAt                                                                                                                                           <  1:3070863/64‑1 (MQ=255)
                                                                                                       tCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCa                                          <  2:2907436/150‑1 (MQ=255)
                                                                                                         gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCaga                                        <  2:2402383/150‑1 (MQ=255)
                                                                                                         gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCaga                                        <  2:2499910/150‑1 (MQ=255)
                                                                                                           cacGCAATGAGTGGTTTCTTCGATCGCTTGATAGCGAACCAGGCGGTGGTCCGGTGGACGGGGCTGCAGGAAGGCGCTTCGCACAAGTTCGCGGCCATAACGTGACACATCGAACAGTTGAAGGCCACTCTCAAAGCCTGCACCCAGagc                                      <  2:2202054/150‑1 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                     >  1:2008442/1‑150 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                     >  2:2560076/1‑150 (MQ=255)
                                                                                                               cAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGcc                                  >  2:475877/1‑150 (MQ=255)
                                                                                                                aaTGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                   >  1:2920278/1‑100 (MQ=255)
                                                                                                                aaTGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                   <  2:2920278/100‑1 (MQ=255)
                                                                                                                     gTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGc                            >  2:2752747/1‑150 (MQ=255)
                                                                                                                      tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaaa                                                                         <  2:2795342/104‑1 (MQ=255)
                                                                                                                      tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaaa                                                                         >  1:2795342/1‑104 (MQ=255)
                                                                                                                      tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCaa                                                     <  2:1907857/124‑1 (MQ=255)
                                                                                                                      tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCaa                                                     >  1:1907857/1‑124 (MQ=255)
                                                                                                                      tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGATCGTCGTGAAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCt                           >  2:2574455/1‑150 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATc                                                                                         >  1:135394/1‑86 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATc                                                                                         <  2:135394/86‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAg                                                    >  2:1787498/1‑123 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAg                                                    <  1:1787498/123‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                         >  2:555281/1‑150 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTCGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                         <  1:326567/150‑1 (MQ=255)
                                                                                                                         tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCa                                                                                           >  2:72427/1‑83 (MQ=255)
                                                                                                                         tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCa                                                                                           <  1:72427/83‑1 (MQ=255)
                                                                                                                         tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCAGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTggg                        >  1:2484587/1‑150 (MQ=255)
                                                                                                                            cttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCt                     <  1:555281/150‑1 (MQ=255)
                                                                                                                             ttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGtt                                                                     <  2:43212/101‑1 (MQ=255)
                                                                                                                             ttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGtt                                                                     >  1:43212/1‑101 (MQ=255)
                                                                                                                               cGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCACCGAGCTGGGCCTgcg                  >  2:2478745/1‑150 (MQ=255)
                                                                                                                                 aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGcgcg                >  1:995562/1‑150 (MQ=255)
                                                                                                                                 aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGcgcg                >  2:781087/1‑150 (MQ=255)
                                                                                                                                    gCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACt             <  1:199995/150‑1 (MQ=255)
                                                                                                                                     ccggATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                    <  1:1109584/91‑1 (MQ=255)
                                                                                                                                     ccggATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                    >  2:1109584/4‑94 (MQ=255)
                                                                                                                                      ttGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc           >  2:1661401/1‑150 (MQ=255)
                                                                                                                                          tGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgcgc                 <  2:1184607/140‑1 (MQ=255)
                                                                                                                                          tGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgcgc                 >  1:1184607/1‑140 (MQ=255)
                                                                                                                                           ggCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACgg      <  1:2478745/150‑1 (MQ=255)
                                                                                                                                             cGGACCCGGCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGc    <  2:1602689/150‑1 (MQ=255)
                                                                                                                                              ggACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGTTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCa   >  1:2971015/1‑150 (MQ=255)
                                                                                                                                              ggACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCa   <  1:781087/150‑1 (MQ=255)
                                                                                                                                               gACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt  <  2:904353/150‑1 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGt                                                                      >  1:1894121/1‑79 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGt                                                                      <  2:1894121/79‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604104‑4604396

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: