Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,872 C→T 17.4% H35Y (CAC→TAC) ‡ mprA → 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,8720CT17.4% 296.6 / 48.0 131H35Y (CAC→TAC) mprAtranscriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base C (50/58);  new base T (12/11);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 161 total aligned reads displayed.

GATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGG  >  NC_000913/2810773‑2810971
                                                                                                   |                                                                                                   
   aGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcaca                                                                                                  <  1:2312551/100‑1 (MQ=255)
   aGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcaca                                                                                                  <  2:142757/100‑1 (MQ=255)
    gTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacat                                                                                                 >  2:3224151/1‑100 (MQ=255)
     ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatg                                                                                                <  1:236760/100‑1 (MQ=255)
     ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatg                                                                                                <  1:224715/100‑1 (MQ=255)
      tCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgc                                                                                               <  2:1938644/100‑1 (MQ=255)
       cGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                              <  1:2659625/100‑1 (MQ=255)
       cGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                              >  1:2551081/1‑100 (MQ=255)
          ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaag                                                                                           <  1:1797757/100‑1 (MQ=255)
          ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaag                                                                                           <  2:550936/100‑1 (MQ=255)
           tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                          <  1:348050/100‑2 (MQ=255)
           tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                          <  2:1284164/100‑2 (MQ=255)
           tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                          <  2:254175/100‑2 (MQ=255)
           tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagc                                                                                          >  2:2754869/1‑100 (MQ=255)
            actcccATTAAACAAATGCTAAAATTTCGCGTCAGCCGCCACGAAGATTTTCCTTATCATTATATCCTTCTTACTCGTCTTTTCATGTGCATGCAaagac                                                                                         <  2:130940/97‑3 (MQ=255)
            aCGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                         <  1:3170817/100‑3 (MQ=255)
             ctcccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                         <  2:3312350/97‑3 (MQ=255)
             cGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                        <  1:43704/100‑1 (MQ=255)
               cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaagc                                                                                      >  2:3587009/1‑100 (MQ=255)
               cccATTGAACAAATGATAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                         <  2:1928807/97‑3 (MQ=255)
                ccATTGAACAAATGCTAAAATCTCGCGCCAGCCGCCACTAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagacgag                                                                                      <  1:2031607/99‑6 (MQ=38)
                    tGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGctgctg                                                                                 <  1:838762/100‑1 (MQ=255)
                     gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                <  2:1037984/100‑1 (MQ=255)
                     gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                >  2:564169/1‑100 (MQ=255)
                      aaCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                               <  2:3161594/100‑1 (MQ=255)
                      aaCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                               <  1:1890378/100‑1 (MQ=255)
                       aCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGgag                                                                              >  1:887295/1‑100 (MQ=255)
                        cAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGgaga                                                                             >  1:509877/1‑100 (MQ=255)
                        cAAATGCTAAAATCTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                               <  2:3326080/98‑1 (MQ=255)
                            tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                         <  1:2925460/100‑1 (MQ=255)
                            tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                         <  2:455739/100‑1 (MQ=255)
                            tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                         <  2:1766079/100‑1 (MQ=255)
                            tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                         <  2:454769/100‑1 (MQ=255)
                            tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                         >  2:3433533/1‑100 (MQ=255)
                             gatAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                        >  1:2454821/3‑100 (MQ=255)
                              cTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                       >  1:700135/1‑100 (MQ=255)
                              cTAAAATATCGTGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                       >  2:1910805/1‑100 (MQ=255)
                               tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                      <  2:1553175/100‑1 (MQ=255)
                               tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                      <  2:2768915/100‑1 (MQ=255)
                                aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAAt                                                                     <  1:1907267/100‑1 (MQ=255)
                                aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAAt                                                                     >  1:1227669/1‑100 (MQ=255)
                                       cgcgCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                              <  1:3419342/100‑1 (MQ=255)
                                       cgcgCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                              <  1:564171/100‑1 (MQ=255)
                                        gcgcCAGCCGCCACGACGATTTTCTTTATCAGGAGATCCTTCTGACTCGACTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                              <  2:843763/99‑1 (MQ=255)
                                        gcgcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCt                                                             >  2:2542705/1‑100 (MQ=255)
                                         cgcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTg                                                            >  2:2470108/1‑100 (MQ=255)
                                          gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGa                                                           <  1:1567469/100‑1 (MQ=255)
                                           ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                          <  2:2095051/100‑1 (MQ=255)
                                           ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                          >  1:1339134/1‑100 (MQ=255)
                                              gccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGc                                                       >  2:729282/1‑100 (MQ=255)
                                               ccgccACGAAGATTTTCCTTATCAGGAGCTCCTTCTGACTCGTCTTTGCATGCACATGAAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAgg                                                        >  2:516966/1‑98 (MQ=255)
                                               ccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAACGCAAGCTGCTGGAGACCCGCAATAAAATGCTAAAGGCt                                                      >  2:2132563/1‑100 (MQ=255)
                                                 gccACGAAGCTTTTCCTTATCAGGAGATCCTTCTGACTCGGCTTTGCATGCACATGCAACGCAAGCTGCTGGGGAACCGCAGTAAAATGCTGAAgg                                                        >  2:1713508/1‑96 (MQ=255)
                                                    aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAggg                                                 >  1:781514/1‑100 (MQ=255)
                                                    aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAggg                                                 <  1:494708/100‑1 (MQ=255)
                                                    aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAggg                                                 >  1:32160/1‑100 (MQ=255)
                                                      gAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                               <  2:2005053/100‑1 (MQ=255)
                                                          aTTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTaa                                           <  1:229832/100‑1 (MQ=255)
                                                           ttttCCTTATCAGGAAATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                          >  1:6901/1‑100 (MQ=255)
                                                                 ttATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                    <  2:2604257/100‑1 (MQ=255)
                                                                 ttATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                    <  2:3478849/100‑1 (MQ=255)
                                                                 ttATCAGGAGATCCTTCTGACTCGTCTTGGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACgag                                       <  2:2986799/97‑1 (MQ=255)
                                                                  tATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                    >  1:969008/1‑99 (MQ=255)
                                                                   aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGACGGCTCAGGGGATTAACGAGACgtt                                  >  1:609341/1‑100 (MQ=255)
                                                                    tgaGGAGATCTTTCTGACTCGTCTTCGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                 <  2:343603/98‑1 (MQ=255)
                                                                     cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                 >  1:1379787/1‑99 (MQ=255)
                                                                       ggAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGttt                              <  1:1007360/100‑1 (MQ=255)
                                                                        gagaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                             <  1:3321386/100‑1 (MQ=255)
                                                                           aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTGACGAGACGTTGTTTATgg                          <  1:3281610/100‑1 (MQ=255)
                                                                           aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATgg                          <  2:1790137/100‑1 (MQ=255)
                                                                           aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATgg                          <  1:3356637/100‑1 (MQ=255)
                                                                            tCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                         <  1:3566129/100‑1 (MQ=255)
                                                                            tCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                         <  2:1082002/100‑1 (MQ=255)
                                                                                tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                     <  2:848423/100‑1 (MQ=255)
                                                                                tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGAAAAAAAGAAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                      <  2:418985/99‑1 (MQ=255)
                                                                                  tGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                   >  2:2271444/1‑100 (MQ=255)
                                                                                  tGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                   >  2:186363/1‑100 (MQ=255)
                                                                                   gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAtt                  >  2:754865/1‑100 (MQ=255)
                                                                                    agcttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTa                 >  1:1005377/7‑100 (MQ=39)
                                                                                     gcttgcTTTGCATGTGCATGCAAAGCAAGCTTCTGGAGAACCGCAATAAAATGCTGAAAGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                >  2:2039749/6‑100 (MQ=255)
                                                                                      cttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGCGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg               >  2:1453021/5‑100 (MQ=255)
                                                                                       ttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc              >  2:3136492/4‑100 (MQ=255)
                                                                                        tgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt             >  1:1074081/3‑100 (MQ=255)
                                                                                        gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt             >  2:2989871/1‑100 (MQ=255)
                                                                                         gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg            >  1:2899649/2‑100 (MQ=255)
                                                                                         gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg            >  1:1546797/2‑100 (MQ=255)
                                                                                         gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg            >  1:1633014/2‑100 (MQ=255)
                                                                                          cTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg           >  2:633858/1‑100 (MQ=255)
                                                                                          cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg           <  1:3286320/100‑1 (MQ=255)
                                                                                          cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg           >  1:1411785/1‑100 (MQ=255)
                                                                                           tttGCATGCACATGCAAAGCAAGCTGCTGGAGATCCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa          <  1:2989165/100‑1 (MQ=255)
                                                                                            ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg         <  1:2271459/100‑1 (MQ=255)
                                                                                                 tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa    <  2:2381828/100‑1 (MQ=255)
                                                                                                 tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa    <  2:3392237/100‑1 (MQ=255)
                                                                                                 tgcacaggtaAAGCAAGCTGCTGGAGAACTGCAACAGAAAGCTGAAGTCTCAGTAGATTATCGAGGCGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa    <  2:1708246/100‑1 (MQ=255)
                                                                                                  gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGCGACGTTGTTTATGGCGTTGATTACGCTGGAGTATCAg   >  2:1135331/1‑100 (MQ=255)
                                                                                                  gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAg   <  2:1076545/100‑1 (MQ=255)
                                                                                                   cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg  >  1:388747/1‑100 (MQ=255)
                                                                                                   cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg  >  1:1785692/1‑100 (MQ=255)
                                                                                                   cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg  >  2:92845/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGG  >  NC_000913/2810773‑2810971

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: