Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,944 G→C 6.1% G78R (GGC→CGC)  moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,9440GC 6.1% 246.1 / 4.0 77G78R (GGC→CGC) moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base G (38/34);  new base C (2/3);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

AGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTC  >  NC_002947/2421796‑2422090
                                                                                                                                                    |                                                                                                                                                   
aGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAgg                                                                                                                                                     <  8:66137/149‑1 (MQ=255)
     cACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAg                                                                                                                                                <  4:164419/149‑1 (MQ=255)
             cgcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGcc                                                                                                                                       >  1:122009/1‑149 (MQ=255)
              gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                      <  1:386366/149‑1 (MQ=255)
               cgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGTCCGAACAGCGCCATTGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCgg                                                                                                                                     <  8:218783/149‑1 (MQ=255)
               cgAGAACGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGCACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTCACTGGCAGCGCCAGCCGCTGCCgg                                                                                                                                     >  6:6493/1‑149 (MQ=255)
                  gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                                  <  5:415748/149‑1 (MQ=255)
                   aCCGAAAACCTGCCGCTGTCCGAAGCCGATGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCTCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTcag                                                                                                                                 <  3:300485/149‑1 (MQ=38)
                   aCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTcag                                                                                                                                 <  3:305748/149‑1 (MQ=39)
                     cGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTcagc                                                                                                                                >  7:228459/1‑148 (MQ=255)
                      gAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTcagc                                                                                                                                >  5:87147/1‑147 (MQ=255)
                      gAAAACCTGCCGCTGGCCGAAGCCGACGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTcagcca                                                                                                                              >  3:104068/1‑149 (MQ=255)
                           ccTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGACCGGcagcggctggc                                                                                                                         <  8:72618/149‑8 (MQ=255)
                             tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtt                                                                                                                       >  4:164733/1‑149 (MQ=255)
                                     gCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGcc                                                                                                                 <  6:255492/147‑1 (MQ=255)
                                       cGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGcc                                                                                                             <  8:228459/149‑1 (MQ=255)
                                         aaGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGAGCTGCGGCCGTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGGCGGTCAGCCAGCGTATTTTCGCCGGCCAt                                                                                                           <  5:6493/149‑1 (MQ=255)
                                         aaGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCGCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCAt                                                                                                           <  6:413420/149‑1 (MQ=255)
                                           gCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATgc                                                                                                         <  2:122009/149‑1 (MQ=255)
                                               aGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGccc                                                                                                     >  2:101502/1‑149 (MQ=255)
                                                gggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGcccc                                                                                                    <  4:45483/149‑1 (MQ=255)
                                                 ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCa                                                                                                   <  8:324494/149‑1 (MQ=255)
                                                 ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCa                                                                                                   <  4:5246/149‑1 (MQ=255)
                                                  gCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAg                                                                                                  >  6:180015/1‑149 (MQ=255)
                                                   ccGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGc                                                                                                 >  7:316862/1‑149 (MQ=255)
                                                    cgcgTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGcc                                                                                                >  2:281724/1‑149 (MQ=255)
                                                             gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCaa                                                                                         <  6:87147/147‑1 (MQ=255)
                                                                aCCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCggg                                                                                    >  2:85812/1‑149 (MQ=255)
                                                                  cGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTa                                                                                  >  7:5820/1‑149 (MQ=255)
                                                                   gACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTAc                                                                                 <  1:8100/149‑1 (MQ=255)
                                                                   gACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTAc                                                                                 <  1:185230/149‑1 (MQ=255)
                                                                      cTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTg                                                                              >  8:319923/1‑149 (MQ=255)
                                                                                 gctcgACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCTGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGAt                                                                     <  7:319923/147‑1 (MQ=255)
                                                                                     gACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCAc                                                               >  8:149791/1‑149 (MQ=255)
                                                                                       ccTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCg                                                             <  2:322215/149‑1 (MQ=255)
                                                                                          gccgccCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGgcg                                                          >  1:397295/1‑149 (MQ=255)
                                                                                           ccgccCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGgcgc                                                         >  1:372737/1‑149 (MQ=255)
                                                                                              cccTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGcc                                                      >  5:88654/1‑149 (MQ=255)
                                                                                               ccTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCt                                                     >  2:132374/1‑149 (MQ=255)
                                                                                                cTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCtg                                                    >  2:290067/1‑149 (MQ=255)
                                                                                                   gCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGc                                                                                              >  1:273521/1‑104 (MQ=255)
                                                                                                   gCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGc                                                                                              <  2:273521/104‑1 (MQ=255)
                                                                                                   gCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGc                                                 >  7:148069/1‑149 (MQ=255)
                                                                                                   gCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGc                                                 <  4:104068/149‑1 (MQ=255)
                                                                                                    ccGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGCCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGTATCTTCACCGGCGCGCCTGTGcc                                                <  5:180015/149‑1 (MQ=255)
                                                                                                     cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                               >  1:362847/1‑149 (MQ=255)
                                                                                                       aaCAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                               >  3:380972/1‑147 (MQ=255)
                                                                                                            cgcCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACAGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAggggg                                        >  4:283098/1‑149 (MQ=255)
                                                                                                                aTGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGa                                    <  1:101502/149‑1 (MQ=255)
                                                                                                                aTGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGa                                    <  4:310953/149‑1 (MQ=255)
                                                                                                                 tGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGAc                                   <  2:362847/149‑1 (MQ=255)
                                                                                                                   gACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACt                                  <  3:164733/148‑1 (MQ=255)
                                                                                                                    aCGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCGTGTGCCGGAGGGGGCCGACTGc                                >  8:67285/1‑149 (MQ=255)
                                                                                                                        ttATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCg                            <  7:149791/149‑1 (MQ=255)
                                                                                                                             ccTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACt                                  >  5:86915/1‑138 (MQ=255)
                                                                                                                             ccTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACt                                  <  6:86915/138‑1 (MQ=255)
                                                                                                                              cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGACTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGc                                       >  1:42246/1‑133 (MQ=255)
                                                                                                                              cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGc                                       <  2:42246/132‑1 (MQ=255)
                                                                                                                              cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                       >  8:386789/1‑148 (MQ=255)
                                                                                                                              cTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGc                      <  8:316862/149‑1 (MQ=255)
                                                                                                                              cTTGCGCCTGGCAGACTTGCAGGGCAAGACCCTGCCGGTCAGCCATCGTATTTTCGACGGCCATGCGCCCCAGACGCTGCAACCGGGTACCTGCGAGCGGATCTTCACAGGCGCGCCTGTGCCTGACGCCGCCGACTGCGTCcagaag                       >  8:106943/1‑146 (MQ=255)
                                                                                                                               ttGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATg                       >  7:11468/1‑147 (MQ=255)
                                                                                                                                 gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCa                     <  8:5820/147‑1 (MQ=255)
                                                                                                                                 gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAg                    >  3:333887/1‑148 (MQ=255)
                                                                                                                                  cgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa                  <  7:386789/149‑1 (MQ=255)
                                                                                                                                   gcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa                  >  4:329488/1‑148 (MQ=255)
                                                                                                                                        gCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAacac            >  1:90236/1‑149 (MQ=255)
                                                                                                                                        gCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAaca             >  7:39373/1‑148 (MQ=255)
                                                                                                                                        gCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAaca             <  8:39373/148‑1 (MQ=255)
                                                                                                                                        gCTGACCGGCAGCGGCTGGCCCTGCCAGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAacac            >  3:359533/1‑149 (MQ=39)
                                                                                                                                          tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc           >  6:102573/1‑148 (MQ=255)
                                                                                                                                            aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGaa        <  1:85812/149‑1 (MQ=255)
                                                                                                                                            aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGaa        <  1:290067/149‑1 (MQ=255)
                                                                                                                                              tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg       >  3:208040/1‑148 (MQ=255)
                                                                                                                                                gCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTAc    >  4:232423/1‑149 (MQ=255)
                                                                                                                                                 cAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGc                      <  5:22835/130‑1 (MQ=255)
                                                                                                                                                 cAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGc                      >  6:22835/1‑130 (MQ=255)
                                                                                                                                                  aGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGTCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTc  >  7:369209/1‑149 (MQ=255)
                                                                                                                                                    ggCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaa                <  2:143719/133‑1 (MQ=255)
                                                                                                                                                    ggCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGaaa                >  1:143719/1‑133 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAGTACTC  >  NC_002947/2421796‑2422090

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: