Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,195 T→G 5.3% intergenic (+23/‑142) pheA → / → PP_1770 bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,1950TG 5.3% 237.5 / 2.9 76intergenic (+23/‑142)pheA/PP_1770bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (32/40);  new base G (2/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

ACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGG  >  NC_002947/1976048‑1976340
                                                                                                                                                   |                                                                                                                                                 
aCGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTg                                                                                                                                                  <  8:217514/149‑1 (MQ=255)
aCGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTg                                                                                                                                                  >  7:220519/1‑149 (MQ=255)
    gtTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGc                                                                                                                                              >  4:177151/1‑149 (MQ=255)
     ttcttcATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGcc                                                                                                                                             >  4:205766/1‑149 (MQ=255)
       cttcATCGACTTCGTCCGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCa                                                                                                                                                >  5:315239/1‑144 (MQ=255)
       cttcATCGACTTCGTCCGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCa                                                                                                                                                <  6:315239/144‑1 (MQ=255)
          cATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTcg                                                                                                                                        >  6:236827/1‑149 (MQ=255)
           aTCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTcgc                                                                                                                                       <  3:241388/149‑1 (MQ=255)
            tCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTcgcg                                                                                                                                      >  8:196477/1‑149 (MQ=255)
             cGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGTAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGGCTTGCAGCCCTTCGCgg                                                                                                                                     >  5:184672/1‑149 (MQ=255)
                cTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGct                                                                                                                                  <  8:276072/149‑2 (MQ=255)
                       ggCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGc                                                                                                                           >  1:127283/1‑149 (MQ=255)
                        gCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCt                                                                                                                          >  5:316753/1‑149 (MQ=255)
                         ccaccaCCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTc                                                                                                                         <  6:8693/149‑1 (MQ=255)
                         ccaccaCCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCAAGGCAGccccaag                                                                                                                                        >  8:264127/1‑130 (MQ=255)
                         ccaccaCCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCAAGGCAGccccaag                                                                                                                                        <  7:264127/134‑5 (MQ=255)
                          caccacCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTcc                                                                                                                        >  8:201200/1‑149 (MQ=255)
                              accGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCcaca                                                                                                                    <  2:127283/149‑1 (MQ=255)
                              accGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCcaca                                                                                                                    <  3:182666/149‑1 (MQ=255)
                              accGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCcaca                                                                                                                    <  8:243379/149‑1 (MQ=255)
                                  cgACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGa                                                                                                                <  5:336185/149‑1 (MQ=255)
                                    aCCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAg                                                                                                               >  6:192283/1‑148 (MQ=255)
                                      ccGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAt                                                                                                            <  4:22256/149‑1 (MQ=255)
                                        gCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGATGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAtgt                                                                                                          <  5:345242/149‑1 (MQ=255)
                                        gCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAtgt                                                                                                          <  7:196477/149‑1 (MQ=255)
                                             tCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGACGGGGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTc                                                                                                                                         <  5:268808/113‑1 (MQ=255)
                                             tCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTc                                                                                                                                         >  6:268808/1‑113 (MQ=255)
                                             tCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGc                                                                                                     <  8:220519/149‑1 (MQ=255)
                                                 ggCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTg                                                                                                                 <  1:65367/133‑1 (MQ=255)
                                                 ggCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTg                                                                                                                 >  2:65367/1‑133 (MQ=255)
                                                           gAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTgg                                                                                         >  3:205441/1‑147 (MQ=255)
                                                              cAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCg                                                                                    <  6:46312/149‑1 (MQ=255)
                                                              cAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCg                                                                                    <  5:236827/149‑1 (MQ=255)
                                                                  tCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGc                                                                                 <  3:205766/148‑1 (MQ=255)
                                                                      ccAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGcc                                                                            >  3:249201/1‑149 (MQ=255)
                                                                             gCTGTGGCGCTGAAGGTGCTGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAt                                                                                                            <  1:25/109‑1 (MQ=255)
                                                                             gCTGTGGCGCTGAAGGTGCTGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAt                                                                                                            >  2:25/1‑109 (MQ=255)
                                                                               tgtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAggg                                                                   >  2:60088/1‑149 (MQ=255)
                                                                               tgtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCCCGAAggg                                                                   >  8:202980/1‑149 (MQ=255)
                                                                                gtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGc                                                                                 <  5:32484/134‑1 (MQ=255)
                                                                                gtgGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGc                                                                                 >  6:32484/1‑134 (MQ=255)
                                                                                   gcgcTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGgccgc                                                               <  4:205441/149‑1 (MQ=255)
                                                                                         aaGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGccgcgccgcg                                                         >  6:17634/1‑149 (MQ=255)
                                                                                           ggTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGc                                                       <  4:249201/149‑1 (MQ=255)
                                                                                              ggtGGGCTCTTAGCCGAGGGCGGTGCGTTGATTCAGGGGCCTTGGGGCTGGCTTGCAGCCCTTCGCGCCCAGGCCCGCTACCACTGTGAGCAGGTGAAGCGGAGCGCTGTGGGAGCGGGCGAGCACGGGAAGGGCCGCCCCGCGGcccc                                                    >  2:73451/3‑149 (MQ=255)
                                                                                                 gggCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAt                                                 <  5:192283/149‑1 (MQ=255)
                                                                                                  ggCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCGGCCTTGCAGCCCTTCGCGGGCATCCCCGCTCCCACAGTGAGCATGTGATCCCGACCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCaa                                                  <  1:73451/147‑1 (MQ=255)
                                                                                                   gCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGc                                                       >  4:170094/1‑141 (MQ=255)
                                                                                                   gCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGc                                                       <  3:170094/141‑1 (MQ=255)
                                                                                                     tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGccc                                                     >  3:135377/1‑141 (MQ=255)
                                                                                                     tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGccc                                                     <  4:135377/141‑1 (MQ=255)
                                                                                                     tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                             >  8:150658/1‑149 (MQ=255)
                                                                                                     tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                             >  2:52210/1‑149 (MQ=255)
                                                                                                       ttATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                             <  7:201200/147‑1 (MQ=255)
                                                                                                         aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTca                                         >  2:132117/1‑149 (MQ=255)
                                                                                                         aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTc                                          <  6:168531/148‑1 (MQ=255)
                                                                                                         aTCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTc                                          >  5:168531/1‑148 (MQ=255)
                                                                                                                       cTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAt                           >  4:152830/1‑149 (MQ=255)
                                                                                                                         ttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAtt                          >  6:198034/1‑148 (MQ=255)
                                                                                                                            aTTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                       <  5:17634/148‑1 (MQ=255)
                                                                                                                            aTTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                       <  1:52210/148‑1 (MQ=255)
                                                                                                                            aTTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                       <  6:184672/148‑1 (MQ=255)
                                                                                                                                 tccgttctGGGGCTGCCTTGCAGCCCTGCGCGGTCATGCCCTCTCCCACAGTGAGCCTGTGAGTCCTATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACa                  <  7:294557/141‑1 (MQ=255)
                                                                                                                                 ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                       >  2:239664/1‑143 (MQ=255)
                                                                                                                                 ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                       <  1:239664/143‑1 (MQ=255)
                                                                                                                                  gggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                <  7:202980/149‑1 (MQ=255)
                                                                                                                                     ccTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGtt             <  1:60088/149‑1 (MQ=255)
                                                                                                                                       gaagggCTGCAAGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACa                  >  6:319933/4‑142 (MQ=255)
                                                                                                                                       gaagggCTGCAAGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACa                  <  5:319933/139‑1 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGc          <  3:177151/148‑1 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGAGCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCg           <  6:355233/147‑1 (MQ=255)
                                                                                                                                           ggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                             <  3:291101/111‑1 (MQ=255)
                                                                                                                                           ggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCAGGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                             >  4:291101/1‑111 (MQ=255)
                                                                                                                                              tGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGt    >  8:142075/1‑149 (MQ=255)
                                                                                                                                                ccTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGt    <  5:198034/147‑1 (MQ=255)
                                                                                                                                                ccTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTgg  <  5:322353/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGCTGAAGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGG  >  NC_002947/1976048‑1976340

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: