Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,103,090 A→T 8.4% intergenic (+22/‑168) 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,0900AT 8.4% 235.3 / 8.3 81intergenic (+22/‑168)phhB/PP_5694pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base T (4/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

AATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACT  >  NC_002947/5102950‑5103235
                                                                                                                                            |                                                                                                                                                   
aaTGGGGCAAGGTCGCCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtt                                                                                                                                             >  6:212090/1‑149 (MQ=255)
 aTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACTCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttt                                                                                                                                            <  3:220414/149‑1 (MQ=255)
 aTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttt                                                                                                                                            >  2:166752/1‑149 (MQ=255)
  tGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttcc                                                                                                                                           <  3:250614/149‑3 (MQ=255)
    gggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          <  2:63131/148‑1 (MQ=255)
    gggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          <  2:290743/148‑1 (MQ=255)
    gggCAAGGTCACCGTGACCTGGTGGAGCCACTAGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                         <  5:228010/149‑1 (MQ=255)
      gCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttccttct                                                                                                                                       <  5:360870/149‑7 (MQ=255)
       cAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCAt                                                                                                                                      <  3:326520/149‑1 (MQ=255)
       cAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCAt                                                                                                                                      <  3:133524/149‑1 (MQ=255)
        caGGTCACGGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                                     >  5:152288/2‑149 (MQ=255)
        aaGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCAAGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                                     >  3:330502/1‑149 (MQ=255)
          ggTCACCGTTACCTGGTGGAGCAACTCTATCAAAGGCCTTCACCGCAAAGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGTAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  2:86735/149‑1 (MQ=255)
          ggTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  2:92351/149‑1 (MQ=255)
              aCCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCCGCGGACGCGTTTTTAAGTGGGCat                                                                                                                               >  4:310250/1‑149 (MQ=255)
                  tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAtgt                                                                                                                           <  3:227507/149‑1 (MQ=255)
                    aCCTGGTGGAGCCACTCGATCAAAGGCCTGCACAGCAACGACTTCATCATGTGTGCACGCGCTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgc                                                                                                                         >  1:53212/1‑149 (MQ=255)
                         gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGt                                                                                                                    >  2:82688/1‑149 (MQ=255)
                            gAGCCACTCGATCAAAGGCCTGCACGGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCATAATGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTc                                                                                                                  <  7:224933/148‑1 (MQ=255)
                              gCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAATGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                                  <  5:284906/130‑1 (MQ=255)
                              gCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                                  >  6:284906/1‑130 (MQ=255)
                                cACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTCTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                             >  8:119475/1‑149 (MQ=255)
                                 aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          <  3:117684/119‑1 (MQ=255)
                                 aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          >  4:117684/1‑119 (MQ=255)
                                  cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGc                                                                                                            >  3:302996/1‑148 (MQ=255)
                                  cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTGTTTCATGGGGCATATGTGCTCTGTTCGCCCGc                                                                                                            >  2:79162/1‑148 (MQ=255)
                                     gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttt                                                                                                                                           <  3:165913/114‑1 (MQ=255)
                                     gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttt                                                                                                                                           >  4:165913/1‑114 (MQ=255)
                                        cAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttc                                                                                                                                            <  6:6875/110‑2 (MQ=255)
                                        cAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttc                                                                                                                                            >  5:6875/1‑109 (MQ=255)
                                                 gCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGGTGTTCGCCCGACAAACCGACAAACGAc                                                                                            >  8:224716/1‑149 (MQ=255)
                                                  cACCTCAACGTCTTCATAATGTGTTCACGCACTGACAAGGTGTCTGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGAcc                                                                                           <  7:285009/149‑1 (MQ=255)
                                                  cACCGCAACGTCTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACTCCTGATGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGc                                                                                                                                 <  5:118096/111‑1 (MQ=255)
                                                  cACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGa                                                                                             <  2:354481/147‑1 (MQ=255)
                                                  cACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTATTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGAcc                                                                                           >  1:336182/1‑149 (MQ=255)
                                                  cACCGCAACGACTTCATCATCTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGc                                                                                                                                 >  6:118096/1‑111 (MQ=255)
                                                    ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCgg                                                                                        >  3:325800/1‑149 (MQ=255)
                                                    ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCgg                                                                                        <  3:172109/149‑1 (MQ=255)
                                                    ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCgg                                                                                        >  3:48540/1‑149 (MQ=255)
                                                          cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGaaaa                                                                                   >  6:37045/1‑148 (MQ=255)
                                                              ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                              >  1:56510/1‑148 (MQ=255)
                                                                catcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                            >  4:275907/1‑148 (MQ=255)
                                                                 atcatGTGTGCACGCGCTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaaa                                                                          <  2:53212/149‑1 (MQ=255)
                                                                  tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                         >  4:78240/1‑149 (MQ=255)
                                                                  tcatGTGTGCACGCAATGACAAGGTGGCGGAAACGGCTGAAGGCCGTAAGTAAAGACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                         <  3:310250/149‑1 (MQ=255)
                                                                    atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTg                                                                       >  3:229632/1‑149 (MQ=255)
                                                                    atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGCCATATTTGCGGGGTTCGCCCGCAAAACCGCCAAACGACCGGGAAAAAACCCataatg                                                                       >  6:135257/1‑144 (MQ=255)
                                                                            acgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                  >  8:235377/1‑146 (MQ=255)
                                                                            acgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                  >  1:290204/1‑146 (MQ=255)
                                                                                   gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGAATtgtg                                                        >  6:396572/1‑149 (MQ=255)
                                                                                    aCAAGGTGGCGTATCCGGCTTAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCAAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCGGCAAAAGCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                         <  1:79162/147‑1 (MQ=255)
                                                                                    aCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                         <  2:56510/147‑1 (MQ=255)
                                                                                        ggTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATcc                                                   <  4:330502/149‑1 (MQ=255)
                                                                                         gTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCt                                                  <  2:36083/149‑1 (MQ=255)
                                                                                           ggCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGc                                                <  4:325800/149‑1 (MQ=255)
                                                                                             cGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAATCCAAACTGTCATCCAGACTTGTGTATCCTTccc                                              >  4:311174/1‑149 (MQ=255)
                                                                                             cGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGccc                                              >  8:99534/1‑149 (MQ=255)
                                                                                               gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                              <  8:247341/115‑1 (MQ=255)
                                                                                               gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                              >  7:247341/1‑115 (MQ=255)
                                                                                                     gCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtttt                                      <  5:212090/149‑1 (MQ=255)
                                                                                                      cTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCATGCCCGAGCttttt                                     >  7:280710/1‑149 (MQ=255)
                                                                                                         aaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAATACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGa                                   <  1:82688/149‑1 (MQ=255)
                                                                                                         aaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                  <  4:302996/149‑1 (MQ=255)
                                                                                                          aGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                  <  7:224716/148‑1 (MQ=255)
                                                                                                          aGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTCTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                  <  7:119475/148‑1 (MQ=255)
                                                                                                            gCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGc                                >  7:290485/1‑148 (MQ=255)
                                                                                                             ccGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                               <  5:37045/148‑1 (MQ=255)
                                                                                                             ccGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTa                              <  5:396572/149‑1 (MQ=255)
                                                                                                              cGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTAc                             <  2:290204/149‑1 (MQ=255)
                                                                                                              cGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTAc                             <  7:99534/149‑1 (MQ=255)
                                                                                                               ggAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACa                            <  3:275907/149‑1 (MQ=255)
                                                                                                                      aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGc                                >  2:21936/1‑138 (MQ=255)
                                                                                                                      aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGc                                <  1:21936/138‑1 (MQ=255)
                                                                                                                            tgaaaaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  3:47121/7‑149 (MQ=255)
                                                                                                                            aGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  6:303387/1‑149 (MQ=255)
                                                                                                                             gaaggaagGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg               >  3:159605/1‑148 (MQ=255)
                                                                                                                                aaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGa           >  4:374909/3‑149 (MQ=255)
                                                                                                                                  aaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACt         >  5:160492/1‑149 (MQ=255)
                                                                                                                                     gggCCAGAGGCCCCGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACGTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACttt       >  8:134262/1‑148 (MQ=255)
                                                                                                                                        ccAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCTTGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtca    >  1:325545/1‑149 (MQ=255)
                                                                                                                                         cAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcact  >  6:10438/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
AATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACT  >  NC_002947/5102950‑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: