Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,994 T→A 10.1% intergenic (+35/‑105) adk → / → tsaB adenylate kinase/endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,9940TA10.1% 346.6 / 22.8 117intergenic (+35/‑105)adk/tsaBadenylate kinase/endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (48/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 = 9.27e-01

Only 100 of 119 total aligned reads displayed.

GGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712845‑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)
    ttCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCgg                                                                                                                                                <  1:154902/150‑1 (MQ=255)
      cTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGc                                                                                                                                              <  1:5055408/150‑1 (MQ=255)
      cTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGc                                                                                                                                              >  2:1785901/1‑150 (MQ=255)
       tACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGcc                                                                                                                                             >  2:3397407/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                                                                                                                                          <  1:4808979/150‑1 (MQ=255)
          cAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                          <  2:4308146/150‑1 (MQ=255)
            gAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGt                                                                                                                                        >  1:5112442/1‑150 (MQ=255)
            gAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGATGt                                                                                                                                        >  1:68132/1‑150 (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)
                       cTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGc                                                                                                                             <  1:865683/150‑1 (MQ=255)
                             aTGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                       <  2:4200547/150‑1 (MQ=255)
                             aTGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGcaca                                                                                                                       <  1:253728/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                                                                                                          <  2:1337651/150‑1 (MQ=255)
                                          gAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAgg                                                                                                          >  1:2394335/1‑150 (MQ=255)
                                           aaGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGa                                                                                                         <  2:1985941/150‑1 (MQ=255)
                                           aaGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGa                                                                                                         <  1:404578/150‑1 (MQ=255)
                                            aGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGac                                                                                                        <  2:3095804/150‑1 (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)
                                                        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                                                                                                                                          <  1:3724890/95‑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)
                                                                            gAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGctct                                                                        <  1:3653526/150‑1 (MQ=255)
                                                                              aGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCACAGGACACCACCCCGCGGCTGCTCGATACGGCCGCTCttt                                                                      >  1:4857072/1‑150 (MQ=255)
                                                                               gCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttt                                                                     <  2:4942573/150‑1 (MQ=255)
                                                                                ccATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCttttt                                                                    <  1:2258522/150‑1 (MQ=255)
                                                                                 cATCACCGCCAAGGTCCTGGCAGCCCTGAGCGGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTc                                                                   <  1:2643432/150‑1 (MQ=255)
                                                                                  aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGGGTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttttaa                                                                  >  2:1650192/1‑148 (MQ=255)
                                                                                     aCCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTc                                                               <  1:4645609/150‑1 (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)
                                                                                                  tGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAATCGCACCTGGATAcc                                                  >  2:3684980/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)
                                                                                                           tGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTTGGCCGCACAGCGCCCCCACAGGACACAACCCCGCGGCTGCGCGATACAGCCGCTCGTTTTGATTGGGACCGGGATAACCCTTCGATg                                         >  1:1384545/1‑150 (MQ=255)
                                                                                                              gCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacc                                      >  1:764420/1‑150 (MQ=255)
                                                                                                               cTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacca                                     <  2:3719579/150‑1 (MQ=255)
                                                                                                                tGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgc                                                                                                 >  1:3341881/1‑89 (MQ=255)
                                                                                                                tGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAAcgc                                                                                                 <  2:3341881/89‑1 (MQ=255)
                                                                                                                     cGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCctgc                               <  1:697766/150‑1 (MQ=255)
                                                                                                                             gACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTg                       <  1:3556114/150‑1 (MQ=255)
                                                                                                                              aCGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTgg                      <  1:2422151/150‑1 (MQ=255)
                                                                                                                                   ccTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATAcc                 <  2:1452900/150‑1 (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)
                                                                                                                                          aaCGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCg          >  2:2665105/1‑150 (MQ=255)
                                                                                                                                          aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                     >  1:1671271/1‑139 (MQ=255)
                                                                                                                                          aaCGGCCCGCAAGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGa                     <  2:1671271/139‑1 (MQ=255)
                                                                                                                                             ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAg       >  1:1963681/1‑150 (MQ=255)
                                                                                                                                             ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTAGCACCAGGATACCCGTTAGATGAACACCCAGCTGGCCCAGGATACCGCCACCCAAg       >  1:1793115/1‑150 (MQ=255)
                                                                                                                                                ccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCt    >  2:4751260/1‑150 (MQ=255)
                                                                                                                                                  gCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGt  <  1:2171292/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGATTTCTACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712845‑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: