Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,308,076 T→C 5.5% R80R (CGT→CGC csiD → carbon starvation induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,308,0760TC 5.5% 235.9 / 2.0 82R80R (CGT→CGCcsiDcarbon starvation induced protein
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base C (2/3);  total (39/43)
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.67e-01

CTGCTGGAGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCAT  >  NC_002947/3307930‑3308225
                                                                                                                                                   |                                                                                                                                                     
ctgctgGAGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCg                                                                                                                                                    <  2:1022742/150‑1 (MQ=255)
  gctgGAGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGt                                                                                                                                                  <  2:1377216/150‑1 (MQ=255)
   ctgGAGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGtg                                                                                                                                                 <  1:1812090/150‑1 (MQ=255)
       aGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGtgctgc                                                                                                                                             <  2:131289/150‑1 (MQ=255)
              cTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAacac                                                                                                                                      >  2:1536602/1‑150 (MQ=255)
                     cggcacaCCGTCTAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTc                                                                                                                               <  2:649941/144‑1 (MQ=255)
                            ccGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACcgcg                                                                                                                        >  2:699551/1‑150 (MQ=255)
                             cGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCgcgc                                                                                                                       >  2:1151726/1‑150 (MQ=255)
                              gTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGcc                                                                                                                      >  2:526023/1‑150 (MQ=255)
                                 gAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCcaca                                                                                                                   >  2:1464953/1‑150 (MQ=255)
                                  aaGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGt                                                                                                                                                  >  1:559533/1‑118 (MQ=255)
                                  aaGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGt                                                                                                                                                  <  2:559533/118‑1 (MQ=255)
                                      gTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGtgctg                                                                                                                                              >  1:394464/1‑118 (MQ=255)
                                      gTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGtgctg                                                                                                                                              <  2:394464/118‑1 (MQ=255)
                                           tGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGTTGCTGCTCAACGCCATCCTCGACCGCGCCACAGGCGGCATGc                                                                                                         <  2:265119/150‑1 (MQ=255)
                                           tGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGc                                                                                                         >  1:785354/1‑150 (MQ=255)
                                            gCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCt                                                                                                        >  1:1327168/1‑150 (MQ=255)
                                              aGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGa                                                                                                      >  1:2129336/1‑150 (MQ=255)
                                                gCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGAtt                                                                                                    >  2:236843/1‑150 (MQ=255)
                                                  cGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTAc                                                                                                  >  2:1309289/1‑150 (MQ=255)
                                                    ttGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTAccc                                                                                                <  1:1151726/150‑1 (MQ=255)
                                                       cccAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTAACCCGa                                                                                             >  1:1560823/1‑150 (MQ=255)
                                                          aaTGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGCGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCg                                                                                          >  1:750993/1‑150 (MQ=255)
                                                              gCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCt                                                                                      <  1:1922390/150‑1 (MQ=255)
                                                                 ggTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGa                                                                                   <  1:778581/150‑1 (MQ=255)
                                                                      aGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACg                                                                              >  1:596098/1‑150 (MQ=255)
                                                                          aTTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAg                                                                          >  1:876753/1‑150 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATc                                                                                                                                  <  2:1390678/90‑1 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATc                                                                                                                                  >  1:1390678/1‑90 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAggcggc                                                                                                             <  2:756222/111‑1 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAggcggc                                                                                                             >  1:756222/1‑111 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAg                                                                      <  1:1464953/150‑1 (MQ=255)
                                                                              gAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGACAg                                                                      <  2:1162992/150‑1 (MQ=255)
                                                                                gTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGc                                                                    <  1:699551/150‑1 (MQ=255)
                                                                                gTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGc                                                                    <  2:596098/150‑1 (MQ=255)
                                                                                gTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGc                                                                    <  1:218876/150‑1 (MQ=255)
                                                                                   cAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGa                                                                 >  1:472376/1‑150 (MQ=255)
                                                                                      gTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGaa                                                                <  2:497031/148‑1 (MQ=255)
                                                                                      gTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGaa                                                                >  1:497031/1‑148 (MQ=255)
                                                                                      gTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGa                                                              >  2:1030887/1‑150 (MQ=255)
                                                                                        cGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGGACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGCTCACGTCATCCAGGCCGatcggg                                                            <  1:1309289/150‑6 (MQ=255)
                                                                                              tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGa                                                                           <  2:141813/129‑1 (MQ=255)
                                                                                              tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGa                                                                           >  1:141813/1‑129 (MQ=255)
                                                                                              tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATg                                                          <  1:624129/146‑1 (MQ=255)
                                                                                              tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCa                                                      >  1:273099/1‑150 (MQ=255)
                                                                                              tGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTCCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGGGAGCCAGGCCGAAGACATGGTCa                                                      >  1:1901103/1‑150 (MQ=255)
                                                                                              tGCGTTGTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATg                                                          >  2:624129/1‑146 (MQ=255)
                                                                                                cGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAg                                                    <  1:1242543/150‑1 (MQ=255)
                                                                                                 gTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGt                                                   >  1:1567193/1‑150 (MQ=255)
                                                                                                    ttCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTca                                                <  2:1751017/150‑1 (MQ=255)
                                                                                                     tCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGCAGACATGGTAAAGTTcac                                               <  2:1560823/150‑1 (MQ=255)
                                                                                                       gtgtAGGCGAAATCCTCACTATTTTGTGCCAGTGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTcacca                                             <  1:121345/150‑1 (MQ=255)
                                                                                                        tgtAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTcaccac                                            >  1:1414627/1‑150 (MQ=255)
                                                                                                        tgtAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTcaccac                                            >  2:1759833/1‑150 (MQ=255)
                                                                                                             gCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCt                                       <  2:472376/150‑1 (MQ=255)
                                                                                                             gCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCt                                       <  2:876753/150‑1 (MQ=255)
                                                                                                             gCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCt                                       <  1:1759833/150‑1 (MQ=255)
                                                                                                                    ccTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACANCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAGGTTCACCACCGCCTGCGCTCa                                >  1:822692/1‑150 (MQ=255)
                                                                                                                     cTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCAc                               <  1:526023/150‑1 (MQ=255)
                                                                                                                      ttGACGAGCTGTGCCAGGGCACGCTGCGCCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCt                                                                                      <  1:1332744/94‑1 (MQ=255)
                                                                                                                      ttGACGAGCTGTGCCAGGGCACGCTGCGCCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCt                                                                                      >  2:1332744/1‑94 (MQ=255)
                                                                                                                      ttGACGAGCTGTGCCAGGGCACGCTGCGCCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCt                                                                                      <  1:1786381/94‑1 (MQ=255)
                                                                                                                      ttGACGAGCTGTGCCAGGGCACGCTGCGCCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCt                                                                                      >  2:1786381/1‑94 (MQ=255)
                                                                                                                            aGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCg                        >  2:330648/1‑150 (MQ=255)
                                                                                                                              gCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCgg                       <  2:1414627/150‑1 (MQ=255)
                                                                                                                               cTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGc                      <  1:1536602/150‑1 (MQ=255)
                                                                                                                                tgcgCCAGGGCACGCTGCGTCCGGTGCCGCTCAACACCATCCTCGGCCGCGCCACAGGCGGTATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGcc                     <  1:495231/147‑1 (MQ=255)
                                                                                                                                  tgCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGt                   >  1:1454786/1‑150 (MQ=255)
                                                                                                                                  tgCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGt                   >  1:1588911/1‑150 (MQ=255)
                                                                                                                                  tgCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGt                   >  1:2022558/1‑150 (MQ=255)
                                                                                                                                   gCCAGGGCACGCTGCGCCCGGTGCTGCTCAACACCATCCTCCACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGtt                  <  2:1901103/150‑1 (MQ=255)
                                                                                                                                    ccAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTc                 <  2:1454786/150‑1 (MQ=255)
                                                                                                                                      aGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCa               >  2:2324082/1‑150 (MQ=255)
                                                                                                                                       gggCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCaa              <  1:1030887/150‑1 (MQ=255)
                                                                                                                                       gggCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCaa              <  2:1588911/150‑1 (MQ=255)
                                                                                                                                           ccGCTGCGTCCGGTGCTGCTCAACACCATCCTCGAGCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCTTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTAt          <  2:822692/149‑1 (MQ=255)
                                                                                                                                             gCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATc                         <  1:1703060/133‑1 (MQ=255)
                                                                                                                                             gCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATc                         >  2:1703060/1‑133 (MQ=255)
                                                                                                                                               tGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCGTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACGGCCTGCGCTCACCTGATGGGCCGTTCCAACTGTGACg      >  1:829784/1‑150 (MQ=255)
                                                                                                                                                 cGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGcc    <  2:2129336/150‑1 (MQ=255)
                                                                                                                                                  gTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACa                                                            <  1:393234/93‑1 (MQ=255)
                                                                                                                                                  gTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACa                                                            >  2:393234/1‑93 (MQ=255)
                                                                                                                                                   tCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCACCTATGACGCCAt  >  1:544633/1‑150 (MQ=255)
                                                                                                                                                   tCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCAt  <  1:330648/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGCTGGAGCTGACCTTTGCCCGGGAAACCGTCGAAGCGTTCGTGCAGGCCGTTGCCCAATGGCCGGTGCAGGCATTGGAGTACAAGTCGTTCCTGCGTTTTCGTGTAGGCGAAATCCTTGACGAGCTGTGCCAGGGCACGCTGCGTCCGGTGCTGCTCAACACCATCCTCGACCGCGCCACAGGCGGCATGCTGATTACCCCGATCGGCCTGGATGACGTGAGCCAGGCCGAAGACATGGTCAAGTTCACCACCGCCTGCGCTCACCTGATCGGCCGTTCCAACTATGACGCCAT  >  NC_002947/3307930‑3308225

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: