Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,873 A→G 17.6% H35R (CAC→CGC) ‡ 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,8730AG17.6% 277.9 / 53.3 131H35R (CAC→CGC) mprAtranscriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base A (50/58);  new base G (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.07e-01

Only 100 of 161 total aligned reads displayed.

ATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGA  >  NC_000913/2810774‑2810972
                                                                                                   |                                                                                                   
aTAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTtgcatgca                                                                                                     <  1:2317097/100‑1 (MQ=255)
 tAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcac                                                                                                    >  1:3428785/1‑100 (MQ=255)
  aGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcaca                                                                                                   <  2:142757/100‑1 (MQ=255)
   gTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCAAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTTACTCGTCTTTGcatgcacat                                                                                                  >  2:1807316/1‑100 (MQ=255)
    ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatg                                                                                                 <  1:224715/100‑1 (MQ=255)
     tCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgc                                                                                                <  2:1938225/100‑1 (MQ=255)
     tCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgc                                                                                                <  1:507506/100‑1 (MQ=255)
      cGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                               >  2:2269481/1‑100 (MQ=255)
      cGTTTACGCCCATTGAACAAATGCTAAAAGTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                               <  2:1315438/100‑1 (MQ=255)
         ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaag                                                                                            <  2:550936/100‑1 (MQ=255)
         ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaag                                                                                            <  1:3388507/100‑1 (MQ=255)
          tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                           <  2:1284164/100‑2 (MQ=255)
          tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                           <  1:348050/100‑2 (MQ=255)
          tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                           <  2:254175/100‑2 (MQ=255)
          tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagc                                                                                           >  2:2754869/1‑100 (MQ=255)
           aCGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagca                                                                                          >  2:1062420/1‑100 (MQ=255)
           aCGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagca                                                                                          <  1:3305468/100‑1 (MQ=255)
            cGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                         >  2:2609559/1‑100 (MQ=255)
            cGCCCATTGAACAAATGCTAAAATTGCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagacg                                                                                         <  2:198526/100‑4 (MQ=255)
              cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaagc                                                                                       >  2:3587009/1‑100 (MQ=255)
               ccagggAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAGTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAgct                                                                                      <  2:3149790/95‑1 (MQ=255)
                    gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                 <  1:619219/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:574933/1‑100 (MQ=255)
                     aaCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                <  1:1890378/100‑1 (MQ=255)
                       cAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGgaga                                                                              >  1:509877/1‑100 (MQ=255)
                       cAAATGCTAAAATCTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                <  2:3326080/98‑1 (MQ=255)
                         aaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAc                                                                            >  1:1999537/1‑100 (MQ=255)
                           tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                          >  2:3433533/1‑100 (MQ=255)
                           tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                          <  2:1766079/100‑1 (MQ=255)
                           tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                          <  2:454769/100‑1 (MQ=255)
                           tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                          <  2:455739/100‑1 (MQ=255)
                            gatAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                         >  1:2454821/3‑100 (MQ=255)
                             cTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        >  1:1835200/1‑100 (MQ=255)
                             cTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        >  2:1745346/1‑100 (MQ=255)
                             cTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        >  2:303744/1‑100 (MQ=255)
                             cTAAAATATCGTGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        >  2:1910805/1‑100 (MQ=255)
                               aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAAt                                                                      >  1:1227669/1‑100 (MQ=255)
                                       gcgcCAGCCGCCACGACGATTTTCTTTATCAGGAGATCCTTCTGACTCGACTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                               <  2:843763/99‑1 (MQ=255)
                                         gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGa                                                            <  1:1567469/100‑1 (MQ=255)
                                          ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                           <  2:2095051/100‑1 (MQ=255)
                                          ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                           >  1:1339134/1‑100 (MQ=255)
                                            aGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGa                                                            >  2:1595046/1‑97 (MQ=255)
                                             gccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGc                                                        >  2:729282/1‑100 (MQ=255)
                                              ccgccACGAAGATTTTCCTTATCAGGAGCTCCTTCTGACTCGTCTTTGCATGCACATGAAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAgg                                                         >  2:516966/1‑98 (MQ=255)
                                                   aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAggg                                                  >  1:781514/1‑100 (MQ=255)
                                                     gAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                                <  2:2005053/100‑1 (MQ=255)
                                                       aGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAtt                                              >  2:1893832/1‑100 (MQ=255)
                                                         aTTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTaa                                            <  1:229832/100‑1 (MQ=255)
                                                          ttttCCTTATCAGGAAATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                           >  1:6901/1‑100 (MQ=255)
                                                               cTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACgaga                                       <  2:2355400/99‑1 (MQ=255)
                                                                ttATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                     <  1:729285/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)
                                                                   tgaGGAGATCTTTCTGACTCGTCTTCGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                  <  2:343603/98‑1 (MQ=255)
                                                                    cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                  >  1:1379787/1‑99 (MQ=255)
                                                                     aGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                >  2:28177/1‑100 (MQ=255)
                                                                     aGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgtt                                >  1:2271708/1‑100 (MQ=255)
                                                                      ggAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGttt                               <  1:1007360/100‑1 (MQ=255)
                                                                       gagaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                              <  1:3321386/100‑1 (MQ=255)
                                                                          aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTGACGAGACGTTGTTTATgg                           <  1:3281610/100‑1 (MQ=255)
                                                                           tCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                          <  2:1082002/100‑1 (MQ=255)
                                                                           tCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                          <  2:3257787/100‑1 (MQ=255)
                                                                              ttCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                       <  1:454314/100‑1 (MQ=255)
                                                                               tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                      <  2:848423/100‑1 (MQ=255)
                                                                               tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGAAAAAAAGAAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                       <  2:418985/99‑1 (MQ=255)
                                                                                cTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                      <  2:2570732/99‑1 (MQ=255)
                                                                                 gcagcttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                    >  1:2931748/9‑100 (MQ=38)
                                                                                  gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAtt                   >  2:2053817/1‑100 (MQ=255)
                                                                                  gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAtt                   >  2:754865/1‑100 (MQ=255)
                                                                                  gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                    <  1:3276561/99‑1 (MQ=255)
                                                                                   agcttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTa                  >  1:1005377/7‑100 (MQ=39)
                                                                                     cttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg                >  1:2594262/5‑100 (MQ=255)
                                                                                      cGTCTTTGCATGCACATGCAAAGCAAGCTGCTAGAGAAACGCAATAAAATGCTGAAGGCTCAGGGGATTAACAAGACATAGTTTATGACGTTAATTa                  >  2:1525366/1‑97 (MQ=255)
                                                                                       tgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt              >  1:1074081/3‑100 (MQ=255)
                                                                                       gTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt              >  1:2368289/1‑100 (MQ=255)
                                                                                        tCTTTGCATGCACATGCAAACCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTCTTTATGGCGTTGATTACg                >  2:1225912/1‑97 (MQ=255)
                                                                                        gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg             >  1:1633014/2‑100 (MQ=255)
                                                                                        gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg             >  1:2899649/2‑100 (MQ=255)
                                                                                        gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg             >  1:1546797/2‑100 (MQ=255)
                                                                                         cTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg            >  2:633858/1‑100 (MQ=255)
                                                                                         cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg            >  1:3180331/1‑100 (MQ=255)
                                                                                         cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg            >  1:1411785/1‑100 (MQ=255)
                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGATCCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa           <  1:2989165/100‑1 (MQ=255)
                                                                                           ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa           >  1:891422/1‑99 (MQ=255)
                                                                                           ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg          <  1:88212/100‑1 (MQ=255)
                                                                                           ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg          >  1:2170693/1‑100 (MQ=255)
                                                                                                tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa     <  2:3392237/100‑1 (MQ=255)
                                                                                                tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa     >  1:2693123/1‑100 (MQ=255)
                                                                                                tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa     <  2:2381828/100‑1 (MQ=255)
                                                                                                tgcacaggtaAAGCAAGCTGCTGGAGAACTGCAACAGAAAGCTGAAGTCTCAGTAGATTATCGAGGCGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa     <  2:1708246/100‑1 (MQ=255)
                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   >  1:1785692/1‑100 (MQ=255)
                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   <  1:383604/100‑1 (MQ=255)
                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   >  1:388747/1‑100 (MQ=255)
                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   >  2:92845/1‑100 (MQ=255)
                                                                                                  cacacgcaAAGCACGCTGCTGGAGAACCGCCATAAAATGCTGAAGGCTGAGGGGATTAACGAGACGTTGCTTATGGTGTTGATTACGCTGGAGtc        <  2:2513566/95‑1 (MQ=255)
                                                                                                   acatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa  >  1:1692650/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGA  >  NC_000913/2810774‑2810972

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: