Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,711,518 C→A 7.4% A8A (GCG→GCT)  gpsA ← glycerol‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,711,5180CA 7.4% 184.2 / 3.7 62A8A (GCG→GCT) gpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base A (3/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGCGGATACCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711370‑4711665
                                                                                                                                                    |                                                                                                                                                   
ggCGGATACCCTTGAGATAGCCCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAAcc                                                                                                                                                     <  3:148961/149‑1 (MQ=255)
       aCCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAg                                                                                                                                              >  7:96864/1‑149 (MQ=255)
         ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  4:84743/1‑149 (MQ=255)
         ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  6:368596/1‑149 (MQ=255)
         ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGt                                                                                                                                            >  2:250292/1‑149 (MQ=255)
         ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACAGCAAcctgt                                                                                                                                            >  5:393854/1‑145 (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                                                                                   >  3:285745/1‑148 (MQ=255)
                                                                   gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAt                                                                                   <  1:38690/148‑1 (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                                                       >  2:252029/1‑139 (MQ=255)
                                                                                                        agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgt                                                       <  1:252029/139‑1 (MQ=255)
                                                                                                        agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg                                                 >  4:251362/1‑145 (MQ=255)
                                                                                                        agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg                                                 <  3:251362/145‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
GGCGGATACCCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711370‑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: ‑