Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,711,531 T→G 6.3% Q4P (CAG→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,5310TG 6.3% 199.6 / 4.5 63Q4P (CAG→CCG) gpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (23/36);  new base G (2/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCAC  >  NC_002947/4711385‑4711668
                                                                                                                                                   |                                                                                                                                         
gaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGt                                                                                                                                         >  7:102119/1‑149 (MQ=255)
 aTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGTCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGtt                                                                                                                                        <  3:241323/149‑1 (MQ=255)
    gcgcGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAg                                                                                                                                     >  5:400119/1‑149 (MQ=255)
      gcgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTc                                                                                                                                   <  5:240465/149‑1 (MQ=255)
      gcgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTc                                                                                                                                   <  6:126910/149‑1 (MQ=255)
      gcgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACAGCAACCTGTTGCGGTTctggg                                                                                                                                   <  7:317807/149‑5 (MQ=255)
       cgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCa                                                                                                                                  >  4:316137/1‑149 (MQ=255)
              tcGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaa                                                                                                                           >  7:21928/1‑149 (MQ=255)
                 cgATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaatc                                                                                                                         >  7:36423/1‑148 (MQ=255)
                       acacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                                  <  8:102119/149‑1 (MQ=255)
                         acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGc                                                                                                                <  7:294927/149‑1 (MQ=255)
                           gCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGcc                                                                                                               >  8:306007/1‑148 (MQ=255)
                              tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                           <  8:21928/149‑1 (MQ=255)
                              tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                           <  7:94923/149‑1 (MQ=255)
                                gCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGt                                                                                                         >  7:84581/1‑149 (MQ=255)
                                 ccTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGt                                                                                                         >  2:177601/1‑148 (MQ=255)
                                  cTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACCCCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCAAAGCCATTAAAGCCAAAAttt                                                                                                        >  6:289543/1‑148 (MQ=255)
                                  cTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCTTCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt                                                                                                        <  1:112283/149‑1 (MQ=255)
                                   tCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTaa                                                                                                      >  6:301959/1‑149 (MQ=255)
                                     ggCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACAGCGTTTTCCGCCAGCATGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTaa                                                                                                      <  3:70470/147‑1 (MQ=255)
                                     ggCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTaa                                                                                                      <  8:117523/147‑1 (MQ=255)
                                            gcCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGAcc                                                                                             >  7:225922/1‑149 (MQ=255)
                                             ccGGGTCGCGCATGCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGAccc                                                                                            <  3:280586/149‑1 (MQ=255)
                                              cGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCg                                                                                           <  4:46026/149‑1 (MQ=255)
                                                 gTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGAGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGcc                                                                                        <  8:84581/149‑1 (MQ=255)
                                                  tCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaatc                                                                                                                         <  1:25955/115‑1 (MQ=255)
                                                  tCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaatc                                                                                                                         >  2:25955/1‑115 (MQ=255)
                                                    gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAtt                                                                                     <  5:301959/149‑1 (MQ=255)
                                                      gcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCTGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATta                                                                                    >  6:335485/1‑148 (MQ=255)
                                                      gcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATta                                                                                    >  2:236341/1‑148 (MQ=255)
                                                            aTTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                           <  4:406551/119‑1 (MQ=255)
                                                                   cACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGc                                                                      <  4:70339/149‑1 (MQ=255)
                                                                     ccGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTa                                                                    <  3:234575/149‑1 (MQ=255)
                                                                        gCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAg                                                                 <  8:36423/149‑1 (MQ=255)
                                                                           cacCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCATTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAAc                                                              <  5:289543/149‑1 (MQ=255)
                                                                            acCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAAc                                                              >  6:262987/1‑148 (MQ=255)
                                                                               gTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                                  <  4:320688/93‑1 (MQ=255)
                                                                               gTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                                  >  3:320688/1‑93 (MQ=255)
                                                                                   tCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTc                                                       <  3:316137/148‑1 (MQ=255)
                                                                                      gCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGa                                                   >  6:207006/1‑149 (MQ=255)
                                                                                       ccAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCt                                                                     >  8:186554/1‑130 (MQ=255)
                                                                                       ccAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCt                                                                     <  7:186554/130‑1 (MQ=255)
                                                                                               tttGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCgg                                          <  7:142905/149‑1 (MQ=255)
                                                                                                   ccACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACt                                      >  6:44946/1‑149 (MQ=255)
                                                                                                    cACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACtt                                     <  1:177601/149‑1 (MQ=255)
                                                                                                        gcggTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCa                                 <  3:168437/149‑1 (MQ=255)
                                                                                                         cggTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCAt                                <  1:236341/149‑1 (MQ=255)
                                                                                                          ggTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATg                               <  5:44946/149‑1 (MQ=255)
                                                                                                           gTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATgg                              <  6:400119/149‑1 (MQ=255)
                                                                                                                      gccgccgc‑cCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg                   >  2:411128/1‑149 (MQ=255)
                                                                                                                       ccgccgc‑cCAGAACCGCAACAGGTTGCTGTTCTGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGa                    <  5:335485/147‑1 (MQ=255)
                                                                                                                       ccgccgc‑cCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGa                    <  7:306007/147‑1 (MQ=255)
                                                                                                                       ccgccgc‑cCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg                   <  5:207006/148‑1 (MQ=255)
                                                                                                                            atgactgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc              >  5:74302/7‑149 (MQ=255)
                                                                                                                               ccAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGc                                                                      >  6:244595/1‑90 (MQ=255)
                                                                                                                               ccAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGc                                                                      <  5:244595/90‑1 (MQ=255)
                                                                                                                                ctgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGAGGTACTTGCCATGGCATTTTCGGAGCCGTc              >  6:269898/3‑145 (MQ=255)
                                                                                                                                cAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCgaga          <  8:225922/149‑1 (MQ=255)
                                                                                                                                   aaCAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATc        <  6:74302/148‑1 (MQ=255)
                                                                                                                                     ccGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAg     <  5:262987/149‑1 (MQ=255)
                                                                                                                                      cGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACTGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCAttt                         >  2:292021/1‑128 (MQ=255)
                                                                                                                                      cGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCgg                                          <  6:174806/111‑1 (MQ=255)
                                                                                                                                        cAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCAc  >  2:373331/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCAC  >  NC_002947/4711385‑4711668

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: ‑