Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,349 G→C 6.7% M1300I (ATG→ATC)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,3490GC 6.7% 182.5 / 4.0 60M1300I (ATG→ATC) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTT  >  NC_002947/1903208‑1903496
                                                                                                                                             |                                                                                                                                                   
ccTTCGGCTTGACACGGCGTTGTCGCGGGACCATGTCGAGAATATCTATGTTCATCAGCGCATGCTCGAGCAGGGTGGGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGGGATGCGCAGCGCATGGCCAggg                                                                                                                                              <  4:174639/149‑1 (MQ=255)
ccTGCGGCTTGACACGGCGTTCTCGCGGGACAAGGTCGAGATGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGGGATGCGCAGCGCATGGCCAggg                                                                                                                                              <  4:168911/149‑1 (MQ=255)
ccTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCCATGCGCTGCGCATCGCCacac                                                                                                                                              <  7:140909/149‑4 (MQ=255)
      gCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGAtgtg                                                                                                                                         >  3:53143/1‑148 (MQ=255)
       cTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGGCAGGGAtgtgca                                                                                                                                       >  6:385896/1‑147 (MQ=255)
         ggACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTTTGTGTGTGGCCATGCGCTGCGCATCGcc                                                                                                                                                  >  8:9476/1‑136 (MQ=255)
         ggACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATgg                                                                                                                                                    <  4:355428/134‑1 (MQ=255)
         ggACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATgg                                                                                                                                                    >  3:355428/1‑134 (MQ=255)
         ggACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCCATGCGCTGCGCATCGcc                                                                                                                                                  <  7:9476/136‑1 (MQ=255)
             acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGGCTGGCGAGGGAGGTGGATgcggc                                                                                                                                 >  8:175686/1‑149 (MQ=255)
             acGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggc                                                                                                                                 >  7:95070/1‑149 (MQ=255)
              cGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATgcggcg                                                                                                                                <  5:312025/149‑1 (MQ=255)
                    tctcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCg                                                                                                                           <  1:355980/148‑1 (MQ=255)
                    tctcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCg                                                                                                                           >  2:355980/1‑148 (MQ=255)
                      tcGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGt                                                                                                                        >  6:204572/1‑149 (MQ=255)
                      tcGCGGGA‑CAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGt                                                                                                                        >  2:14209/1‑148 (MQ=255)
                        gcgGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCggg                                                                                                                         <  7:298181/146‑1 (MQ=255)
                        gcgGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGa                                                                                                                      >  2:370941/1‑149 (MQ=255)
                        gcgGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCAAGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCggg                                                                                                                         <  5:324137/146‑1 (MQ=255)
                              ccAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCg                                                                                                                <  5:204572/149‑1 (MQ=255)
                                  gTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGa                                                                                                            <  3:265500/149‑1 (MQ=255)
                                  gTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGa                                                                                                            >  5:37048/1‑149 (MQ=255)
                                   tCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAg                                                                                                           >  3:366337/1‑149 (MQ=255)
                                     gagaAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGAGGCGGCGTTGCGGGTGATTGTCGCCGAGCa                                                                                                         >  8:10357/1‑149 (MQ=255)
                                         agaTCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGc                                                                                                     <  8:95070/149‑1 (MQ=255)
                                                tgtgCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCggat                                                                                                <  1:370941/147‑1 (MQ=255)
                                                     agcagcGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATTGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATTTGGATCCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATgcg                                                                                         >  8:86411/1‑149 (MQ=255)
                                                          cgcATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATgcggcggc                                                                                    >  8:342722/1‑149 (MQ=255)
                                                          cgcATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATgcggcggc                                                                                    <  7:342722/149‑1 (MQ=255)
                                                            cATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGcggcggcgg                                                                                  >  3:127802/1‑149 (MQ=255)
                                                               gCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGtt                                                                                                                              >  5:119728/1‑102 (MQ=255)
                                                               gCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGtt                                                                                                                              <  6:119728/102‑1 (MQ=255)
                                                                     gCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGGGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCtt                                                                         <  6:249112/149‑1 (MQ=255)
                                                                     gCAGGGTGCGCTGTTGTGGCAATGGCGGGCGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGAGGTGGATCCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCtt                                                                         <  5:385896/149‑1 (MQ=255)
                                                                        gggTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATg                                                                      <  1:14209/149‑1 (MQ=255)
                                                                         ggTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGt                                                                     >  4:3305/1‑149 (MQ=255)
                                                                              gcTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAg                                                                <  4:284833/149‑1 (MQ=255)
                                                                                tgttgtGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGc                                                              >  4:73297/1‑149 (MQ=255)
                                                                                  ttgtGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCt                                                             <  2:171928/148‑1 (MQ=255)
                                                                                    gtgGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCtt                                                            >  3:319217/1‑147 (MQ=255)
                                                                                     tgGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGa                                                          <  5:378736/148‑1 (MQ=255)
                                                                                      ggCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAg                                                         >  3:369596/1‑148 (MQ=255)
                                                                                                 gAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAgg                                               <  3:73297/147‑1 (MQ=255)
                                                                                                 gAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGc                                              <  2:206262/148‑1 (MQ=255)
                                                                                                   ggCGGGGGCGTATTTCTATGTGGGTGGCGATGCGGAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGtt                                           <  7:175686/149‑1 (MQ=255)
                                                                                                     cGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTTGAAGGCTTGAGCAAGGCCAGGCGTTAt                                         <  6:37048/149‑1 (MQ=255)
                                                                                                         ggCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGgc                                     <  3:3305/149‑1 (MQ=255)
                                                                                                          gCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGgcg                                    <  4:366337/149‑1 (MQ=255)
                                                                                                               tttCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGAtg                               >  7:198505/1‑149 (MQ=255)
                                                                                                                ttCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGAtgt                              >  5:25944/1‑149 (MQ=255)
                                                                                                                       tgtgtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTg                        >  6:313755/1‑148 (MQ=255)
                                                                                                                         tgtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAg                     <  4:319217/149‑1 (MQ=255)
                                                                                                                          gtgGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGt                    >  8:64531/1‑149 (MQ=255)
                                                                                                                           tgcccaTGCGCTGCGCATCGACAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTAACGGCGCGATGTGTACTGAAGTg                   >  6:8973/6‑149 (MQ=255)
                                                                                                                                aTGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCtt                                                                         >  7:7932/1‑90 (MQ=255)
                                                                                                                                aTGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCtt                                                                         <  8:7932/90‑1 (MQ=255)
                                                                                                                                    gcagcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGAt          >  4:154165/1‑149 (MQ=255)
                                                                                                                                      agcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATgg        <  7:86411/149‑1 (MQ=255)
                                                                                                                                      agcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCGGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATgg        >  2:112731/1‑149 (MQ=255)
                                                                                                                                       gcgcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATggg       <  7:64531/149‑1 (MQ=255)
                                                                                                                                         gcaTGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGgt     >  2:48324/1‑149 (MQ=255)
                                                                                                                                            tGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGtgtt  >  2:18567/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGTGTGGCGATGCGCAGCGCATGGCCAGGGATGTGGATGCGGCGTTGCGGGTGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTT  >  NC_002947/1903208‑1903496

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