Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,120,470 T→A 6.4% intergenic (‑212/‑90) pepA ← / → PP_0982 cytosol aminopeptidase/YjgP/YjgQ family permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,120,4700TA 6.4% 290.2 / 6.8 94intergenic (‑212/‑90)pepA/PP_0982cytosol aminopeptidase/YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base T (36/52);  new base A (2/4);  total (38/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

CCGCGTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCCGTCAGCG  >  NC_002947/1120323‑1120617
                                                                                                                                                   |                                                                                                                                                   
ccGCGTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCttt                                                                                                                                                    <  1:86300/149‑1 (MQ=255)
 cgcgTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTg                                                                                                                                                   <  4:207050/149‑1 (MQ=255)
  gcgTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGTCAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTgg                                                                                                                                                  <  8:42924/149‑1 (MQ=255)
  gcgTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTgg                                                                                                                                                  <  7:248878/149‑1 (MQ=255)
     tcgtcgCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCg                                                                                                                                                >  1:178715/1‑148 (MQ=255)
       gtcgCGCACCAGCGGAATAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCGCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTGCGCCTTTGGCGAAc                                                                                                                                             <  5:36919/149‑1 (MQ=255)
       gtcgCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAAc                                                                                                                                             <  8:90098/149‑1 (MQ=255)
          gcgcACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCgg                                                                                                                                          >  4:343930/1‑149 (MQ=255)
          gcgcACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCg                                                                                                                                           <  3:117666/148‑1 (MQ=255)
            gcACCAGCGGAAAAGCTGCAAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctcca                                                                                                                                        <  8:137289/149‑5 (MQ=255)
             cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCAt                                                                                                                                       <  2:114659/149‑1 (MQ=255)
             cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctccat                                                                                                                                       <  6:133634/149‑6 (MQ=255)
             cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctccat                                                                                                                                       <  6:70655/149‑6 (MQ=255)
              aCCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCAtt                                                                                                                                      >  3:92679/1‑149 (MQ=255)
               ccAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTg                                                                                                                                     <  6:163669/149‑1 (MQ=255)
                cAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTgg                                                                                                                                    >  8:96101/1‑149 (MQ=255)
                cAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTg                                                                                                                                     >  6:329165/1‑148 (MQ=255)
                cAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTg                                                                                                                                     <  5:329165/148‑1 (MQ=255)
                    ggAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                <  3:349223/149‑1 (MQ=255)
                    ggAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                <  1:148940/149‑1 (MQ=255)
                    ggAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                >  2:296090/1‑149 (MQ=255)
                    ggAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGGCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                <  5:232450/149‑1 (MQ=255)
                           cTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTa                                                                                                                         >  6:110766/1‑149 (MQ=255)
                            tGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTAc                                                                                                                        >  6:219374/1‑149 (MQ=255)
                             gCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTAc                                                                                                                        >  3:73932/1‑148 (MQ=255)
                                aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                                    <  6:373579/149‑1 (MQ=255)
                                aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGCCGGTTCGCCAAAGGCGAAcctccat                                                                                                                                       <  1:391208/130‑6 (MQ=255)
                                aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGCCGGTTCGCCAAAGGCGAAcctccat                                                                                                                                       >  2:391208/1‑125 (MQ=255)
                                 gTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                                    >  1:316664/1‑148 (MQ=255)
                                  tttGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGt                                                                                                                   >  3:207017/1‑148 (MQ=255)
                                      aGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCAtt                                                                                                               <  4:207017/148‑1 (MQ=255)
                                      aGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCAtt                                                                                                               <  2:178715/148‑1 (MQ=255)
                                       gCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGt                                                                                                             <  8:212724/149‑1 (MQ=255)
                                        ccTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGtt                                                                                                            <  8:161606/149‑1 (MQ=255)
                                         cTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGttt                                                                                                           >  3:128058/1‑149 (MQ=255)
                                            tgGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGgc                                                                                                        <  4:73932/149‑1 (MQ=255)
                                              gCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGgcgc                                                                                                      >  7:314552/1‑149 (MQ=255)
                                                aaGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCa                                                                                                    <  5:219374/149‑1 (MQ=255)
                                                       cTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATg                                                                                             <  2:316664/149‑1 (MQ=255)
                                                       cTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGGTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCAt                                                                                              <  7:226213/148‑1 (MQ=255)
                                                       cTGACAACCACGATGTGCGCATTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACATGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCAt                                                                                              <  7:80224/148‑1 (MQ=255)
                                                        tGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                <  8:309868/113‑1 (MQ=255)
                                                        tGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                >  7:309868/1‑113 (MQ=255)
                                                         gACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCAGTGGTCTTGGAGGTACTAGGTCATTGTTTGGCGCCCTTGCATgg                                                                                            >  8:327953/1‑148 (MQ=255)
                                                          aCAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                >  4:247356/1‑111 (MQ=255)
                                                          aCAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                                <  3:247356/111‑1 (MQ=255)
                                                            aaCCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATgg                                                                                            <  7:21908/145‑1 (MQ=255)
                                                            aaCCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATgg                                                                                            >  8:21908/1‑145 (MQ=255)
                                                                 cGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGcc                                                                                   <  4:92679/149‑1 (MQ=255)
                                                                  gATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGccc                                                                                  <  2:299291/149‑1 (MQ=255)
                                                                  gATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGccc                                                                                  <  7:8371/149‑1 (MQ=255)
                                                                   aTGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCt                                                                                 <  1:296090/149‑1 (MQ=255)
                                                                    tgtgCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTg                                                                                <  5:110766/149‑1 (MQ=255)
                                                                     gtgCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATgg                                                                                            <  2:299782/136‑1 (MQ=255)
                                                                     gtgCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATgg                                                                                            >  1:299782/1‑136 (MQ=255)
                                                                       gcgcCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGAc                                                                              >  4:334321/1‑148 (MQ=255)
                                                                        cgcCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACaa                                                                            <  3:334321/149‑1 (MQ=255)
                                                                         gcCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAc                                                                           >  3:349521/1‑149 (MQ=255)
                                                                         gcCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAc                                                                           >  2:307003/1‑149 (MQ=255)
                                                                          ccTTCGTACGCGGCCAAGTTACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAcc                                                                          <  7:75191/149‑1 (MQ=255)
                                                                          ccTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAcc                                                                          <  7:218371/149‑1 (MQ=255)
                                                                           cTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCAt                                                                                                   <  3:234642/123‑1 (MQ=255)
                                                                           cTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCAt                                                                                                   >  4:234642/1‑123 (MQ=255)
                                                                           cTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAcc                                                                          >  5:73315/1‑148 (MQ=255)
                                                                           cTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAcc                                                                          <  6:73315/148‑1 (MQ=255)
                                                                            ttCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAcc                                                                          >  4:344949/1‑147 (MQ=255)
                                                                              cGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCa                                                                        <  3:171859/147‑1 (MQ=255)
                                                                              cGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCa                                                                        <  7:240558/147‑1 (MQ=255)
                                                                              cGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAg                                                                       <  7:96101/148‑1 (MQ=255)
                                                                               gTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGa                                                                     <  1:175702/149‑1 (MQ=255)
                                                                               gTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGa                                                                     <  5:113736/149‑1 (MQ=255)
                                                                                         aaGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAg                                                                       >  3:158782/1‑137 (MQ=255)
                                                                                         aaGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAg                                                                       <  4:158782/137‑1 (MQ=255)
                                                                                         aaGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGtt                                                           >  2:391329/1‑149 (MQ=255)
                                                                                            tGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGa                                                        >  6:129507/1‑149 (MQ=255)
                                                                                            tGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGa                                                        <  4:128058/149‑1 (MQ=255)
                                                                                                aGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGt                                                    >  5:315027/1‑149 (MQ=255)
                                                                                                       gTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGt                                             <  7:327953/149‑1 (MQ=255)
                                                                                                               caGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTAATAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTc                                     >  5:283126/1‑149 (MQ=255)
                                                                                                                  gATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCcg                                  >  7:261679/1‑149 (MQ=255)
                                                                                                                   aTAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCcgc                                 >  3:267435/1‑149 (MQ=255)
                                                                                                                          gcgcACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTcc                          <  6:315027/149‑1 (MQ=255)
                                                                                                                            gcACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTg                        <  8:314552/149‑1 (MQ=255)
                                                                                                                                 ttACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGt                                             >  8:124430/1‑123 (MQ=255)
                                                                                                                                 ttACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGt                                             <  7:124430/123‑1 (MQ=255)
                                                                                                                                   aCATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCAcc                  <  1:391329/148‑1 (MQ=255)
                                                                                                                                      tGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTAGTCACCTTGa              <  5:129507/149‑1 (MQ=255)
                                                                                                                                       ggAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAg             >  8:27331/1‑149 (MQ=255)
                                                                                                                                       ggAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAg             >  1:261829/1‑149 (MQ=255)
                                                                                                                                           gTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGcc         <  1:307003/149‑1 (MQ=255)
                                                                                                                                           gTTCGCCAAAGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGcc         >  7:34943/1‑149 (MQ=255)
                                                                                                                                               gCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCCGt       >  7:339307/1‑147 (MQ=255)
                                                                                                                                               gCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCCGTCa     <  3:344949/149‑1 (MQ=255)
                                                                                                                                                 cTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGa              <  7:366965/138‑1 (MQ=255)
                                                                                                                                                 cTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGa              >  8:366965/1‑138 (MQ=255)
                                                                                                                                                  tttGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGCCACCTTGAGCGCCGTCAgcg  >  3:199201/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGCGTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCCGTCAGCG  >  NC_002947/1120323‑1120617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: