Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,829,821 G→C 8.2% G142A (GGC→GCC)  PP_2482 → molybdenum cofactor biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,829,8210GC 8.2% 180.1 / 4.9 60G142A (GGC→GCC) PP_2482molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base C (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGG  >  NC_002947/2829676‑2829967
                                                                                                                                                 |                                                                                                                                                  
accaccAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCAt                                                                                                                                                 <  6:163931/149‑1 (MQ=255)
 ccaccaATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGAGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCAtt                                                                                                                                                <  7:270275/149‑1 (MQ=255)
    ccaATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAg                                                                                                                                             >  3:189374/1‑149 (MQ=255)
      aaTGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAg                                                                                                                                             >  7:131745/1‑147 (MQ=255)
       aTGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAATGCCGTTGAGCac                                                                                                                                          <  5:106086/149‑2 (MQ=255)
        tGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATCCGCTCGACCCGCGGGCCTTCCTCACGATCGCCAAGGTGGGCGACCTGCAGACTGTTCTCAATGCCGTTTAGCacc                                                                                                                                         <  2:93147/149‑3 (MQ=255)
          gCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGa                                                                                                                                              >  8:49660/1‑142 (MQ=255)
          gCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGa                                                                                                                                              <  7:49660/142‑1 (MQ=255)
          gCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGgcg                                                                                                                                       >  3:24114/1‑149 (MQ=255)
                        ggCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGAGGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCAt                                                                                                                         >  2:197810/1‑149 (MQ=255)
                           aaaCTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAg                                                                                                                      <  5:6946/149‑1 (MQ=255)
                           aaaCTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAg                                                                                                                      <  4:189374/149‑1 (MQ=255)
                              cTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCAt                                                                                                                   >  2:163056/1‑149 (MQ=255)
                                 ccGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGCTACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATc                                                                                                                  >  6:106427/1‑147 (MQ=255)
                                 ccGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATTCGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCaa                                                                                                                <  5:106427/149‑1 (MQ=255)
                                         gCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACa                                                                                                        >  3:40259/1‑149 (MQ=255)
                                                 ccGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCg                                                                                                <  4:24114/149‑1 (MQ=255)
                                                       tCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGt                                                                                          <  5:165813/149‑1 (MQ=255)
                                                        ccAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTg                                                                                         >  8:1788/1‑149 (MQ=255)
                                                          aGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTgg                                                                                        >  1:180396/1‑148 (MQ=255)
                                                           gCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTgg                                                                                        >  6:27015/1‑147 (MQ=255)
                                                                tCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAAc                                                                                 >  3:277141/1‑149 (MQ=255)
                                                                        tCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATc                                                                          >  6:47857/1‑148 (MQ=255)
                                                                         cccTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATc                                                                          >  2:324100/1‑147 (MQ=255)
                                                                          ccTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAg                                                                       <  7:1788/149‑1 (MQ=255)
                                                                           cTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGt                                                                      >  7:283924/1‑149 (MQ=255)
                                                                            tCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGtt                                                                     >  4:172931/1‑149 (MQ=255)
                                                                              gATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGtttt                                                                   <  2:180396/149‑1 (MQ=255)
                                                                                taggCTCGACCCGAGTGCCTTCCGCACGATCGGCTAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGc                                                                 <  4:277141/146‑1 (MQ=255)
                                                                                               tGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCt                                                  <  4:40259/149‑1 (MQ=255)
                                                                                                 ccTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGc                                                <  1:163056/149‑1 (MQ=255)
                                                                                                   ttCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGtt                                                                     >  6:249784/1‑126 (MQ=255)
                                                                                                   ttCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGtt                                                                     <  5:249784/126‑1 (MQ=255)
                                                                                                     ccGCACGATCGCCAAGGGCGGCGACCTGCAGAGGGTGCTCAACGTCATTGAGGAGGCGCGCCTAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCATGTTTTGCCCATGTTGGCGTTCTGCCtt                                             <  1:197810/148‑1 (MQ=255)
                                                                                                      cGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCtt                                             >  8:240796/1‑147 (MQ=255)
                                                                                                       gCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcag                                          >  5:16753/1‑149 (MQ=255)
                                                                                                        cACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcag                                          >  1:333543/1‑148 (MQ=255)
                                                                                                         aCGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcag                                          >  1:108045/1‑147 (MQ=255)
                                                                                                             tCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGTCCATGTTGGCGTTCTGCCTTCAGCAGGGc                                    <  8:131745/149‑1 (MQ=255)
                                                                                                                ccAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 <  8:283924/149‑1 (MQ=255)
                                                                                                                ccAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 >  1:318580/1‑149 (MQ=255)
                                                                                                                  aaGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 >  5:104988/1‑147 (MQ=255)
                                                                                                                       cggcgACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGc                          <  7:240796/149‑1 (MQ=255)
                                                                                                                          cgACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGtt                       <  2:318580/149‑1 (MQ=255)
                                                                                                                           gACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGttt                      >  3:105959/1‑149 (MQ=255)
                                                                                                                               tGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATc                  >  6:152697/1‑149 (MQ=255)
                                                                                                                                gCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCGTGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 >  6:214339/1‑133 (MQ=255)
                                                                                                                                gCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCGTGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 <  5:214339/133‑1 (MQ=255)
                                                                                                                                gCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcagcag                                       >  2:351376/1‑127 (MQ=255)
                                                                                                                                  aGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTTTCGaa               <  2:108045/149‑1 (MQ=255)
                                                                                                                                   gACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAAc              >  5:93852/1‑149 (MQ=255)
                                                                                                                                   gACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAAc              <  2:333543/149‑1 (MQ=255)
                                                                                                                                       ctGCTCAATGCCGTTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTAt          >  3:104665/2‑149 (MQ=255)
                                                                                                                                         gCTCAATGCCGTTGAGCAGGCGTGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcagcag                                       <  5:352834/118‑1 (MQ=255)
                                                                                                                                         gCTCAATGCCGTTGAGCAGGCGAGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTcagcag                                       >  6:352834/1‑118 (MQ=39)
                                                                                                                                          cTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATgc        <  1:324100/148‑1 (MQ=255)
                                                                                                                                           tCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATgc        >  6:221575/1‑147 (MQ=255)
                                                                                                                                           tCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATgc        <  5:221575/147‑1 (MQ=255)
                                                                                                                                             aaCGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCAt    >  8:242289/1‑149 (MQ=255)
                                                                                                                                               cGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATgg  <  6:16753/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGG  >  NC_002947/2829676‑2829967

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