Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,756 G→A 100% intergenic (+77/‑14) 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,7560GA99.4% 534.2 / ‑4.9 162intergenic (+77/‑14)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 (0/1);  new base A (84/76);  total (85/77)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 162 total aligned reads displayed.

AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810659‑2810852
                                                                                                 |                                                                                                
aaaGTGATGGCGATTATATAACCCACACGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTCCTATATCGGCTGAAATTAATAAg                                                                                                <  2:3922035/100‑1 (MQ=255)
 aaGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAgg                                                                                               <  1:1036309/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)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCAt                                                                                           >  2:2841378/1‑100 (MQ=255)
     gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCAt                                                                                           >  2:1827723/1‑100 (MQ=255)
       tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAc                                                                                         >  2:3392141/1‑100 (MQ=255)
        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAcc                                                                                        <  1:2853958/100‑1 (MQ=255)
        ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAc                                                                                         >  2:884785/1‑99 (MQ=255)
         gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAccc                                                                                       >  1:3254263/1‑100 (MQ=255)
         gCGATCATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAccc                                                                                       >  1:559819/1‑100 (MQ=255)
           gATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACCATATCGGCTGAAATTAATAAGGTCATACCCaa                                                                                     >  1:3056299/1‑100 (MQ=255)
             ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAAt                                                                                   <  2:1815481/100‑1 (MQ=255)
              tatataACCCACCAGAATCATTTTTCTTAAAGAATACGTTTAGTTTCTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATATAGTCATACCCAAATg                                                                                  <  2:46243/100‑1 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGCCATACCCAAATgg                                                                                 >  1:187165/1‑100 (MQ=255)
               atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTAGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATgg                                                                                 >  1:2404449/1‑100 (MQ=255)
                 ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGAt                                                                               >  2:2746775/1‑100 (MQ=255)
                  taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATa                                                                              >  1:1484016/1‑100 (MQ=255)
                  taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                              >  2:2069108/1‑100 (MQ=255)
                  taACCCACAAGAATAATTTTTCTAAAACAATAAATTAAATTTATTGGTTACGGTCGTTACTATATAGTATGATATTAATAAGGTAATAGGCAAATGAATa                                                                              >  2:557222/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)
                        acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGGAATTAATAAGGTCATACCCAAATGGATAGTTCGt                                                                        <  1:2528177/100‑1 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGtt                                                                       <  1:158423/100‑1 (MQ=255)
                         caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGtt                                                                       >  1:1774817/1‑100 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                    >  1:64517/1‑100 (MQ=255)
                            gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTAc                                                                    <  1:1984897/100‑1 (MQ=255)
                              aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACg                                                                   <  2:3963643/99‑1 (MQ=255)
                              aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                                  <  1:534295/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)
                                  ttttttTAAAATAATTCATTTAATTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCTTTTATGctaat                                                              >  2:536911/1‑96 (MQ=255)
                                   ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                             <  2:1650256/100‑1 (MQ=255)
                                    tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                            <  1:3579766/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                                                       >  2:3285113/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                      >  1:2885604/1‑100 (MQ=255)
                                          aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                      >  2:2939418/1‑100 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                    <  2:2282062/100‑1 (MQ=255)
                                            aCAATACATTTACTTTATTTGTCACTGGCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCATATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                    <  1:1552657/100‑1 (MQ=255)
                                              aaTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCt                                                  >  2:197950/1‑100 (MQ=255)
                                                  cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                              <  2:3104834/100‑1 (MQ=255)
                                                  cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGCACAAATGATaaaa                                              >  2:35030/1‑100 (MQ=255)
                                                   aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             >  2:26491/1‑100 (MQ=255)
                                                   aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             >  2:112946/1‑100 (MQ=255)
                                                   aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             >  1:3431580/1‑100 (MQ=255)
                                                     ttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAttt                                           <  2:1634870/100‑1 (MQ=255)
                                                      tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                          >  1:3016961/1‑100 (MQ=255)
                                                      tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                          >  1:3876909/1‑100 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  1:1975762/100‑1 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  2:3248862/100‑1 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  1:2841394/100‑1 (MQ=255)
                                                       aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         >  2:3693799/1‑100 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGGAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  2:3334438/1‑100 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTGATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTACAATTTcgc                                        >  2:2475992/1‑100 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        <  1:893179/100‑1 (MQ=255)
                                                        cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  2:3484718/1‑100 (MQ=255)
                                                           tATTTGTCACTGTCGTTACTATATTGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                     <  2:267696/100‑1 (MQ=255)
                                                             tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCTCCAg                                   >  2:951127/1‑100 (MQ=255)
                                                              ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                                  <  2:1628200/100‑1 (MQ=255)
                                                               tGTCACTGTCGTTACTATATCGGCTGCAATTAATACGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                 >  2:2687513/1‑100 (MQ=255)
                                                               tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                 >  1:3773848/1‑100 (MQ=255)
                                                               tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                 >  1:2238989/1‑100 (MQ=255)
                                                                  caaggTCGTTACTATATCTGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGGTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                              <  2:3883231/96‑1 (MQ=255)
                                                                   aCTGTCGTTACTATATCGGCTGAAATTAATACGCTCATACCCAAATGGATAGTTCGTTTACGCTCATTGAACAAATGCTAACATTTCGCGCCGGCCCCCa                             >  2:205929/1‑100 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                            <  1:2161340/100‑1 (MQ=255)
                                                                    cTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGCTAGTTCGTGTACGCCCATCGAACAAAGGATAAAATTTCGCGCCAgccgc                               >  2:3087919/1‑97 (MQ=255)
                                                                     tGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                            >  2:2829694/1‑99 (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:2551364/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        <  2:3773824/100‑1 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        >  2:3636898/1‑100 (MQ=255)
                                                                        cGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        >  1:3831053/1‑100 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                       <  1:1842083/100‑1 (MQ=255)
                                                                         gTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                       >  2:1452452/1‑100 (MQ=255)
                                                                          ttACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAt                      <  1:3472783/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:474360/100‑1 (MQ=255)
                                                                               atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:78700/100‑1 (MQ=255)
                                                                                tatCGGCTGAAATTAATAAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                <  2:1776932/100‑1 (MQ=255)
                                                                                  tCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa              <  2:2312770/100‑1 (MQ=255)
                                                                                   cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt             >  2:903722/1‑100 (MQ=255)
                                                                                     gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa           <  2:1773301/100‑1 (MQ=255)
                                                                                      cTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg          <  1:3241578/100‑1 (MQ=255)
                                                                                       tGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg          <  1:2963023/99‑1 (MQ=255)
                                                                                        gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga        >  1:1238832/1‑100 (MQ=255)
                                                                                          aaTTAATAAGGCAATACCCAAATGAATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga      <  2:3334132/100‑1 (MQ=255)
                                                                                           aTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     >  1:1742443/1‑100 (MQ=255)
                                                                                           aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     <  1:1997008/100‑1 (MQ=255)
                                                                                           aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     <  1:332554/100‑1 (MQ=255)
                                                                                              aaTAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt  <  2:1254878/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AAAGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810659‑2810852

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: