Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,453,492 C→G 7.3% G1032G (GGC→GGG nicB → nicotinate dehydrogenase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,453,4920CG 7.3% 214.8 / 4.9 67G1032G (GGC→GGGnicBnicotinate dehydrogenase subunit B
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (3/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGCGAGC  >  NC_002947/4453346‑4453640
                                                                                                                                                   |                                                                                                                                                    
cTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGccc                                                                                                                                                    >  1:126372/1‑149 (MQ=255)
  tGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATTCCCTGGGCCAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCACCTACATCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCgg                                                                                                                                                  >  7:305989/1‑149 (MQ=255)
  tGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCgg                                                                                                                                                  >  6:251487/1‑149 (MQ=255)
   gCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGTACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCATGGGGCCggc                                                                                                                                                 <  3:265766/148‑2 (MQ=255)
    cTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                                <  5:306087/149‑1 (MQ=255)
    cTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCggtg                                                                                                                                                <  6:181445/149‑1 (MQ=255)
    cTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCATGGGGCCggcc                                                                                                                                                <  2:101756/148‑3 (MQ=255)
        cggcATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCa                                                                                                                                            <  2:84496/149‑1 (MQ=255)
          gcATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGt                                                                                                                                          >  1:114222/1‑149 (MQ=255)
            aTGGGCGATGGCGGGGAGGCCCCGGCCTTGTATGCCCTGGGCAAGTCATCGACGCCGGGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGTCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGa                                                                                                                                        <  8:305989/149‑1 (MQ=255)
               ggcGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGGCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCCTCTCCGAACAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCt                                                                                                                                     <  4:194207/147‑1 (MQ=255)
                tCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGAGAGCtt                                                                                                                                    >  6:246969/1‑149 (MQ=255)
                                cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCaaaaa                                                                                                                    <  5:251487/149‑1 (MQ=255)
                                cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCaaa                                                                                                                      >  1:205161/1‑147 (MQ=255)
                                cccGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCaaa                                                                                                                      >  1:201260/1‑147 (MQ=255)
                                    gCCTTTAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGa                                                                                                                <  4:280661/149‑1 (MQ=255)
                                    gCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTGAGTGAGCTTGCCACCCTGCCAAAAAGCGa                                                                                                                <  2:126372/149‑1 (MQ=255)
                                       ttGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGt                                                                                                             >  6:17719/1‑149 (MQ=255)
                                           aTGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTgcgc                                                                                                         <  5:17719/149‑1 (MQ=255)
                                                  gggCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCg                                                                                                  >  7:58454/1‑149 (MQ=255)
                                                   ggCATGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTGCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGc                                                                                                 <  5:331244/149‑1 (MQ=255)
                                                    gCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGcc                                                                                                <  1:282401/149‑1 (MQ=255)
                                                           aTCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTAcc                                                                                         >  6:30802/1‑149 (MQ=255)
                                                            tCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCt                                                                                        >  6:156860/1‑149 (MQ=255)
                                                              gACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCTGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCa                                                                                      <  6:94669/149‑1 (MQ=255)
                                                                    gTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTccc                                                                                <  7:87853/149‑1 (MQ=255)
                                                                       gAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTg                                                                              >  8:116503/1‑148 (MQ=255)
                                                                       gAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTg                                                                              >  5:10119/1‑148 (MQ=255)
                                                                        aGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGa                                                                            <  8:58454/149‑1 (MQ=255)
                                                                         gCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAg                                                                           <  1:328096/149‑1 (MQ=255)
                                                                         gCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAg                                                                           <  2:114222/149‑1 (MQ=255)
                                                                           gAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGa                                                                            <  4:198159/146‑1 (MQ=255)
                                                                           gAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGa                                                                            >  3:198159/1‑146 (MQ=255)
                                                                            aaGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGc                                                                        >  7:25492/1‑149 (MQ=255)
                                                                                cGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGc                                                                    <  5:30802/149‑1 (MQ=255)
                                                                                 gAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCGCAGTTCCCTGGAGGGCGAGcc                                                                   >  4:173005/1‑149 (MQ=255)
                                                                                 gAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCACATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGcc                                                                   >  1:41124/1‑149 (MQ=255)
                                                                                    ttgttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGc                                                                 <  6:59234/148‑1 (MQ=255)
                                                                                      gttCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAgg                                                              <  5:255035/148‑1 (MQ=255)
                                                                                            cTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTgg                                                        <  1:300475/149‑1 (MQ=255)
                                                                                               ccTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgccg                                                     >  8:13928/1‑149 (MQ=255)
                                                                                               ccTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgccg                                                     >  8:4102/1‑149 (MQ=255)
                                                                                               ccTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgcc                                                      <  2:205161/148‑1 (MQ=255)
                                                                                                cacAGCACCGGCGGGTAAGAAAAACACGGCGTGGCGGCCCGCCCCCAGGGGGCCGGGGGCAGGGGGGGTGGCCCCCGGCCAAAAACCGACCTGGGCGCCCACGCCCCCAACCTCACGTACCATGAGGGCGCGCCCCCGCCACAGgccgc                                                    >  6:109512/3‑149 (MQ=255)
                                                                                                cTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGgccgc                                                    >  2:97905/1‑149 (MQ=255)
                                                                                                     cACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATg                                               >  6:295205/1‑149 (MQ=255)
                                                                                                     cACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATg                                               >  3:118068/1‑149 (MQ=255)
                                                                                                       ccGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGcc                                             >  6:21746/1‑149 (MQ=255)
                                                                                                            ttCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGcc                                        >  3:178265/1‑149 (MQ=255)
                                                                                                                 cGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAgg                                   >  6:150938/1‑149 (MQ=255)
                                                                                                                  gAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGt                                  <  8:25492/149‑1 (MQ=255)
                                                                                                                  gAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGt                                  >  7:159766/1‑149 (MQ=255)
                                                                                                                   aaaaaCACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTg                                 <  7:13928/149‑1 (MQ=255)
                                                                                                                   aaaaaCACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTg                                 <  5:295205/149‑1 (MQ=255)
                                                                                                                      aaCACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGac                              <  2:201260/149‑1 (MQ=255)
                                                                                                                       acacGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGaca                             >  4:137592/1‑149 (MQ=255)
                                                                                                                       acacGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGaca                             >  6:172896/1‑149 (MQ=255)
                                                                                                                          cGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGAcacgc                          <  5:109512/149‑1 (MQ=255)
                                                                                                                                tgaccaccgggcCCATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACgtg                     >  3:154674/7‑149 (MQ=255)
                                                                                                                                  cggccggcCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCg                  <  7:116503/149‑1 (MQ=255)
                                                                                                                                    gccggccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCt                >  2:227818/1‑149 (MQ=255)
                                                                                                                                    gccggccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCt                <  5:156860/149‑1 (MQ=255)
                                                                                                                                    caccgggcCCATGGGGCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGc                 >  4:120958/3‑149 (MQ=255)
                                                                                                                                       cggccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGa              <  3:137592/149‑1 (MQ=255)
                                                                                                                                       cggccCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGa              <  4:131926/149‑1 (MQ=255)
                                                                                                                                              aTGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACgg       <  1:83151/149‑1 (MQ=255)
                                                                                                                                               tGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGAGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGc      <  8:130894/149‑1 (MQ=255)
                                                                                                                                                gggCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGCg     >  3:57783/1‑149 (MQ=255)
                                                                                                                                                   cccGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGCGAGc  >  4:153455/1‑149 (MQ=255)
                                                                                                                                                   cccGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGCGAGc  <  4:118068/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTGCTGGCGGCATGGTCGATGGCTGGGAGGCCCCGGCCTTGAATGCCCTGGGCAAGTCATCGACGCCGTGGAGCGAAGACGAGTTGTTCAACTACCTCAGCACCGGCTTCTCCGAAAAACACGGCGTGGCGGCCGGCCCCATGGGCCCGGTGGTCAGTGAGCTTGCCACCCTGCCAAAAAGCGACGTGCGCGCCATCGCCCACTACCTCAGTTCCCTGGAGGGCGAGCCGCAGGCACTGGCCGCCAATGCCGCGCCCCAGGTGGACACGCACGTGTCGCTGAGCAACGGCGAGC  >  NC_002947/4453346‑4453640

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: