Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,759 G→A 100% intergenic (+80/‑11) 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,7590GA99.1% 358.4 / ‑4.4 108intergenic (+80/‑11)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 (50/57);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 108 total aligned reads displayed.

GTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810662‑2810854
                                                                                                 |                                                                                               
gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATc                                                                                               <  1:3144715/100‑1 (MQ=255)
gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATc                                                                                               >  2:348047/1‑100 (MQ=255)
  gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCAt                                                                                             >  1:2878479/1‑100 (MQ=255)
    tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAc                                                                                           <  1:2831073/100‑1 (MQ=255)
    tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAc                                                                                           <  1:2893963/100‑1 (MQ=255)
    tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAc                                                                                           <  1:113165/100‑1 (MQ=255)
    tGGCGATTATATAACCCACAAGAATCATTTCTCTAAAACCATACCTATACTTTACTTGTCAATGTAGTTACCAAATCGGCTGCAATTAATGAGATCATAc                                                                                           <  2:1442239/100‑1 (MQ=255)
     ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAcc                                                                                          <  2:1940698/100‑1 (MQ=255)
     ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAcc                                                                                          <  1:3337660/100‑1 (MQ=255)
      gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAccc                                                                                         <  2:2911147/100‑1 (MQ=255)
      gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAccc                                                                                         >  2:3400376/1‑100 (MQ=255)
       cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCa                                                                                        >  1:745812/1‑100 (MQ=255)
       cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCa                                                                                        >  1:1037992/1‑100 (MQ=255)
             tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATg                                                                                    <  1:1888085/98‑1 (MQ=255)
             tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCAAAATGGa                                                                                  >  2:1907257/1‑100 (MQ=255)
              ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGAt                                                                                 <  2:2300830/100‑1 (MQ=255)
              ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGAt                                                                                 >  1:2355416/1‑100 (MQ=255)
              ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGAt                                                                                 <  2:1023036/100‑1 (MQ=255)
               taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGACATTAATGAAATCATACCCAAATGGATa                                                                                >  2:3622223/1‑100 (MQ=255)
               taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATTAGATCATACCCAAATGGAt                                                                                 >  2:724409/1‑99 (MQ=255)
               taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATa                                                                                <  2:363708/100‑1 (MQ=255)
               taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATa                                                                                <  1:3466007/100‑1 (MQ=255)
               taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATa                                                                                <  2:245398/100‑1 (MQ=255)
                aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAg                                                                               <  1:1880883/100‑1 (MQ=255)
                 aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGt                                                                              >  1:985466/1‑100 (MQ=255)
                 aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTCTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGt                                                                              <  2:1986598/100‑1 (MQ=255)
                   ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTc                                                                            >  1:1538113/1‑100 (MQ=255)
                    cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                           <  1:1161201/100‑1 (MQ=255)
                    cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                           <  2:1452976/100‑1 (MQ=255)
                     acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTACGt                                                                          >  1:1247628/1‑100 (MQ=255)
                       aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGttt                                                                        <  2:2183398/100‑1 (MQ=255)
                       aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGttt                                                                        <  1:1687024/100‑1 (MQ=255)
                         gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTAc                                                                      >  1:1708254/1‑100 (MQ=255)
                         gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTAc                                                                      <  2:259676/100‑1 (MQ=255)
                         gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTAc                                                                      >  1:563884/1‑100 (MQ=255)
                         gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTAc                                                                      <  1:2813995/100‑1 (MQ=255)
                          aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACg                                                                     <  1:155134/100‑1 (MQ=255)
                           aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGc                                                                    <  1:1706652/100‑1 (MQ=255)
                           aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGc                                                                    <  1:2817878/100‑1 (MQ=255)
                            tCATTTTTCTAAAACAATACATTTACTTTATGTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGcc                                                                   <  1:2891466/100‑1 (MQ=255)
                              aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCa                                                                 >  2:383602/1‑100 (MQ=255)
                              aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCa                                                                 <  2:531855/100‑1 (MQ=255)
                               tctttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCACAt                                                                <  2:2203454/98‑1 (MQ=255)
                                ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                                <  1:2958706/99‑1 (MQ=255)
                                  ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                             <  1:3035869/100‑1 (MQ=255)
                                  ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                             <  1:930491/100‑1 (MQ=255)
                                  ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                             <  1:1719407/100‑1 (MQ=255)
                                   tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGTGATCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                               >  2:619217/1‑97 (MQ=255)
                                    cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           >  2:2178285/1‑100 (MQ=255)
                                    cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           >  1:2185532/1‑100 (MQ=255)
                                     tAATACAATACATGTACTTTATATATCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAAAGGATACTTCGTTTACGCCCATTGAAc                                                           <  2:2001672/99‑1 (MQ=255)
                                          ccaTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                     >  1:2465751/3‑100 (MQ=255)
                                           aaTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCt                                                    >  2:1186395/1‑100 (MQ=255)
                                            aTACATTTACTTTATTTGTCACCGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTa                                                   >  1:1736821/1‑100 (MQ=255)
                                               cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                <  2:470795/100‑1 (MQ=255)
                                               cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:3044206/1‑100 (MQ=255)
                                               cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:426251/1‑100 (MQ=255)
                                                aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                <  1:524771/99‑1 (MQ=255)
                                                aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                               >  1:2821700/1‑100 (MQ=255)
                                                  ttACTTTATTGGTCACTGGGATTACTATGTCTGCTGAAATTAATGAGATCATACCGAAATGGATAGAGCGTTGCGGCCCATTGAAGAAATGGTaaga                                                >  2:2220731/1‑95 (MQ=255)
                                                   tACTTTAATTGTCACTGTCGTTAGTATATCGGCTGACATCAATGAGATCAGACCCAACTGGATAGTTCGTTTACGCCCATTGAACAAATGCTgaaa                                                <  2:909322/96‑1 (MQ=255)
                                                    aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                           <  2:2613516/100‑1 (MQ=255)
                                                     cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          >  2:1501793/1‑100 (MQ=255)
                                                        tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                       <  2:1828227/100‑1 (MQ=255)
                                                        tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                       >  1:2445271/1‑100 (MQ=255)
                                                          tttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                     <  1:271080/100‑1 (MQ=255)
                                                              tCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgc                                 <  1:2932794/100‑1 (MQ=255)
                                                               cccTGTCGTTACTATATCGGCTGAAATGAATGGTATTATACCCTAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                <  2:3550460/98‑1 (MQ=255)
                                                               cACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                >  1:760760/1‑100 (MQ=255)
                                                               cACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                <  2:515903/100‑1 (MQ=255)
                                                               cACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCACTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                >  2:2638149/1‑100 (MQ=255)
                                                               cACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTAGTTTACGCCCATTGAACAAATGATAAAATTTCGCGCCAgccgcc                                >  2:2752794/1‑100 (MQ=255)
                                                                aCTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                               >  1:3047910/1‑100 (MQ=255)
                                                                 cTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                              >  1:3257811/1‑100 (MQ=255)
                                                                 cTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                              <  2:2079338/100‑1 (MQ=255)
                                                                 cTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                              >  2:229832/1‑100 (MQ=255)
                                                                     cGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                          >  1:254175/1‑100 (MQ=255)
                                                                      gTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                         <  2:942971/100‑1 (MQ=255)
                                                                        tACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                       >  1:315970/1‑100 (MQ=255)
                                                                         aCTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                      >  2:1047657/1‑100 (MQ=255)
                                                                          cTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTCCGTTTACGCCCAGCGAACAAATGCAAAAATTTCGCGCAAGCCGCCACGAAGa                         >  2:2056311/1‑96 (MQ=255)
                                                                           tctatCTGTTGAAATTCATGAAATCATACCCAAATGAATAGTTCGTTTAAGCCCATTGAACAAATGCTAAAATTTCGCTCCAGCCGCCACTAATAttt                      <  2:2814520/96‑1 (MQ=255)
                                                                           tataTCGGCTGAAATTAATGAGATCATACCCTAGTGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAGATTTCGCGCCAGCCGTCACGAAGATTTTc                    <  1:744995/100‑1 (MQ=255)
                                                                           tataTCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                    >  2:3551762/1‑100 (MQ=255)
                                                                            atatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcn                   <  1:1374344/100‑2 (MQ=255)
                                                                               tCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa                >  2:2700738/1‑100 (MQ=255)
                                                                               tCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa                <  1:194431/100‑1 (MQ=255)
                                                                                 ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              >  1:959915/1‑100 (MQ=255)
                                                                                 ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              >  1:1675322/1‑100 (MQ=255)
                                                                                   cTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg            >  1:418987/1‑100 (MQ=255)
                                                                                    tGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg           >  1:3318341/1‑100 (MQ=255)
                                                                                       aaTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga        >  1:3290016/1‑100 (MQ=255)
                                                                                        aTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       >  1:1938654/1‑100 (MQ=255)
                                                                                        aTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       >  1:1938235/1‑100 (MQ=255)
                                                                                          tAATTAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTACCAGGAGATcc     <  2:165341/100‑1 (MQ=255)
                                                                                          tAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATcc     >  1:415897/1‑100 (MQ=255)
                                                                                            aTGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt   >  1:2429447/1‑100 (MQ=255)
                                                                                            aTGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCtt   >  1:2095064/1‑100 (MQ=255)
                                                                                             tGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc  <  2:3179668/100‑1 (MQ=255)
                                                                                             tGAGATCATACCCAAATGGATAGTTCGTTGACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCTGCCACGAAGATTTTCCTTATCAGGAGATCCTTc  <  2:2185518/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810662‑2810854

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: