Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,829,819 C→T 8.2% N141N (AAC→AAT)  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,8190CT 8.2% 171.1 / 2.6 60N141N (AAC→AAT) PP_2482molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGG  >  NC_002947/2829671‑2829967
                                                                                                                                                    |                                                                                                                                                    
gCCTGACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAAc                                                                                                                                                      >  4:263125/1‑149 (MQ=255)
 ccTGACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACg                                                                                                                                                     <  3:263125/149‑1 (MQ=255)
     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                                 <  5:214339/133‑1 (MQ=255)
                                                                                                                                     gCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCGTGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTAc                                 >  6:214339/1‑133 (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)
                                                                                                                                                    |                                                                                                                                                    
GCCTGACCACCAATGGCCAGTTGCTGGCCGGCAAACTGCCGTTATTGCTCAGCGCCGGCATCCAGCGGCTCAATGTGTCCCTCGATACGCTCGACCCGAGTGCCTTCCGCACGATCGCCAAGGGCGGCGACCTGCAGACGGTGCTCAACGGCATTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGG  >  NC_002947/2829671‑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: ‑