Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,329,743 G→T D261Y (GAT→TAT)  dacB → D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,7430GT97.4% 1318.7 / 20.3 377D261Y (GAT→TAT) dacBD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base T (194/173);  total (198/179)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 379 total aligned reads displayed.

TGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGA  >  NC_000913/3329595‑3329890
                                                                                                                                                    |                                                                                                                                                   
       gcgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaa                                                                                                                                            <  3:849551/151‑1 (MQ=255)
       gcgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaa                                                                                                                                            <  1:977910/151‑1 (MQ=255)
        cgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaa                                                                                                                                           <  3:1365861/151‑1 (MQ=255)
           aaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAc                                                                                                                                         >  7:943299/1‑150 (MQ=255)
            aTACTTCCACCTGGATGGTGGGCCAGGCGACCGGACCCGCATTCCGCTGACGGGACGCCTGCCTCAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAg                                                                                                                                       <  2:1460867/151‑1 (MQ=255)
             tACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAgg                                                                                                                                      >  4:766507/1‑151 (MQ=255)
              aCTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAgg                                                                                                                                      <  1:24405/150‑1 (MQ=255)
               cTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAACACAGGCg                                                                                                                                    >  4:1310525/1‑151 (MQ=255)
                tGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                   >  5:1425436/1‑151 (MQ=255)
                 gCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                   <  7:852166/150‑1 (MQ=255)
                 gCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCg                                                                                                                                    <  4:1270905/149‑1 (MQ=255)
                    aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTAt                                                                                                                               >  1:1517924/1‑151 (MQ=255)
                    aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTAt                                                                                                                               >  8:1696426/1‑151 (MQ=255)
                        ggATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCa                                                                                                                             <  4:912611/149‑1 (MQ=255)
                         gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCt                                                                                                                          <  8:1074603/151‑1 (MQ=255)
                           tgtgGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGCAATATCAGTTAATACAGGAGGGTATCACCTgg                                                                                                                        >  3:1440347/1‑151 (MQ=255)
                                tgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCgg                                                                                                                   >  3:705083/1‑151 (MQ=255)
                                  ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGaa                                                                                                                 <  7:322510/151‑1 (MQ=255)
                                   cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAAGAGGCGGGTATCACCTGGAGCGGAac                                                                                                                <  1:134884/151‑1 (MQ=255)
                                    aGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAac                                                                                                                >  5:662528/1‑150 (MQ=255)
                                     ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAacac                                                                                                              <  2:693678/151‑1 (MQ=255)
                                      gCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACAct                                                                                                             <  1:632149/151‑1 (MQ=255)
                                        gACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                           <  4:537999/151‑1 (MQ=255)
                                        gACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                           >  5:1663210/1‑151 (MQ=255)
                                        gACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                           >  3:182931/1‑151 (MQ=255)
                                        gACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                           >  7:459839/1‑151 (MQ=255)
                                          ccTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgctg                                                                                                         <  2:1069211/151‑1 (MQ=255)
                                           cTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACtgctgc                                                                                                        >  1:1259595/1‑151 (MQ=255)
                                           cTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGACAGCTATGCCGGGGCAATTCTGAAATATGAGTTAACACTGGCGGGTATCACCTGGCGCGGTACAGtgctgc                                                                                                        >  1:1052210/1‑151 (MQ=255)
                                            tGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTgcg                                                                                                       >  8:114803/1‑151 (MQ=255)
                                             gAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTgcg                                                                                                       <  6:1002434/150‑1 (MQ=255)
                                              aaCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGcc                                                                                                     >  6:1671332/1‑151 (MQ=255)
                                                     ttagtcTGACTGGATGACTGCCACAACGTTCTGAGCCGCTCCCGTTTGCTTTTGCCGTGCAGGATGTAGCCAGCTATGCCGGTGGAATTCTGAAATATGAGTTAAAAGAGCCGGTAATCTCCTGGAGCGGAACACTGCTGCTCCAGACt                                                                                                <  2:83451/144‑1 (MQ=255)
                                                     ttACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTc                                                                                               >  3:365382/1‑150 (MQ=255)
                                                     ttACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTc                                                                                               <  8:775483/150‑1 (MQ=255)
                                                      tACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCa                                                                                              <  7:1270910/150‑1 (MQ=255)
                                                       aCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCCGAACACTGCTGCGCCAGACTCAgg                                                                                            >  7:1445895/1‑151 (MQ=255)
                                                        cGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGt                                                                                           >  8:1469424/1‑151 (MQ=255)
                                                          cTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGtt                                                                                          >  6:666271/1‑150 (MQ=255)
                                                            gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAAc                                                                                       >  5:1152430/1‑151 (MQ=255)
                                                                ggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAc                                                                                   <  8:929629/151‑1 (MQ=255)
                                                                 gATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAcc                                                                                  <  3:646753/151‑1 (MQ=255)
                                                                  aTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCt                                                                                 <  6:357572/151‑1 (MQ=255)
                                                                  aTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCt                                                                                 <  1:269143/151‑1 (MQ=255)
                                                                  aTGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCt                                                                                 <  3:1602835/151‑1 (MQ=255)
                                                                   tgcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTg                                                                                >  7:1053001/1‑151 (MQ=255)
                                                                   tgcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTg                                                                                >  8:475503/1‑151 (MQ=255)
                                                                    gcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTgg                                                                               <  2:1241779/151‑1 (MQ=255)
                                                                      ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGACCACTGCTGCGCCAGACTCAGGTTAACGAACCTGGa                                                                              >  1:578876/1‑150 (MQ=255)
                                                                        gccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                            <  6:1512251/150‑1 (MQ=255)
                                                                        gccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACg                                                                           <  6:994094/151‑1 (MQ=255)
                                                                          cacaACGTTCTGAGCCGCTCCCGTTGGCTTTTGGCGTGCAGGATGGAGCCAGCTATGCCGGGGCAATTGTGAAAGATGAGTTGAAACAGGGGGGATTCAGCTGGAGCGGAAGACTGCTGCGCCAGCCTGAGGGTAAGGAAGATAGAACGGt                                                                         >  2:288982/1‑151 (MQ=255)
                                                                            caACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGt                                                                         >  8:1238022/1‑149 (MQ=255)
                                                                              aCGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAACCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGtt                                                                     <  1:1590772/151‑1 (MQ=255)
                                                                                    tGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTTGTTGCCAGt                                                               >  3:1156100/1‑151 (MQ=255)
                                                                                     gAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAg                                                                <  4:65461/149‑1 (MQ=255)
                                                                                        ccGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAAc                                                           <  1:470167/151‑1 (MQ=255)
                                                                                        ccGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGAGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAAc                                                           >  6:457015/1‑151 (MQ=255)
                                                                                         cGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACa                                                          <  4:223141/151‑1 (MQ=255)
                                                                                              ccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCgg                                                     >  7:531739/1‑151 (MQ=255)
                                                                                                  tGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGCCTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGcccc                                                 >  7:771620/1‑151 (MQ=255)
                                                                                                   ggCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCg                                                <  5:279469/151‑1 (MQ=255)
                                                                                                   ggCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCAGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGcccc                                                 <  4:1043588/150‑1 (MQ=255)
                                                                                                       tttGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCt                                              <  4:877608/149‑1 (MQ=255)
                                                                                                       tttGCCGTGCAGGAGGGAGCCAGCTATGCCGTTGAAATTCTGAAATATGAGTTATATGTGGCGTTTATCTCCTGGAGCGTAACACTGCTGCGGCAGACTCTGGTTAACGACCCTGGCACGGTAGTTGCCAGTAAACAGTCGGCCCCGCt                                              <  1:785160/149‑1 (MQ=255)
                                                                                                         tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAGCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCa                                           <  6:554885/150‑1 (MQ=255)
                                                                                                         tGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCa                                           <  1:1114056/150‑1 (MQ=255)
                                                                                                          gCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCAc                                          <  4:1517074/150‑1 (MQ=255)
                                                                                                          gCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCAc                                          <  5:1534628/150‑1 (MQ=255)
                                                                                                            cGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                       <  2:657163/151‑1 (MQ=255)
                                                                                                             gTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCgg                                                     >  7:1136804/1‑136 (MQ=255)
                                                                                                              tGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                     >  4:301488/1‑151 (MQ=255)
                                                                                                               gCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTg                                    <  5:1209247/151‑1 (MQ=255)
                                                                                                                 cggATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                  <  8:63016/150‑1 (MQ=255)
                                                                                                                 aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                  >  3:550101/1‑151 (MQ=255)
                                                                                                                 aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                  >  4:177569/1‑151 (MQ=255)
                                                                                                                 aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCt                                  >  5:718882/1‑151 (MQ=255)
                                                                                                                    aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAGGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTa                                >  5:573238/1‑150 (MQ=255)
                                                                                                                      gtaGCCAGCTATGCAGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAg                              <  2:1642439/148‑1 (MQ=255)
                                                                                                                         gCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTa                          <  7:602944/151‑1 (MQ=255)
                                                                                                                           cAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATg                        >  5:309370/1‑151 (MQ=255)
                                                                                                                            aGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGc                       <  2:193465/151‑1 (MQ=255)
                                                                                                                             gCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCt                      <  8:41926/151‑1 (MQ=255)
                                                                                                                              cTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTg                     >  5:1013649/1‑151 (MQ=255)
                                                                                                                                aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGCTCTGCTTAAGATTATGCTGa                    >  4:1563023/1‑150 (MQ=255)
                                                                                                                                aTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTg                     <  1:927301/149‑1 (MQ=255)
                                                                                                                                   ccGGGGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCACCGCTGCACGATCTGCTTAAGATTATGCTGaaaaa                <  5:850423/151‑1 (MQ=255)
                                                                                                                                    cGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaaaa                 >  6:1248997/1‑149 (MQ=255)
                                                                                                                                    cGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAg               <  6:496618/151‑1 (MQ=255)
                                                                                                                                      gTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTc             >  1:366256/1‑151 (MQ=255)
                                                                                                                                       tGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCg            >  6:289573/1‑151 (MQ=255)
                                                                                                                                          aaTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGac         >  4:1403873/1‑151 (MQ=255)
                                                                                                                                           aTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGaca        <  7:400278/151‑1 (MQ=255)
                                                                                                                                           aTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGac         >  5:1288985/1‑150 (MQ=255)
                                                                                                                                               tGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAt    >  7:1582143/1‑151 (MQ=255)
                                                                                                                                               tGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAt    >  2:505607/1‑151 (MQ=255)
                                                                                                                                                gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAt    <  2:1387247/150‑1 (MQ=255)
                                                                                                                                                gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATg   <  7:210414/151‑1 (MQ=255)
                                                                                                                                                gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATg   <  6:815467/151‑1 (MQ=255)
                                                                                                                                                 aaaTATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATg   >  7:1627477/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGA  >  NC_000913/3329595‑3329890

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: