Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,103,090 A→T 7.2% 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 7.2% 191.9 / 3.2 64intergenic (+22/‑168)phhB/PP_5694pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base T (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5102947‑5103237
                                                                                                                                               |                                                                                                                                                     
ctgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGccc                                                                                                                                                >  4:298742/1‑149 (MQ=255)
   aaTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtt                                                                                                                                             >  4:979/1‑149 (MQ=255)
        ggCAAGGTCACCGTGACCTGGTGGAGTCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                         >  2:550/1‑148 (MQ=255)
        ggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttccttc                                                                                                                                        <  1:148799/149‑6 (MQ=255)
         gCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGccccatgaaaa                                                                                                                                        >  5:306249/1‑141 (MQ=255)
            aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgg                                                                                                                                    <  3:979/149‑1 (MQ=255)
               gcACCGTGACCTTGTGGAGCCACTCCATCTCACGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGc                                                                                                                                 <  8:160753/148‑1 (MQ=255)
                 aCCGTGACCTGGTGGAGCCACTCTATCAAAGCCCTGCACCGCAACGACATCATCATGTGTGCACGCACTGACAATGTGGCGGAAACGGCTGAAGGCCGTAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCat                                                                                                                               <  4:300378/149‑1 (MQ=255)
                 aCCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCat                                                                                                                               <  6:61249/149‑1 (MQ=255)
                    gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttcctt                                                                                                                                         <  6:111290/136‑5 (MQ=255)
                    gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttcctt                                                                                                                                         >  5:111290/1‑132 (MQ=255)
                        ccTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcg                                                                                                                        <  3:298742/149‑1 (MQ=255)
                            gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGt                                                                                                                    <  7:153709/149‑1 (MQ=255)
                            gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGt                                                                                                                    >  3:138691/1‑149 (MQ=255)
                               gAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTc                                                                                                                  <  3:319194/148‑1 (MQ=255)
                               gAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCg                                                                                                                 >  2:293651/1‑149 (MQ=255)
                                  ccACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          <  3:194066/121‑1 (MQ=255)
                                  ccACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          <  3:198931/121‑1 (MQ=255)
                                  ccACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          >  4:194066/1‑121 (MQ=255)
                                  ccACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                          >  4:198931/1‑121 (MQ=255)
                                    aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGc                                                                                                            <  8:29869/149‑1 (MQ=255)
                                         aTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAc                                                                                                       <  4:138691/149‑1 (MQ=255)
                                                gCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaaa                                                                                                >  4:95374/1‑149 (MQ=255)
                                                 ccTTCACAGCAACGACTTCATCTTGTGTGCACGCACTGACATGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCTGATGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaaa                                                                                                <  4:238049/148‑1 (MQ=255)
                                                    gCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGAc                                                                                            >  2:184420/1‑149 (MQ=255)
                                                     cACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGAGCGCTGTTCGCCCGCAAAACCGACAAACGAcc                                                                                           <  2:299275/149‑1 (MQ=255)
                                                            aCGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGc                                                                                                                                                    >  6:214746/1‑86 (MQ=255)
                                                                 ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                                >  7:211858/1‑148 (MQ=255)
                                                                   catcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              >  6:184664/1‑148 (MQ=255)
                                                                     tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                           >  2:338338/1‑149 (MQ=255)
                                                                     taatttgtgCAAGCACTGACAATGTGGCGGAAACGGCTTAACGCCGGAATTAAACACCATAAGGAATGGGCCATTTGCTCACCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaa                                                                             <  1:293651/144‑1 (MQ=255)
                                                                      catGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACt                                                                          >  1:308540/1‑149 (MQ=255)
                                                                          tgtgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCa                                                                      >  7:76936/1‑149 (MQ=255)
                                                                                cgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGa                                                                <  8:41778/149‑1 (MQ=255)
                                                                                    cTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtg                                                            >  8:160068/1‑149 (MQ=255)
                                                                                     tGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                           <  8:215570/149‑1 (MQ=255)
                                                                                     tGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                           <  3:95374/149‑1 (MQ=255)
                                                                                                cctatctggCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGcc                                                 <  7:58951/141‑1 (MQ=255)
                                                                                                     aCGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCt                                           >  2:267130/1‑149 (MQ=255)
                                                                                                        gCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGcc                                                 >  8:58951/1‑140 (MQ=255)
                                                                                                          tGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCCCAAAACCGACAAACCCCCGGCAAAAAAACCAACGTGTCATCCAGACTAGTGTAACCTCCCCACCCttttt                                       >  7:337919/1‑148 (MQ=255)
                                                                                                             aGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                    <  1:184420/148‑1 (MQ=255)
                                                                                                                 cGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                      <  8:60515/126‑1 (MQ=255)
                                                                                                                 cGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                      >  7:60515/1‑126 (MQ=255)
                                                                                                                       tAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcgc                         >  7:274721/1‑149 (MQ=255)
                                                                                                                         aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGTGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACgcgc                       <  8:274721/149‑1 (MQ=255)
                                                                                                                         aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACgcgc                       >  2:2911/1‑149 (MQ=255)
                                                                                                                           cacCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCt                     >  1:214473/1‑149 (MQ=255)
                                                                                                                            acCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTg                    >  8:41710/1‑149 (MQ=255)
                                                                                                                            acCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTg                    >  2:69095/1‑149 (MQ=255)
                                                                                                                              cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATGCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                  <  1:550/149‑1 (MQ=255)
                                                                                                                              cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                  <  6:67639/149‑1 (MQ=255)
                                                                                                                              cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                  >  5:67639/1‑149 (MQ=255)
                                                                                                                              cAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc                  <  8:211858/149‑1 (MQ=255)
                                                                                                                               tgaaaaaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg                 >  7:269906/7‑149 (MQ=255)
                                                                                                                               aGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg                 <  8:76936/149‑1 (MQ=255)
                                                                                                                                  aggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCa               <  2:214473/148‑1 (MQ=255)
                                                                                                                                    gaaggGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATACAGACTTGTGTATCCTGCCCGCGCATTTTGAAGCTACAAACGCGCCTGGCGCAGAc            >  7:292795/1‑149 (MQ=255)
                                                                                                                                    gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGAc            <  2:308540/149‑1 (MQ=255)
                                                                                                                                      aggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGCCGCAGACtt          >  2:136119/1‑149 (MQ=255)
                                                                                                                                       ggggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                                >  5:154563/1‑78 (MQ=255)
                                                                                                                                       ggggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAAc                                                                                <  6:154563/78‑1 (MQ=255)
                                                                                                                                          gCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTAc                               >  7:70592/1‑124 (MQ=255)
                                                                                                                                          gCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTAc                               <  8:70592/124‑1 (MQ=255)
                                                                                                                                             agagGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcact    >  5:75779/1‑148 (MQ=255)
                                                                                                                                              gagGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                      >  3:88945/1‑97 (MQ=255)
                                                                                                                                              gagGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                      <  4:88945/97‑1 (MQ=255)
                                                                                                                                              gagGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcactca  >  7:250081/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5102947‑5103237

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: