Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,758 G→A 18.4% intergenic (+79/‑12) 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,7580GA18.4% 358.6 / 65.1 167intergenic (+79/‑12)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (73/63);  new base A (16/15);  total (89/78)
Fisher's exact test for biased strand distribution p-value = 8.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 172 total aligned reads displayed.

AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGA  >  NC_000913/2810659‑2810857
                                                                                                   |                                                                                                    
aaaGTGATGGCGATTATATAACCCACACGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTCCTATATCGGCTGAAATTAATAAg                                                                                                      <  2:3922035/100‑1 (MQ=255)
 aaGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAgg                                                                                                     >  2:594593/1‑100 (MQ=255)
 aaGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAgg                                                                                                     <  1:1036309/100‑1 (MQ=255)
 aaGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAgg                                                                                                     <  1:3821764/100‑1 (MQ=255)
  aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                                    <  2:1005467/100‑1 (MQ=255)
  aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                                    <  1:3570747/100‑1 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGT‑c                                                                                                  <  2:1393335/100‑1 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑c                                                                                                  >  1:1047279/1‑100 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑c                                                                                                  >  1:1734242/1‑100 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑c                                                                                                  <  1:3622203/100‑1 (MQ=255)
    tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑Ca                                                                                                 >  2:1579664/1‑100 (MQ=255)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CAt                                                                                                >  2:2841378/1‑100 (MQ=255)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CAt                                                                                                >  2:1827723/1‑100 (MQ=255)
      aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATa                                                                                               <  2:2022255/100‑1 (MQ=255)
        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATAcc                                                                                             <  1:2853958/100‑1 (MQ=255)
         gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATAccc                                                                                            >  2:1758751/1‑100 (MQ=255)
         gCGATCATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATAccc                                                                                            >  1:559819/1‑100 (MQ=255)
             ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAAt                                                                                        <  2:1815481/100‑1 (MQ=255)
              tatataACCCACCAGAATCATTTTTCTTAAAGAATACGTTTAGTTTCTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATATAGT‑CATACCCAAATg                                                                                       <  2:46243/100‑1 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATgg                                                                                      >  1:3251422/1‑100 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGC‑CATACCCAAATgg                                                                                      >  1:187165/1‑100 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTAGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATgg                                                                                      >  1:2404449/1‑100 (MQ=255)
                 ataACCCACAAGAATCATTTTTCTAATACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGAt                                                                                    <  2:2059660/100‑1 (MQ=255)
                 ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGAt                                                                                    >  2:2746775/1‑100 (MQ=255)
                  taACCCACAAGAATCATTTTTTTTTTTTTTTTCATTTACTTTATTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAGTT‑CATACCCAAATGtata                                                                                   >  2:1386372/1‑100 (MQ=255)
                   aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATATGGT‑CATACCCAAATGGATAg                                                                                  <  2:2521125/100‑1 (MQ=255)
                    aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGt                                                                                 >  1:3369972/1‑100 (MQ=255)
                       cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCg                                                                              >  1:2733603/1‑100 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGtt                                                                            >  1:1774817/1‑100 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACAGTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTGCGtt                                                                            <  2:1914563/100‑1 (MQ=255)
                           aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTAATATAGCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGttt                                                                           <  2:3452605/99‑1 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTAc                                                                         <  1:1984897/100‑1 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTAc                                                                         >  1:64517/1‑100 (MQ=255)
                             aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACCATCGTTACTATATCGGCTGAAATTAATAAGGT‑CGTACCCAAATGGATAGTTCGTTTACg                                                                        <  2:3132199/100‑1 (MQ=255)
                              aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACg                                                                        <  2:3963643/99‑1 (MQ=255)
                               tCATTTTTCTAAGACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:3380078/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:255381/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:292246/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGcc                                                                      >  1:1564659/1‑100 (MQ=255)
                                cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGccc                                                                     >  1:3180289/1‑100 (MQ=255)
                                 aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCa                                                                    >  1:497473/1‑100 (MQ=255)
                                  ttttttTAAAATAATTCATTTAATTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCTTTTATGctaat                                                                   >  2:536911/1‑96 (MQ=255)
                                   ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                                  <  2:3570614/100‑1 (MQ=255)
                                     ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                                <  1:919562/100‑1 (MQ=255)
                                       cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTCAAc                                                              >  2:1104783/1‑100 (MQ=255)
                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCTGCTGAAATTATTAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                             >  2:3219767/1‑100 (MQ=255)
                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                             >  1:497407/1‑100 (MQ=255)
                                         aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                            >  1:3393417/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                           >  1:2885604/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAAT‑AGGTCCATACCCAACTGGATAGTTCGTTTACGCCCATTGAACaaa                                                           >  2:2668907/1‑100 (MQ=255)
                                           aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                          <  1:2414034/100‑1 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                         <  2:2282062/100‑1 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGGCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCATATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                         <  1:1552657/100‑1 (MQ=255)
                                                tACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaa                                                     <  1:3075332/100‑1 (MQ=255)
                                                  cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGCACAAATGATaaaa                                                   >  2:35030/1‑100 (MQ=255)
                                                     ttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAttt                                                <  2:1634870/100‑1 (MQ=255)
                                                      tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                               >  1:3016961/1‑100 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                              <  1:1975762/100‑1 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                              >  2:3693799/1‑100 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGGAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                             >  2:3334438/1‑100 (MQ=255)
                                                         tttATTTGTCACTGTCGTTACTAGATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCACATTGAACAAATGCTAAAATTTcgcg                                            <  2:3289675/100‑1 (MQ=255)
                                                           tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTGGTTTAAGCCCATTGAACAAATGCTAAAATTTCGCGcc                                          >  2:2793522/1‑100 (MQ=255)
                                                             tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCTCCAg                                        >  2:951127/1‑100 (MQ=255)
                                                               tGTCACTGTCGTTACTATATCGGCTGCAATTAATACGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                      >  2:2687513/1‑100 (MQ=255)
                                                                  caaggTCGTTACTATATCTGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGGTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                   <  2:3883231/96‑1 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                 <  1:2161340/100‑1 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGCTAGTTCGTGTACGCCCATCGAACAAAGGATAAAATTTCGCGCCAgccgc                                    >  2:3087919/1‑97 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                 <  2:121574/100‑1 (MQ=255)
                                                                      gTCGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                               >  1:2970383/1‑100 (MQ=255)
                                                                       tCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                              >  2:3024366/1‑100 (MQ=255)
                                                                       tCGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTGCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                              <  2:2039344/100‑1 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3444389/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3570887/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             <  2:3773824/100‑1 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:2551364/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3831053/1‑100 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                            <  1:1842083/100‑1 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                            >  2:237319/1‑100 (MQ=255)
                                                                            aCTATATCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                         <  2:3282265/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  2:78700/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  2:474360/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  1:1542059/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      >  2:2989306/1‑100 (MQ=255)
                                                                               atatCGCCTCAAATTAATAAGGT‑CATACCCAAATTGATAGTTCTTTTACGCCCATTGAACAAATGCTAAAATTTTGCGCCAGCCGCCACGAAGATTTTcc                      <  2:86467/100‑1 (MQ=255)
                                                                                tatCGGCTGAAATTAATAAGGC‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                     <  2:1776932/100‑1 (MQ=255)
                                                                                   cGGCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt                  >  2:903722/1‑100 (MQ=255)
                                                                                     gCTGAAATTAATAAAGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa                <  2:1643306/100‑1 (MQ=255)
                                                                                        gAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             >  1:2075251/1‑100 (MQ=255)
                                                                                        gAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             <  1:2964891/100‑1 (MQ=255)
                                                                                        gAAATTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             >  1:1238832/1‑100 (MQ=255)
                                                                                          aaTTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga           <  1:330871/100‑1 (MQ=255)
                                                                                           aTTAATAAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt          >  1:1742443/1‑100 (MQ=255)
                                                                                               ataaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:3673290/4‑100 (MQ=255)
                                                                                               ataaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:2817254/4‑100 (MQ=255)
                                                                                               ataaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      <  2:217654/97‑1 (MQ=255)
                                                                                               ataaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:1601093/4‑100 (MQ=255)
                                                                                                taaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc     <  1:4012371/98‑1 (MQ=255)
                                                                                                taaGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc     >  2:2275216/3‑100 (MQ=255)
                                                                                                  aGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg   <  1:2331750/100‑1 (MQ=255)
                                                                                                   ggT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGa  <  1:2829710/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGT‑CATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGA  >  NC_000913/2810659‑2810857

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: ‑