Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,993 T→A 10.2% intergenic (+34/‑106) adk → / → tsaB adenylate kinase/endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,9930TA10.2% 313.2 / 22.9 117intergenic (+34/‑106)adk/tsaBadenylate kinase/endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (46/57);  new base A (5/7);  total (53/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 120 total aligned reads displayed.

TGGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712844‑1713140
                                                                                                                                                     |                                                                                                                                                    
 ggATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCtt                                                                                                                                                     <  2:104362/150‑1 (MQ=255)
 ggATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCtt                                                                                                                                                     <  2:2841541/150‑1 (MQ=255)
   aTTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGc                                                                                                                                                   <  1:101199/150‑1 (MQ=255)
    tttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCg                                                                                                                                                  >  2:3437879/1‑150 (MQ=255)
    tttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCAACAACGGCCCGCTTGCg                                                                                                                                                  >  2:4819495/1‑150 (MQ=255)
     ttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCgg                                                                                                                                                 <  1:322774/150‑1 (MQ=255)
       cTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGc                                                                                                                                               <  1:5055408/150‑1 (MQ=255)
        tACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGcc                                                                                                                                              >  2:3397407/1‑150 (MQ=255)
         aCCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCg                                                                                                                                             >  1:1964578/1‑150 (MQ=255)
          ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGGCGt                                                                                                                                            >  2:4246964/1‑150 (MQ=255)
          ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGt                                                                                                                                            <  2:4487404/150‑1 (MQ=255)
          ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCGTGCGGGCCGt                                                                                                                                            >  2:3779581/1‑150 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                           <  1:2301815/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                           <  1:2760886/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                           <  2:5038651/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                           <  2:1975169/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                           <  2:4308146/150‑1 (MQ=255)
           cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                           <  1:4808979/150‑1 (MQ=255)
             gAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGt                                                                                                                                         >  1:5112442/1‑150 (MQ=255)
              aaGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTc                                                                                                                                        <  2:1711931/150‑1 (MQ=255)
               aGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTCg                                                                                                                                       <  2:1062751/150‑1 (MQ=255)
                 cTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGtt                                                                                                                                     <  2:1801044/150‑1 (MQ=255)
                     cGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTa                                                                                                                                 <  1:4224798/150‑1 (MQ=255)
                      ggCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGTTCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAg                                                                                                                                <  2:4332008/150‑1 (MQ=255)
                         tGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGcc                                                                                                                             <  2:1272730/150‑1 (MQ=255)
                              aTGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                        <  1:253728/150‑1 (MQ=255)
                              aTGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                        <  2:4200547/150‑1 (MQ=255)
                              aTGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                        >  1:836862/1‑150 (MQ=255)
                                  cggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgcgc                                                                                                                    <  1:2017405/150‑1 (MQ=255)
                                  cggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgcgc                                                                                                                    <  1:4151363/150‑1 (MQ=255)
                                          cGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAg                                                                                                            <  2:864293/150‑1 (MQ=255)
                                           gAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAgg                                                                                                           >  1:2394335/1‑150 (MQ=255)
                                           gAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAgg                                                                                                           <  2:1337651/150‑1 (MQ=255)
                                            aaGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGa                                                                                                          <  1:404578/150‑1 (MQ=255)
                                            aaGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGa                                                                                                          <  2:1985941/150‑1 (MQ=255)
                                             aGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGac                                                                                                         <  2:3095804/150‑1 (MQ=255)
                                                aCAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGAcaca                                                                                                      >  2:3272944/1‑150 (MQ=255)
                                                  aGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACCACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAc                                                                                                    <  2:4435775/150‑1 (MQ=255)
                                                  aGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAc                                                                                                    <  2:1052528/150‑1 (MQ=255)
                                                    ccACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgc                                                                                                  >  1:1124678/1‑150 (MQ=255)
                                                     cacaTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgcc                                                                                                 >  1:1231481/1‑150 (MQ=255)
                                                      acaTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgccg                                                                                                >  2:3644027/1‑150 (MQ=255)
                                                         tCGAAGGTGTCGGCTCGGTTGAAGTCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctg                                                                                             <  1:1406011/150‑1 (MQ=255)
                                                         tCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctg                                                                                             <  1:3669791/150‑1 (MQ=255)
                                                         tCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACAACGCCgctg                                                                                             >  2:3545116/1‑150 (MQ=255)
                                                          cGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgc                                                                                            <  1:2762349/150‑1 (MQ=255)
                                                           gAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACCGCCCGCGTGCGGGCCGTTTTGGTTTTTAGGCCGCACAGCCCCCCCACAGGACACAACCCCcctcct                                                                                           >  1:1460809/1‑150 (MQ=255)
                                                               gtgtCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGctc                                                                                       <  1:1734868/150‑1 (MQ=255)
                                                               gtgtCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGctc                                                                                       <  2:2099642/150‑1 (MQ=255)
                                                                  tCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCACGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                           >  2:3724890/1‑95 (MQ=255)
                                                                    ggCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATAc                                                                                  >  2:4899978/1‑150 (MQ=255)
                                                                     gCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCtataat                                                                                 >  2:847038/1‑148 (MQ=255)
                                                                       tCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGc                                                                               >  1:2396983/1‑150 (MQ=255)
                                                                       tCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACAACGCCGCTGCTGCTCTATACTGc                                                                               >  1:3272405/1‑150 (MQ=255)
                                                                       tCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACAACGCCGCTGCTGCTCGATACGGc                                                                               >  1:1147778/1‑150 (MQ=255)
                                                                          gTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGc                                                                            <  1:1785901/150‑1 (MQ=255)
                                                                            tGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctc                                                                          <  1:249035/150‑1 (MQ=255)
                                                                             gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctct                                                                         <  1:3653526/150‑1 (MQ=255)
                                                                               aGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACCACCCCGCGGCTGCTCGATACGGCCGCTCttt                                                                       >  1:4857072/1‑150 (MQ=255)
                                                                                gCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttt                                                                      <  2:4942573/150‑1 (MQ=255)
                                                                                   aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGGTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttttaa                                                                   >  2:1650192/1‑148 (MQ=255)
                                                                                    tCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGGCGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAt                                                                  >  1:836332/1‑150 (MQ=255)
                                                                                      aCCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTc                                                                >  2:4079768/1‑150 (MQ=255)
                                                                                         gCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCa                                                             <  1:3784765/150‑1 (MQ=255)
                                                                                         gCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCa                                                             <  1:2306560/150‑1 (MQ=255)
                                                                                            aaGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTAGCACCt                                                          >  1:2461911/1‑150 (MQ=255)
                                                                                             aGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTg                                                         <  1:2077673/150‑1 (MQ=255)
                                                                                              ggTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTgg                                                        <  2:1674568/150‑1 (MQ=255)
                                                                                                tCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCGATACTGCCGCTCTTTTTAAATCGCCCCAGGAt                                                      >  1:3584370/1‑150 (MQ=255)
                                                                                                 ccTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATa                                                    >  2:351878/1‑150 (MQ=255)
                                                                                                 ccTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATa                                                    <  1:70176/150‑1 (MQ=255)
                                                                                                    ggCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCCCAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATAccc                                                 >  1:1596236/1‑150 (MQ=255)
                                                                                                    ggCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATAccc                                                 <  1:131048/150‑1 (MQ=255)
                                                                                                       aGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGtt                                              <  1:2450678/150‑1 (MQ=255)
                                                                                                       aGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTACACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGtt                                              <  1:1525748/150‑1 (MQ=255)
                                                                                                         cccTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCg                                            >  2:4062884/1‑150 (MQ=255)
                                                                                                          ccTGAGCTGATCCGGCGCCTGGCGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGa                                           <  2:836332/150‑1 (MQ=255)
                                                                                                           cTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGAt                                          <  2:1657847/150‑1 (MQ=255)
                                                                                                            tGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTTGGCCGCCCAGCGCCCCCGCAGGACACAACCCCGCTGCTGCGCTATACTGGCGCGCTTTTTAAATAGCACCAGGATAACCCGtagatg                                          >  1:1384605/1‑145 (MQ=255)
                                                                                                               gCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacc                                      >  1:764420/1‑150 (MQ=255)
                                                                                                                cTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacca                                     <  2:3719579/150‑1 (MQ=255)
                                                                                                                      cGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctgc                               <  2:2396983/150‑1 (MQ=255)
                                                                                                                      cGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctgc                               <  1:697766/150‑1 (MQ=255)
                                                                                                                              gACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTg                       <  1:3556114/150‑1 (MQ=255)
                                                                                                                                cGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                     <  1:4195866/150‑1 (MQ=255)
                                                                                                                                    aaacGACAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTGATTAGCACCTGGATACCCGTTAGCTGACCACCCTGCAGGGCCGGGATAcc                 >  1:2393617/4‑150 (MQ=255)
                                                                                                                                    aaacGACAACGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATAcc                 >  1:249953/4‑150 (MQ=255)
                                                                                                                                     cTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCg                >  2:3194398/1‑150 (MQ=255)
                                                                                                                                      tCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGc               >  2:3251640/1‑150 (MQ=255)
                                                                                                                                        gACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTATGCCGCACAGCGCCCCCTCAGGACACAACGCCGCAGCAGCTCTATACTGCCCCTCTTTTTAATTAGCACCTGGATACCCCGTCGATGACCACCCAGCTGGCCCTGGATACCGCCa             >  1:1473931/1‑150 (MQ=255)
                                                                                                                                           aaCGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCg          >  2:64635/1‑150 (MQ=255)
                                                                                                                                           aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                     <  2:1671271/139‑1 (MQ=255)
                                                                                                                                           aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                     >  1:1671271/1‑139 (MQ=255)
                                                                                                                                           aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCg          >  1:5070507/1‑150 (MQ=255)
                                                                                                                                              ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAg       >  1:1963681/1‑150 (MQ=255)
                                                                                                                                              ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTAGCACCTGGATAACCGTTCGATGACCACCCTGCTGGCCCGGGATACCGCCACCGaac       >  1:2798090/1‑149 (MQ=255)
                                                                                                                                              ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTAGCACCAGGATACCCGTTAGATGAACACCCAGCTGGCCCAGGATACCGCCACCCAAg       >  1:1793115/1‑150 (MQ=255)
                                                                                                                                               gCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGc      <  1:3644027/150‑1 (MQ=255)
                                                                                                                                                 ccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCt    >  2:4751260/1‑150 (MQ=255)
                                                                                                                                                   gCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGt  <  1:2171292/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712844‑1713140

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