Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,105 C→T 6.6% A662V (GCA→GTA)  ydhK → putative efflux protein (PET) component of YdhJK efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,1050CT 6.6% 193.3 / 2.7 60A662V (GCA→GTA) ydhKputative efflux protein (PET) component of YdhJK efflux pump
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCC  >  NC_000913/1723961‑1724248
                                                                                                                                                |                                                                                                                                               
gCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGcacc                                                                                                                                              <  3:145704/148‑1 (MQ=255)
       gCAAAAATCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGc                                                                                                                                         >  3:136900/1‑146 (MQ=255)
         aaaaaTCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAgg                                                                                                                                     <  8:135825/148‑1 (MQ=255)
           aaaTCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTa                                                                                                                                   >  4:99942/1‑148 (MQ=255)
           aaaTCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTa                                                                                                                                   >  3:234983/1‑148 (MQ=255)
           aaaTCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTa                                                                                                                                   >  1:232474/1‑148 (MQ=255)
               gacttggggcCAGACTTGAAGATTTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCTCGCTGCGGGCTGTGTGCTGGCGGGTATTGTCTGGCGGCTGTAGTGCTCGCTTTGGCTACTTGAGCAAGCTCGACCGCAAGGTACGc                                                                                                                                <  2:49112/140‑1 (MQ=255)
                   ggcggcCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGcc                                                                                                                           >  5:343475/1‑148 (MQ=255)
                        ccACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTctta                                                                                                                      <  1:158998/148‑1 (MQ=255)
                        ccACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTctta                                                                                                                      <  2:232474/148‑1 (MQ=255)
                        ccACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTctta                                                                                                                      <  4:136900/148‑1 (MQ=255)
                          acacTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaat                                                                                                                    <  7:19719/148‑1 (MQ=255)
                           cacTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCaa                                                                                                                                       >  2:24470/1‑128 (MQ=255)
                           cacTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCaa                                                                                                                                       <  1:24470/128‑1 (MQ=255)
                            acTAGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCCACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaatta                                                                                                                  <  6:356205/148‑1 (MQ=255)
                             cTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaattac                                                                                                                 >  8:233048/1‑148 (MQ=255)
                              ttGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGagcaagca                                                                                                                                                >  2:181259/1‑116 (MQ=255)
                              ttGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGagcaagca                                                                                                                                                <  1:181259/116‑1 (MQ=255)
                                    aaTTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaattacttaatt                                                                                                           <  4:234983/147‑1 (MQ=255)
                                     aTTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCttaattacttaattac                                                                                                         <  2:156868/148‑1 (MQ=255)
                                      ttACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTaca                                                                                                        <  6:343475/148‑1 (MQ=255)
                                        aCAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGTACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTacac                                                                                                       >  6:379243/1‑147 (MQ=255)
                                           gCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTAcacca                                                                                                     >  4:184447/1‑146 (MQ=255)
                                           gCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTAcacca                                                                                                     >  7:318197/1‑146 (MQ=255)
                                                tttGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCAt                                                                                              <  2:178432/148‑1 (MQ=255)
                                                 ttGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCa                                                                                               >  6:197945/1‑146 (MQ=255)
                                                      aCAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGTAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaatt                                                                                                                   >  6:351071/1‑121 (MQ=255)
                                                      aCAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGTAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTcttaatt                                                                                                                   <  5:351071/121‑1 (MQ=255)
                                                       cAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTc                                                                                       <  8:318197/148‑1 (MQ=255)
                                                         gCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTAAATTACACCACAGGCATAGCGTTcac                                                                                     <  8:192194/148‑1 (MQ=255)
                                                                    catcatcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACATGCATAGCGTTCaccgccaccgcc                                                                           >  2:103898/1‑147 (MQ=255)
                                                                    catcatcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCaccgccaccgcc                                                                           >  7:296249/1‑147 (MQ=255)
                                                                    catcatcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCaccgccaccgcc                                                                           >  2:47119/1‑147 (MQ=255)
                                                                       catcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTCCGCTGGCCTCTTAATTACTTAATTACACCACAGGCCTAGCGTTCACCGCCACCGCCCAGc                                                                       >  3:282910/1‑148 (MQ=255)
                                                                       catcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGc                                                                       >  4:174566/1‑148 (MQ=255)
                                                                       catcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGc                                                                       >  3:74174/1‑148 (MQ=255)
                                                                        atcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCCCCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACAACAGGCATAGCGTTCACCGCCACCGCCCAGCg                                                                      >  5:199045/1‑148 (MQ=255)
                                                                        atcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGc                                                                       <  8:325712/147‑1 (MQ=255)
                                                                         tcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGc                                                                                                                                <  8:328313/89‑1 (MQ=255)
                                                                         tcaACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGc                                                                                                                                >  7:328313/1‑89 (MQ=255)
                                                                           aaCCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGtt                                                                   >  7:382426/1‑148 (MQ=255)
                                                                            aCCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGttt                                                                  <  7:76128/148‑1 (MQ=255)
                                                                            aCCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATCACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGttt                                                                  <  6:66648/148‑1 (MQ=255)
                                                                              cTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAg                                                                >  6:103387/1‑148 (MQ=255)
                                                                               tGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAgg                                                               >  8:303725/1‑148 (MQ=255)
                                                                                       gTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCgg                                                       <  7:303725/148‑1 (MQ=255)
                                                                                        tGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGa                                                      <  1:317074/148‑1 (MQ=255)
                                                                                           gCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCATCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCCGCGGTTTAGGTTGATCGGACa                                                    >  8:198092/1‑147 (MQ=255)
                                                                                               gcggcAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACAtatta                                               <  7:233048/148‑1 (MQ=255)
                                                                                               gcggcAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACAtatta                                               <  1:47119/148‑1 (MQ=255)
                                                                                               gcggcAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACAtatt                                                <  5:103387/147‑1 (MQ=255)
                                                                                               gcggcAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACAtatt                                                <  1:103898/147‑1 (MQ=255)
                                                                                                     aTTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATcgccg                                         >  5:305626/1‑148 (MQ=255)
                                                                                                          cTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAAc                                    >  6:124103/1‑148 (MQ=255)
                                                                                                            ggcggcTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGt                                  >  8:279970/1‑148 (MQ=255)
                                                                                                             gcggcTGTACTGCTCGCTTTCGCAACTTGAGCAAGTACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTg                                 <  5:379243/148‑1 (MQ=255)
                                                                                                                          tCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCt                    >  6:212804/1‑148 (MQ=255)
                                                                                                                            gCTTTCGCAACTTGAGCAAGTACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCttt                  >  7:218546/1‑148 (MQ=255)
                                                                                                                               ttCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGttt                                                                  >  3:161433/1‑97 (MQ=255)
                                                                                                                               ttCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGttt                                                                  <  4:161433/97‑1 (MQ=255)
                                                                                                                                  gCAACTTGAGCAAGTACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCtt             <  8:218546/147‑1 (MQ=255)
                                                                                                                                     aCTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGAt         <  3:184447/148‑1 (MQ=255)
                                                                                                                                     aCTTGAGCAAGCACCACCGCAAGGTACGATGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGAt         <  5:366280/148‑1 (MQ=255)
                                                                                                                                          agaaagcaCCACCGCAAGGTACGCTGGCCTCTTTATTACTTAATTACACCTCAGGCATAGCGTTCACCGCCACCTCCCTGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCa     <  7:198092/144‑1 (MQ=255)
                                                                                                                                            caagcaCCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATcc  <  8:382426/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACAGCGGATTTGCGACAGCCTTGCCCGTCATCATCAACCTGCCGCCCGTGAGCTGGCGGCAATTGTCTGGCGGCTGTACTGCTCGCTTTCGCAACTTGAGCAAGCACCACCGCAAGGTACGCTGGCCTCTTAATTACTTAATTACACCACAGGCATAGCGTTCACCGCCACCGCCCAGCGGTTTAGGTTGATCGGACATATTATCGCCGCCAACGTGGACCATCAGCGCTTTGTCTTTGATTTCATCC  >  NC_000913/1723961‑1724248

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