Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,711,518 C→A 6.0% 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 6.0% 202.0 / 3.6 67A8A (GCG→GCT) gpsAglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base A (2/2);  total (29/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 = 8.58e-02

CCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711379‑4711665
                                                                                                                                            |                                                                                                                                                   
ccTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTCCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAgga                                                                                                                                            >  4:70470/1‑148 (MQ=255)
 cTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGtt                                                                                                                                           >  7:117523/1‑149 (MQ=255)
  ttGAGATAGCGCGGATTCTCGCGTTTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTg                                                                                                                                          >  4:234575/1‑149 (MQ=255)
   tGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGc                                                                                                                                         >  3:70339/1‑149 (MQ=255)
      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)
                              cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTg                                                                                                                                          <  3:178487/121‑1 (MQ=255)
                              cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTg                                                                                                                                          >  4:178487/1‑121 (MQ=255)
                               acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGc                                                                                                             <  7:294927/149‑1 (MQ=255)
                                 gCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGcc                                                                                                            >  8:306007/1‑148 (MQ=255)
                                   aTGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAgg                                                                                                                                             >  5:336288/1‑113 (MQ=255)
                                   aTGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAgg                                                                                                                                             <  6:336288/113‑1 (MQ=255)
                                    tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                        <  7:94923/149‑1 (MQ=255)
                                    tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa                                                                                                        <  8:21928/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                                                                                                               >  3:320688/1‑93 (MQ=255)
                                                                                     gTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa                                                                                                               <  4:320688/93‑1 (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)
                                                                                                                                            |                                                                                                                                                   
CCTTGAGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGC‑CCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAG  >  NC_002947/4711379‑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: ‑