Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,760 T→C 100% intergenic (+81/‑10) ygaH → / → mprA putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,7600TC99.4% 499.3 / ‑5.5 145intergenic (+81/‑10)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (84/60);  total (85/60)
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.12e-01

Only 100 of 145 total aligned reads displayed.

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTG  >  NC_000913/2810663‑2810856
                                                                                                 |                                                                                                
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCa                                                                                                <  2:1165731/100‑1 (MQ=255)
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCa                                                                                                <  2:1795791/100‑1 (MQ=255)
 gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCAt                                                                                               >  2:2732972/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                             <  1:1160597/100‑1 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                             <  1:1383623/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACTATACATTTACTTTATTTGTCACTGTCGTTACTATATCTTCTGAAATTAATGAGGCCATa                                                                                              >  1:728915/1‑98 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAcc                                                                                            <  2:769149/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAcc                                                                                            <  2:269247/100‑1 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAccc                                                                                           >  2:1702530/1‑100 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACGTTGTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAccc                                                                                           <  2:2105589/100‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCa                                                                                          >  1:441670/1‑100 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                      >  2:1006309/1‑100 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                      <  2:2644713/100‑1 (MQ=255)
            tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGa                                                                                    <  2:1004482/100‑1 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGTAt                                                                                   >  1:2342567/1‑100 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGAt                                                                                   <  1:1719154/100‑1 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGACATTAATGAGGCCATACCCAAATGGATa                                                                                  >  1:3573371/1‑100 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATa                                                                                  >  2:2352835/1‑100 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATa                                                                                  >  2:1305509/1‑100 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGt                                                                                >  1:3181490/1‑100 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGt                                                                                >  2:1607900/1‑100 (MQ=255)
                 cctacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGACCATACACAAATGGATAg                                                                                 >  1:148184/4‑98 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGtt                                                                               >  2:3576264/1‑100 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGtt                                                                               >  1:3569981/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             >  2:1753742/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             >  2:2239329/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             >  2:2196745/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             <  1:3549033/100‑1 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             >  1:3276212/1‑100 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGACGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGtt                                                                           >  1:2075184/1‑100 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTa                                                                         >  1:3390489/1‑100 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTa                                                                         <  1:160319/100‑1 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGCGGCCATACCCAAATGGATAGTTCGTTTAc                                                                        >  1:1694952/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACg                                                                       >  1:2028090/1‑100 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGcc                                                                     <  2:3215792/100‑1 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGcc                                                                     <  1:1838944/100‑1 (MQ=255)
                           tCATTTTTCTAAAACAATACAGTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                          <  1:645504/95‑1 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGccc                                                                    >  1:2782002/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGccc                                                                    >  2:1606739/1‑100 (MQ=255)
                             aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCTTTACTATATCGGCTGAAATTAATGAGGCAATACCAAAATGGATAGTTCGTTTACGcc                                                                     >  1:1915840/1‑98 (MQ=255)
                              tttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                                  >  2:1490372/1‑100 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                                <  1:627241/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                                <  1:1283743/100‑1 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                               >  1:2842625/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGaa                                                              >  2:3491390/1‑99 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             <  1:3656595/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             >  1:76205/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             <  1:228416/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             <  2:3493033/100‑1 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                           >  1:1934199/1‑100 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                           <  2:2875291/100‑1 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                          >  2:3412025/1‑100 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                        <  2:3594019/100‑1 (MQ=255)
                                         aaaTACATTTACTTTATTTGTCACGGTCGTTACTATATAGGGTGAATTTAATGAGGCCATAACCACATGGATAGTTCGTTGACGCGCATTGAAc                                                             >  1:1737476/2‑94 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCCTACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                  >  1:1230590/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                  >  1:458803/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                  >  1:2563490/1‑100 (MQ=255)
                                               aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                                 <  2:2342580/100‑1 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                             <  2:1694959/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                            <  1:449398/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                             >  2:1615958/1‑99 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                           <  2:559603/100‑1 (MQ=255)
                                                      ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                          <  2:849492/100‑1 (MQ=255)
                                                       tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                         >  2:1625938/1‑100 (MQ=255)
                                                        aTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCa                                        >  1:2798384/1‑100 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                     <  1:2753860/100‑1 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                     >  1:2096389/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                  <  2:894815/100‑1 (MQ=255)
                                                                cTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                >  1:3171071/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                              >  1:3309574/1‑100 (MQ=255)
                                                                  gTCATTAATATATCGGCTGAAATTAATGAGCCCACACCCCAATGGATAGTTCGTTTACGCCCCTTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACg                               <  2:2980003/99‑1 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                            >  2:2582135/1‑100 (MQ=255)
                                                                     cttACTACATCGGCTGAATTTCATGAGGCCATACCTAAATGGATAGTTCGTTTACGCCCAGTTAACAAATGCTAAAATTTCTCGCCAGCCGCCACGAAg                            >  1:1617035/2‑99 (MQ=255)
                                                                       tACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                         <  2:148646/100‑1 (MQ=255)
                                                                       tACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                         <  1:114599/100‑1 (MQ=255)
                                                                         cTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt                       >  2:1416593/1‑100 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  1:3443840/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     >  1:2948296/1‑100 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  2:441672/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  2:1816149/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                    >  1:1980513/1‑100 (MQ=255)
                                                                             atCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                   >  1:414454/1‑100 (MQ=255)
                                                                             atCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                   >  1:2629689/1‑100 (MQ=255)
                                                                             atCGACTGAAATTAATGAGGCCATACCCAAATGGCTACTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCCCGAATATCTTCCtt                   >  1:371165/1‑100 (MQ=255)
                                                                              tCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAAGAAGATTTTACTTa                  >  1:2855109/1‑100 (MQ=255)
                                                                               cGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt                 >  2:2565592/1‑100 (MQ=255)
                                                                                ggCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                >  1:2232165/1‑100 (MQ=255)
                                                                                ggCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                >  2:3097291/1‑100 (MQ=255)
                                                                                  cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg              <  2:1258652/100‑1 (MQ=255)
                                                                                      aaTTAATGAGGCTATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGTCAGCCGCCACGAAGATTTTCCTTAGCAGgaga          >  1:112703/1‑100 (MQ=255)
                                                                                      aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga          >  1:2662951/1‑100 (MQ=255)
                                                                                      aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga          <  1:2582120/100‑1 (MQ=255)
                                                                                      aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag           >  1:895107/1‑99 (MQ=255)
                                                                                       aTTAATGAGGCAATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt         <  1:3012225/100‑1 (MQ=255)
                                                                                        ttAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCATTATCAGGAGAt         >  1:1805068/1‑99 (MQ=255)
                                                                                            tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc    >  2:439014/1‑100 (MQ=255)
                                                                                            tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc    >  2:1185604/1‑100 (MQ=255)
                                                                                            tGAGGCCAGACCCAAATGGATAGTTCTTTTACGCCCATTGAACAAATGCTAAAAATTCGCGCCAGCCGCCACTAAAATTTTCCCTAACAGGAGATcactc    >  1:871446/1‑96 (MQ=255)
                                                                                              aGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg  >  1:531622/1‑100 (MQ=255)
                                                                                              aGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg  <  2:2051510/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTG  >  NC_000913/2810663‑2810856

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