Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,801,767 G→A 6.8% G51S (GGC→AGC)  rbsA‑I → ribose ABC transporter ‑ ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,801,7670GA 6.8% 174.4 / 2.6 59G51S (GGC→AGC) rbsA‑Iribose ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CTGCCAGGAGAATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTCATGCAG  >  NC_002947/2801620‑2801907
                                                                                                                                                    |                                                                                                                                            
ctgcCAGGAGAATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACgg                                                                                                                                             <  5:65922/149‑1 (MQ=255)
   ccAGGAGAATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGgcg                                                                                                                                           <  2:281535/148‑1 (MQ=255)
      ggAGAATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGCCCGGCGAGAACGGCGCAgg                                                                                                                                       >  5:342034/1‑149 (MQ=255)
         gaATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGTCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCaa                                                                                                                                    <  1:150170/149‑1 (MQ=255)
         gaATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAgg                                                                                                                                       >  1:256750/1‑146 (MQ=255)
         gaATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAgg                                                                                                                                       <  2:256750/146‑1 (MQ=255)
         gaATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCaa                                                                                                                                    <  3:288886/149‑1 (MQ=255)
           aTCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAaga                                                                                                                                  <  7:163103/149‑1 (MQ=255)
           aTCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAaga                                                                                                                                  <  3:174342/149‑1 (MQ=255)
            tCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAagag                                                                                                                                 <  8:117660/149‑1 (MQ=255)
                  aTGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGtt                                                                                                                           <  5:151394/149‑1 (MQ=255)
                  aTGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGtt                                                                                                                           <  4:275985/149‑1 (MQ=255)
                  aTGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGACCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGtt                                                                                                                           >  3:303713/1‑149 (MQ=255)
                       tGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTcc                                                                                                                       >  1:54716/1‑148 (MQ=255)
                               ccAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTAt                                                                                                               >  4:85545/1‑148 (MQ=255)
                                cAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTAt                                                                                                               >  3:88596/1‑147 (MQ=255)
                                 aaTGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCGTGTTGTCGAGGCAGAGCTAGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAg                                                                                                            >  1:305357/1‑149 (MQ=255)
                                   tGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCg                                                                                                          <  6:107395/149‑1 (MQ=255)
                                      aGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGATGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGc                                                                                                        <  6:216664/148‑1 (MQ=255)
                                             cTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGt                                                                                                <  4:88596/149‑1 (MQ=255)
                                              ttGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTg                                                                                               >  3:259164/1‑149 (MQ=255)
                                              ttGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTg                                                                                               <  5:296199/149‑1 (MQ=255)
                                                gccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCGGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGcc                                                                                             >  1:49839/1‑149 (MQ=255)
                                                gccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGc                                                                                                        <  7:352179/138‑1 (MQ=255)
                                                gccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGc                                                                                                        >  8:352179/1‑138 (MQ=255)
                                                 ccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGccc                                                                                            <  1:233016/149‑1 (MQ=255)
                                                 ccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACAGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGcc                                                                                             >  8:158793/1‑148 (MQ=255)
                                                 ccgccAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACAGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGcc                                                                                             >  8:163745/1‑148 (MQ=255)
                                                          gCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCTGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTgg                                                                                                    <  6:29361/132‑1 (MQ=255)
                                                          gCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCTGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTgg                                                                                                    >  5:29361/1‑132 (MQ=255)
                                                            cTGGGCAAGACCTATGCCTAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTca                                                                                 >  1:37787/1‑149 (MQ=255)
                                                                  aaGACCTATGCCTAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGAc                                                                           <  2:37787/149‑1 (MQ=255)
                                                                    gACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCt                                                                         <  3:85545/149‑1 (MQ=255)
                                                                       cTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGTCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTAcc                                                                      >  7:243311/1‑149 (MQ=255)
                                                                       cTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGAc                                                                           >  5:125173/1‑144 (MQ=255)
                                                                       cTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGAc                                                                           <  6:125173/144‑1 (MQ=255)
                                                                       cTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTAcc                                                                      >  1:185475/1‑149 (MQ=255)
                                                                            cccAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGc                                                                 <  6:294367/149‑1 (MQ=255)
                                                                             ccAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGcc                                                                <  4:259164/149‑1 (MQ=255)
                                                                             ccAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGcc                                                                >  8:333410/1‑149 (MQ=255)
                                                                             ccAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGc                                                                 >  2:331337/1‑148 (MQ=255)
                                                                              cAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCa                                                               <  8:243311/149‑1 (MQ=255)
                                                                                 cccGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGc                                                            >  1:84357/1‑149 (MQ=255)
                                                                                     tACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTAc                                                        <  4:303713/149‑1 (MQ=255)
                                                                                                gTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTAt                                                                                                               >  4:371414/1‑83 (MQ=255)
                                                                                                gTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTAt                                                                                                               <  3:371414/83‑1 (MQ=255)
                                                                                                         agcctgCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCgcggc                                      <  1:27449/147‑1 (MQ=255)
                                                                                                         agcctgCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGa                                    <  2:185475/149‑1 (MQ=255)
                                                                                                            ctgCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGc                                 <  2:305357/149‑1 (MQ=255)
                                                                                                            ctgCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACAGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGc                                 <  7:158793/149‑1 (MQ=255)
                                                                                                            ctgCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACAGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGc                                 <  7:163745/149‑1 (MQ=255)
                                                                                                                 gcgcgGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCGGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGaa                             <  2:49839/149‑1 (MQ=255)
                                                                                                                         gAAGTGCTGGCCCTGACCGGCGAGTACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTgg                     >  7:64036/1‑149 (MQ=255)
                                                                                                                                  gCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCAt            >  8:321490/1‑149 (MQ=255)
                                                                                                                                    ccTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATgg          <  6:342034/149‑1 (MQ=255)
                                                                                                                                    ccTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCCCGGCGAGGCCGAACGCCTTGGGGTGCGCATgg          >  3:237896/1‑149 (MQ=255)
                                                                                                                                        aCCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTc        <  7:333410/147‑1 (MQ=255)
                                                                                                                                         ccGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTCATg     <  2:54716/149‑1 (MQ=255)
                                                                                                                                          cGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCAACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTCATGc    >  7:93492/1‑149 (MQ=255)
                                                                                                                                            gCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTCATGCAg  >  1:153623/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CTGCCAGGAGAATCGGCCATGCCTGCATTGGCCAATGAAGTGGTGCTTGCCGCCAGCGGCCTGGGCAAGACCTATGCCCAGCCCGTACTCGGCGAGGTCAGCCTGAGCCTGCGCGCGGGCGAAGTGCTGGCCCTGACCGGCGAGAACGGCGCAGGCAAGAGTACGTTGTCCAAGCTTATCAGCGGCCTGGAGGTGCCCACCACCGGTCACATGACCTACCGCGGCCAGGCCTACGCGCCCGGCAGCCGCGGCGAGGCCGAACGCCTTGGGGTGCGCATGGTCATGCAG  >  NC_002947/2801620‑2801907

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: