Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,103,088 A→T 6.7% intergenic (+20/‑170) phhB → / → PP_5694 pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,103,0880AT 6.7% 257.1 / 5.1 84intergenic (+20/‑170)phhB/PP_5694pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base T (4/2);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CTGACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACT  >  NC_002947/5102943‑5103235
                                                                                                                                                 |                                                                                                                                                   
ctgactgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCagag                                                                                                                                                  <  6:150356/149‑1 (MQ=255)
     tgaATGGGGCAAGGTCACCGTGACCTGGTGGAGACACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGcccc                                                                                                                                             >  2:241692/1‑149 (MQ=255)
       aaTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCtt                                                                                                                                           <  3:297031/149‑1 (MQ=255)
           gggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGGAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                       >  8:83452/1‑149 (MQ=255)
               aaGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                                   <  6:151562/149‑1 (MQ=255)
                aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCCGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgg                                                                                                                                  <  6:13703/149‑1 (MQ=255)
                aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgg                                                                                                                                  <  2:66609/149‑1 (MQ=255)
                 ggTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                 <  7:196229/149‑1 (MQ=255)
                 ggTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                 <  2:115954/149‑1 (MQ=255)
                  gTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                                >  4:350886/1‑149 (MQ=255)
                   tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCCGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                                >  4:284002/1‑148 (MQ=255)
                       cGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCatat                                                                                                                           <  1:241692/149‑1 (MQ=255)
                       cGTGACCTGGTGGAGCCACTCGATCAAAGGCATGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGACGCCCCTTTTTTCATGGGGCatat                                                                                                                           >  6:238860/1‑149 (MQ=255)
                                gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                                    <  5:289728/147‑1 (MQ=255)
                                   gAGCCACTCGATCAAAGGCCCGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTc                                                                                                                <  3:284002/148‑1 (MQ=255)
                                     gCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGcc                                                                                                             >  7:101651/1‑149 (MQ=255)
                                        aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCAAGCACTGACAAGGTGTCGGAAACGGCTGAAGGCCTGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGAGCATATGTGCGCTGTTCGacct                                                                                                           >  8:134854/1‑147 (MQ=255)
                                         cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGc                                                                                                          >  5:48862/1‑148 (MQ=255)
                                          tCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCGAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCa                                                                                                         >  2:353494/1‑148 (MQ=255)
                                            gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaaaa                                                                                                      <  3:350886/149‑1 (MQ=255)
                                             aTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAc                                                                                                     >  3:107642/1‑149 (MQ=255)
                                              tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCgaggcc                                                                                                      >  1:277261/1‑142 (MQ=255)
                                              tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAcc                                                                                                    <  6:48862/149‑1 (MQ=255)
                                              tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGAGCATATGTGCGCTGTACGCCCGCAAAAcc                                                                                                    >  5:335463/1‑149 (MQ=255)
                                                      cTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAAc                                                                                             >  7:181181/1‑148 (MQ=255)
                                                       tGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGa                                                                                           >  5:42725/1‑149 (MQ=255)
                                                        gCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGAc                                                                                          <  7:43454/149‑1 (MQ=255)
                                                        gCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCATTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCCACAAACGAc                                                                                          >  3:192706/1‑149 (MQ=255)
                                                         cACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGAcc                                                                                         <  1:350870/149‑1 (MQ=255)
                                                             gCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACCACCggg                                                                                      >  7:232448/1‑148 (MQ=255)
                                                              cAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGa                                                                                     >  6:265659/1‑148 (MQ=255)
                                                                  gACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGaaaaa                                                                                 >  4:99915/1‑148 (MQ=255)
                                                                   aCTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGaaaaaaa                                                                               <  4:21996/149‑1 (MQ=255)
                                                                         tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                         >  6:220118/1‑149 (MQ=255)
                                                                         tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                         >  1:349464/1‑149 (MQ=255)
                                                                           atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAACACCGACAAACGCCCGGGAAAAAAACCAAACTg                                                                       >  8:273240/1‑149 (MQ=255)
                                                                           atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTg                                                                       <  2:359023/149‑1 (MQ=255)
                                                                           atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTg                                                                       >  5:283307/1‑149 (MQ=255)
                                                                              tgtgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCa                                                                    >  1:24434/1‑149 (MQ=255)
                                                                                   acgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                  >  8:278790/1‑146 (MQ=255)
                                                                                    cgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGa                                                              <  8:101651/149‑1 (MQ=255)
                                                                                     gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGAc                                                             >  5:317847/1‑149 (MQ=255)
                                                                                         tGACAAGGTGGCGGAAAAGTCTGAAGGCCGTAAGTAAACACCAGAAGTAAGGTGCCAGAGGACACTTTTTTCATGGGGCATATGTGCGCTGTTAGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                         <  5:238860/149‑1 (MQ=255)
                                                                                          gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgtg                                                        >  5:144501/1‑149 (MQ=255)
                                                                                           aCAAGGTGGCGGAAACGGCTGGAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                         <  7:83452/147‑1 (MQ=255)
                                                                                           aCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                  <  6:231719/138‑1 (MQ=255)
                                                                                           aCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                  >  5:231719/1‑138 (MQ=255)
                                                                                           aCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                         <  5:265659/147‑1 (MQ=255)
                                                                                              aGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                    <  6:317847/149‑1 (MQ=255)
                                                                                              aGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                    <  6:144501/149‑1 (MQ=255)
                                                                                                 tGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCggg                                                                                      <  7:76770/112‑1 (MQ=255)
                                                                                                 tGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCggg                                                                                      >  8:76770/1‑112 (MQ=255)
                                                                                                   gCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaaa                                                                          <  2:4075/122‑1 (MQ=255)
                                                                                                    cGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGccc                                              >  6:51715/1‑149 (MQ=255)
                                                                                                       aaaCGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAg                                           <  8:181181/149‑1 (MQ=255)
                                                                                                             cTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCttttt                                     <  6:42725/149‑1 (MQ=255)
                                                                                                             cTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCttttt                                     <  6:283307/149‑1 (MQ=255)
                                                                                                              tGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCttttt                                     >  7:326341/1‑148 (MQ=255)
                                                                                                                aaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGa                                   <  3:99915/148‑1 (MQ=255)
                                                                                                                aaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGa                                   <  4:107642/148‑1 (MQ=255)
                                                                                                                      ggAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCCAAAAACCACCGGGAAAAAAACCAAACTCTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTAc                             >  7:281221/1‑148 (MQ=255)
                                                                                                                         aGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAc                         >  3:87808/1‑149 (MQ=255)
                                                                                                                          gTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGccc                                                                                                            <  6:117832/65‑1 (MQ=255)
                                                                                                                          gTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGccc                                                                                                            >  5:117832/1‑65 (MQ=255)
                                                                                                                          gTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                        >  7:46371/1‑149 (MQ=255)
                                                                                                                             aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACgcgc                     >  6:290412/1‑149 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                               >  6:55416/1‑134 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                               <  5:55416/134‑1 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                >  7:1188/1‑149 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                <  8:46371/149‑1 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                >  6:361667/1‑149 (MQ=255)
                                                                                                                                  cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTACGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAATGTGATACAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                >  1:305554/1‑149 (MQ=255)
                                                                                                                                  atgaaaaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                >  6:350105/8‑149 (MQ=255)
                                                                                                                                   tgaaaaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAATGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  2:86231/7‑149 (MQ=255)
                                                                                                                                     aaggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  4:165257/1‑147 (MQ=255)
                                                                                                                                     aaaaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  4:359641/5‑147 (MQ=255)
                                                                                                                                       aaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGa           >  6:209916/3‑149 (MQ=255)
                                                                                                                                        gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCatat                                                                                                                           >  6:229630/1‑36 (MQ=255)
                                                                                                                                        gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCatat                                                                                                                           <  5:229630/36‑1 (MQ=255)
                                                                                                                                         aaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCttt                                       >  6:329561/1‑119 (MQ=255)
                                                                                                                                         aaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCttt                                       <  5:329561/119‑1 (MQ=255)
                                                                                                                                          aggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt        <  5:51715/149‑1 (MQ=255)
                                                                                                                                              gCCTGAGGCCGCTTTTTTCATGGGGAATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtca    <  4:192706/149‑1 (MQ=255)
                                                                                                                                              gCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtca    >  3:352012/1‑149 (MQ=255)
                                                                                                                                                gagagGCCCCTTTTATCATGGTGCATAGGTGCGCTGTGCGCCCGCAAAACCCACACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcact  <  8:147148/148‑1 (MQ=255)
                                                                                                                                                cAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcact  <  5:220118/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACT  >  NC_002947/5102943‑5103235

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: