Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,711,524 A→C 6.3% P6P (CCT→CCG)  gpsA ← glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,711,5240AC 6.3% 199.2 / 4.7 63P6P (CCT→CCG) gpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base C (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

ACCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711377‑4711665
                                                                                                                                                   |                                                                                                                                             
aCCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAg                                                                                                                                              >  7:96864/1‑149 (MQ=255)
  ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  6:368596/1‑149 (MQ=255)
  ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  4:84743/1‑149 (MQ=255)
  ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  2:250292/1‑149 (MQ=255)
   cTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGtt                                                                                                                                           >  2:38690/1‑149 (MQ=255)
    ttGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTg                                                                                                                                          >  6:381120/1‑149 (MQ=255)
      gagaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTg                                                                                                                                          <  3:229867/147‑1 (MQ=255)
      gagaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTg                                                                                                                                          <  3:224994/147‑1 (MQ=255)
       agaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACGGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTg                                                                                                                                       <  8:96864/149‑1 (MQ=255)
        gaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGt                                                                                                                                      <  2:301369/149‑1 (MQ=255)
         aTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGtt                                                                                                                                     <  1:250292/149‑1 (MQ=255)
             cgcgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACAGCAACCTGTTGCGGTTctgg                                                                                                                                 <  3:241696/149‑4 (MQ=255)
               cgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCa                                                                                                                               >  5:338891/1‑149 (MQ=255)
                ggATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCa                                                                                                                               >  7:165103/1‑148 (MQ=255)
                ggATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACAGCAACCTGTTGCGGTTctgggcg                                                                                                                              <  6:393854/149‑7 (MQ=255)
                      tcGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaa                                                                                                                        >  5:214055/1‑149 (MQ=255)
                        gcgATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaatc                                                                                                                      <  1:302137/149‑1 (MQ=255)
                             tCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCAtt                                                                                                                  >  5:349742/1‑148 (MQ=255)
                              cacaCGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTa                                                                                                                 <  7:24405/148‑1 (MQ=255)
                               acacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                               <  1:326790/149‑1 (MQ=255)
                                cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                               <  1:260991/148‑1 (MQ=255)
                                cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTa                                                                                                                 >  2:260991/1‑146 (MQ=255)
                                 acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaat                                                                                                                       <  6:202361/139‑1 (MQ=255)
                                 acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaat                                                                                                                       >  5:202361/1‑139 (MQ=255)
                                 acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGc                                                                                                             >  1:238677/1‑149 (MQ=255)
                                   gCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGcc                                                                                                            >  7:360339/1‑148 (MQ=255)
                                     aTGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaa                                                                                                          >  7:21867/1‑148 (MQ=255)
                                      tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                        <  5:381120/149‑1 (MQ=255)
                                        gCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAg                                                                                                       <  6:134623/148‑1 (MQ=255)
                                         ccTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt                                                                                                     <  6:164589/149‑1 (MQ=255)
                                         ccTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt                                                                                                     <  5:368596/149‑1 (MQ=255)
                                          cTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt                                                                                                     >  4:23352/1‑148 (MQ=255)
                                                        ggTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGc                                                                                      <  6:349742/149‑1 (MQ=255)
                                                          tCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCa                                                                                    <  8:114467/149‑1 (MQ=255)
                                                            gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAt                                                                                   <  1:38690/148‑1 (MQ=255)
                                                            gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAt                                                                                   >  3:285745/1‑148 (MQ=255)
                                                                    aTTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCt                                                                          >  7:362992/1‑149 (MQ=255)
                                                                        gcgcACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGcc                                                                      >  2:139526/1‑149 (MQ=255)
                                                                              cGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGCCCCGGCCATTATACGGCTCGCCTGCTa                                                                 >  7:186501/1‑148 (MQ=255)
                                                                              cGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTAc                                                                <  3:23352/149‑1 (MQ=255)
                                                                                 cacaCCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGa                                                             <  8:165103/149‑1 (MQ=255)
                                                                                    acCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAAcc                                                          <  6:43190/149‑1 (MQ=255)
                                                                                          ttCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgtgt                                                     <  5:240472/148‑1 (MQ=255)
                                                                                             cGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTTTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTg                                                  <  8:360339/148‑1 (MQ=255)
                                                                                             cGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTg                                                  <  3:84743/148‑1 (MQ=255)
                                                                                               ccAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                               >  4:321979/1‑85 (MQ=255)
                                                                                                 agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgt                                                       <  1:252029/139‑1 (MQ=255)
                                                                                                 agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgt                                                       >  2:252029/1‑139 (MQ=255)
                                                                                                 agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg                                                 <  3:251362/145‑1 (MQ=255)
                                                                                                 agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg                                                 >  4:251362/1‑145 (MQ=255)
                                                                                                     ggTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGc                                         <  6:338891/149‑1 (MQ=255)
                                                                                                       tttGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCgg                                       <  8:21867/149‑1 (MQ=255)
                                                                                                               ggcggTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGcc                               <  4:285745/149‑1 (MQ=255)
                                                                                                                     gCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGCCCCGGCCTTTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCa                         >  3:279159/1‑149 (MQ=255)
                                                                                                                                   cgccCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc           <  7:9368/149‑1 (MQ=255)
                                                                                                                                   cgccCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc           <  1:139526/149‑1 (MQ=255)
                                                                                                                                     gactgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc           >  5:219096/5‑147 (MQ=255)
                                                                                                                                     cccAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc           >  3:111350/1‑147 (MQ=255)
                                                                                                                                       cAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg                >  1:69565/1‑140 (MQ=255)
                                                                                                                                       cAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg                <  2:69565/140‑1 (MQ=255)
                                                                                                                                        tgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc           >  1:50663/2‑144 (MQ=255)
                                                                                                                                            ccGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAg  >  4:369956/1‑149 (MQ=255)
                                                                                                                                             cGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGACATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATagag  >  3:374837/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ACCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711377‑4711665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑