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.

AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGA  >  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)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTc                                                                                                  <  2:1393335/100‑1 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                                  >  1:1047279/1‑100 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTc                                                                                                  >  1:1734242/1‑100 (MQ=255)
   gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTc                                                                                                  <  1:3622203/100‑1 (MQ=255)
    tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCa                                                                                                 >  2:1579664/1‑100 (MQ=255)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCAt                                                                                                >  2:2841378/1‑100 (MQ=255)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCAt                                                                                                >  2:1827723/1‑100 (MQ=255)
      aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATa                                                                                               <  2:2022255/100‑1 (MQ=255)
        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAcc                                                                                             <  1:2853958/100‑1 (MQ=255)
         gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAccc                                                                                            >  2:1758751/1‑100 (MQ=255)
         gCGATCATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAccc                                                                                            >  1:559819/1‑100 (MQ=255)
             ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAAt                                                                                        <  2:1815481/100‑1 (MQ=255)
              tatataACCCACCAGAATCATTTTTCTTAAAGAATACGTTTAGTTTCTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATATAGTCATACCCAAATg                                                                                       <  2:46243/100‑1 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATgg                                                                                      >  1:3251422/1‑100 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGCCATACCCAAATgg                                                                                      >  1:187165/1‑100 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTAGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATgg                                                                                      >  1:2404449/1‑100 (MQ=255)
                 ataACCCACAAGAATCATTTTTCTAATACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGAt                                                                                    <  2:2059660/100‑1 (MQ=255)
                 ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGAt                                                                                    >  2:2746775/1‑100 (MQ=255)
                  taACCCACAAGAATCATTTTTTTTTTTTTTTTCATTTACTTTATTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAGTTCATACCCAAATGtata                                                                                   >  2:1386372/1‑100 (MQ=255)
                   aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATATGGTCATACCCAAATGGATAg                                                                                  <  2:2521125/100‑1 (MQ=255)
                    aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGt                                                                                 >  1:3369972/1‑100 (MQ=255)
                       cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCg                                                                              >  1:2733603/1‑100 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGtt                                                                            >  1:1774817/1‑100 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACAGTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTGCGtt                                                                            <  2:1914563/100‑1 (MQ=255)
                           aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTAATATAGCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGttt                                                                           <  2:3452605/99‑1 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                         <  1:1984897/100‑1 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                         >  1:64517/1‑100 (MQ=255)
                             aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACCATCGTTACTATATCGGCTGAAATTAATAAGGTCGTACCCAAATGGATAGTTCGTTTACg                                                                        <  2:3132199/100‑1 (MQ=255)
                              aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACg                                                                        <  2:3963643/99‑1 (MQ=255)
                               tCATTTTTCTAAGACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:3380078/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:255381/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                      <  2:292246/100‑1 (MQ=255)
                               tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                      >  1:1564659/1‑100 (MQ=255)
                                cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGccc                                                                     >  1:3180289/1‑100 (MQ=255)
                                 aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCa                                                                    >  1:497473/1‑100 (MQ=255)
                                  ttttttTAAAATAATTCATTTAATTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCTTTTATGctaat                                                                   >  2:536911/1‑96 (MQ=255)
                                   ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                                  <  2:3570614/100‑1 (MQ=255)
                                     ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                                <  1:919562/100‑1 (MQ=255)
                                       cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTCAAc                                                              >  2:1104783/1‑100 (MQ=255)
                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCTGCTGAAATTATTAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                             >  2:3219767/1‑100 (MQ=255)
                                        tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                             >  1:497407/1‑100 (MQ=255)
                                         aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                            >  1:3393417/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                           >  1:2885604/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAGGTCCATACCCAACTGGATAGTTCGTTTACGCCCATTGAACaaa                                                           >  2:2668907/1‑100 (MQ=255)
                                           aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                          <  1:2414034/100‑1 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                         <  2:2282062/100‑1 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGGCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCATATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                         <  1:1552657/100‑1 (MQ=255)
                                                tACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaa                                                     <  1:3075332/100‑1 (MQ=255)
                                                  cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGCACAAATGATaaaa                                                   >  2:35030/1‑100 (MQ=255)
                                                     ttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAttt                                                <  2:1634870/100‑1 (MQ=255)
                                                      tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                               >  1:3016961/1‑100 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                              <  1:1975762/100‑1 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                              >  2:3693799/1‑100 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGGAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                             >  2:3334438/1‑100 (MQ=255)
                                                         tttATTTGTCACTGTCGTTACTAGATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCACATTGAACAAATGCTAAAATTTcgcg                                            <  2:3289675/100‑1 (MQ=255)
                                                           tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTGGTTTAAGCCCATTGAACAAATGCTAAAATTTCGCGcc                                          >  2:2793522/1‑100 (MQ=255)
                                                             tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCTCCAg                                        >  2:951127/1‑100 (MQ=255)
                                                               tGTCACTGTCGTTACTATATCGGCTGCAATTAATACGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                      >  2:2687513/1‑100 (MQ=255)
                                                                  caaggTCGTTACTATATCTGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGGTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                   <  2:3883231/96‑1 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                 <  1:2161340/100‑1 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGCTAGTTCGTGTACGCCCATCGAACAAAGGATAAAATTTCGCGCCAgccgc                                    >  2:3087919/1‑97 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                 <  2:121574/100‑1 (MQ=255)
                                                                      gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                               >  1:2970383/1‑100 (MQ=255)
                                                                       tCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                              >  2:3024366/1‑100 (MQ=255)
                                                                       tCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTGCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                              <  2:2039344/100‑1 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3444389/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3570887/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             <  2:3773824/100‑1 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:2551364/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                             >  1:3831053/1‑100 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                            <  1:1842083/100‑1 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                            >  2:237319/1‑100 (MQ=255)
                                                                            aCTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                         <  2:3282265/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  2:78700/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  2:474360/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      <  1:1542059/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                      >  2:2989306/1‑100 (MQ=255)
                                                                               atatCGCCTCAAATTAATAAGGTCATACCCAAATTGATAGTTCTTTTACGCCCATTGAACAAATGCTAAAATTTTGCGCCAGCCGCCACGAAGATTTTcc                      <  2:86467/100‑1 (MQ=255)
                                                                                tatCGGCTGAAATTAATAAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                     <  2:1776932/100‑1 (MQ=255)
                                                                                   cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt                  >  2:903722/1‑100 (MQ=255)
                                                                                     gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa                <  2:1643306/100‑1 (MQ=255)
                                                                                        gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             >  1:2075251/1‑100 (MQ=255)
                                                                                        gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             <  1:2964891/100‑1 (MQ=255)
                                                                                        gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga             >  1:1238832/1‑100 (MQ=255)
                                                                                          aaTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga           <  1:330871/100‑1 (MQ=255)
                                                                                           aTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt          >  1:1742443/1‑100 (MQ=255)
                                                                                               ataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:3673290/4‑100 (MQ=255)
                                                                                               ataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:2817254/4‑100 (MQ=255)
                                                                                               ataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      <  2:217654/97‑1 (MQ=255)
                                                                                               ataaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt      >  1:1601093/4‑100 (MQ=255)
                                                                                                taaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc     <  1:4012371/98‑1 (MQ=255)
                                                                                                taaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc     >  2:2275216/3‑100 (MQ=255)
                                                                                                  aGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg   <  1:2331750/100‑1 (MQ=255)
                                                                                                   ggTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGa  <  1:2829710/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGA  >  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: