Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,160,444 C→G 5.6% L86V (CTG→GTG) ‡ gtsC → mannose/glucose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,160,4440CG 5.6% 233.4 / 4.1 71L86V (CTG→GTG) ‡gtsCmannose/glucose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base G (2/2);  total (34/37)
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.73e-01

CTGACCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCG  >  NC_002947/1160297‑1160588
                                                                                                                                                   |                                                                                                                                                
ctgACCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGc                                                                                                                                                  <  4:82263/148‑1 (MQ=255)
ctgACCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCt                                                                                                                                                 <  2:382740/149‑1 (MQ=255)
   aCCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGTAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTCCCCGCGGTGCTGa                                                                                                                                               >  4:379133/1‑148 (MQ=255)
   aCCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGAt                                                                                                                                              <  8:64975/149‑1 (MQ=255)
    ccAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGCCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGAtc                                                                                                                                             >  8:62708/1‑149 (MQ=255)
                cccccGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGGAGATCAGCACcgcgggc                                                                                                                                 <  6:280163/149‑4 (MQ=255)
                 ccccGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTg                                                                                                                                <  8:331328/149‑1 (MQ=255)
                   ccGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGAt                                                                                                                                              >  6:157179/1‑133 (MQ=255)
                   ccGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGAt                                                                                                                                              <  5:157179/133‑1 (MQ=255)
                   ccGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGt                                                                                                                              >  3:80541/1‑149 (MQ=255)
                    cGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGATCTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTg                                                                                                                             >  4:274748/1‑149 (MQ=255)
                    cGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTg                                                                                                                             <  3:5946/149‑1 (MQ=255)
                     agcagcATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGc                                                                                                                            <  4:310606/144‑1 (MQ=255)
                       aGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCAc                                                                                                                          <  2:358893/149‑1 (MQ=255)
                        gACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGGAGATCAGCACcgcgggc                                                                                                                                 >  8:79051/1‑138 (MQ=255)
                        gACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGGAGATCAGCACcgcgggc                                                                                                                                 <  7:79051/141‑4 (MQ=255)
                           aTCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGaa                                                                                                                      >  8:108947/1‑149 (MQ=255)
                           aTCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGaa                                                                                                                      >  5:79641/1‑149 (MQ=255)
                              accaccGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACgg                                                                                                                   >  6:325976/1‑149 (MQ=255)
                              accaccGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACg                                                                                                                    <  4:80541/148‑1 (MQ=255)
                                 accGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTa                                                                                                                >  8:35816/1‑149 (MQ=255)
                                       aaCCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCt                                                                                                          >  4:360164/1‑149 (MQ=255)
                                        aCCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGGCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATATGCTg                                                                                                         <  8:357425/149‑1 (MQ=255)
                                         ccTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGt                                                                                                        >  2:125842/1‑149 (MQ=255)
                                          ctgctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGt                                                                                                        <  6:298946/148‑1 (MQ=255)
                                          ctgctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGt                                                                                                        >  5:298946/1‑148 (MQ=255)
                                           tgctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCt                                                                                                                 >  1:144399/1‑138 (MQ=255)
                                           tgctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCt                                                                                                                 <  2:144399/138‑1 (MQ=255)
                                           tgctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTc                                                                                                       >  2:213116/1‑148 (MQ=255)
                                            gctgAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGGTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGa                                                                                                     <  5:219510/149‑1 (MQ=255)
                                                aGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGAtgtg                                                                                                 <  2:15831/149‑1 (MQ=255)
                                                aGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGAtg                                                                                                   <  2:217403/147‑1 (MQ=255)
                                                   tGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGgcg                                                                                              >  5:37235/1‑149 (MQ=255)
                                                       cTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTc                                                                                          >  8:296831/1‑149 (MQ=255)
                                                          cAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGa                                                                                       <  1:213116/149‑1 (MQ=255)
                                                            gTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAgg                                                                                     <  8:259030/149‑1 (MQ=255)
                                                              gATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGtt                                                                                   >  8:154651/1‑149 (MQ=255)
                                                              gATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGtt                                                                                   >  5:287727/1‑149 (MQ=255)
                                                                 tACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCg                                                                                 >  5:158467/1‑148 (MQ=255)
                                                                   ccGGCATTGGCGGGGTAATGGCGTGGGGAGCGGTGATCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCaa                                                                              <  3:379133/149‑1 (MQ=255)
                                                                        aTAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGtt                                                                         <  1:364009/149‑1 (MQ=255)
                                                                             cTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCg                                                                    <  3:360164/149‑1 (MQ=255)
                                                                                  tCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCctg                                                               >  2:309118/1‑149 (MQ=255)
                                                                                     aGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCctgctg                                                            <  6:79641/149‑1 (MQ=255)
                                                                                     aGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCctgctg                                                            <  8:224200/149‑1 (MQ=255)
                                                                                      ggCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCctgctgc                                                           >  2:39087/1‑149 (MQ=255)
                                                                                           ggggAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCctgcttgaatgc                                                      >  2:141044/1‑142 (MQ=255)
                                                                                            gggAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGtt                                                       <  7:296831/147‑1 (MQ=255)
                                                                                            gggAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTACTTCGGCCTGCTGCTGtt                                                       <  7:154651/147‑1 (MQ=255)
                                                                                             cgAGCCGTCACCGGTTACTTCTCCAACTCGATCATGATCACCGTGCCCTCGCTGCTCATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTCTCGATGTGCCGCTTCCGAGGTTCGCAACTGATCTTCGGCCTGCTGCTGTTCgg                                                    <  4:83283/148‑1 (MQ=255)
                                                                                              gAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGc                                                   >  1:39468/1‑149 (MQ=255)
                                                                                              gAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGc                                                   <  7:108947/149‑1 (MQ=255)
                                                                                                  ggTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGtt                                               <  2:39468/149‑1 (MQ=255)
                                                                                                   gTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGttt                                              >  4:250196/1‑149 (MQ=255)
                                                                                                    tCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTc                                             <  1:309118/149‑1 (MQ=255)
                                                                                                    tCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCACGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGATGATCGGCTGttta                                             >  6:105060/1‑148 (MQ=255)
                                                                                                      aGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGAt                                                                                                    <  5:272614/92‑1 (MQ=255)
                                                                                                      aGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGAt                                                                                                    >  6:272614/1‑92 (MQ=255)
                                                                                                      aGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCt                                           <  6:287727/149‑1 (MQ=255)
                                                                                                       gCGGTTACTTCTGGAACTAGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTg                                          >  2:258315/1‑149 (MQ=255)
                                                                                                              cTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTc                                    <  7:35816/148‑1 (MQ=255)
                                                                                                                          gATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGtgctg                        <  5:325976/148‑1 (MQ=255)
                                                                                                                              aTGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGcc                   <  6:37235/149‑1 (MQ=255)
                                                                                                                              aTGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGcc                   >  4:188326/1‑149 (MQ=255)
                                                                                                                                      cGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGt           >  6:365941/1‑149 (MQ=255)
                                                                                                                                      cGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCAGATGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGt           >  1:316521/1‑149 (MQ=255)
                                                                                                                                        caatcgcgGTGGAGATCAGCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTTCCGGCCTCGtt          >  4:253943/5‑148 (MQ=255)
                                                                                                                                         gCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCa        >  3:266049/1‑149 (MQ=255)
                                                                                                                                         gCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGGTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCa        <  3:274748/149‑1 (MQ=255)
                                                                                                                                          cccGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCa        <  3:188326/148‑1 (MQ=255)
                                                                                                                                               ggTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCg  >  4:23633/1‑149 (MQ=255)
                                                                                                                                               ggTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCGGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCg  <  6:158467/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGACCAGCTTCAAGACCCCCGAAGACATCACCACCGGCAACCTGCTGAGCTGGCCTGCAGTGATTACCGGCATAGGCTGGGTCAAGGCCTGGGGAGCGGTCAGCGGTTACTTCTGGAACTCGATCATGATCACCGTGCCCGCGGTGCTGATCTCCACCGCGATTGGTGCACTGAACGGCTATGTGCTGTCGATGTGGCGCTTCCGAGGTTCGCAACTGTTCTTCGGCCTGCTGCTGTTCGGCTGTTTCCTGCCGTTCCAGACCGTGCTGCTGCCGGCCTCGTTCACCCTCG  >  NC_002947/1160297‑1160588

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