Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,811,274 G→T G169C (GGT→TGT)  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,811,2740GT99.0% 682.0 / ‑4.5 207G169C (GGT→TGT) mprAtranscriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (107/98);  total (108/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

CATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTGG  >  NC_000913/2811178‑2811369
                                                                                                |                                                                                               
 aTCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACtgtg                                                                                              >  2:3228522/1‑99 (MQ=255)
   cAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCTAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACtgtgtg                                                                                            <  1:4477664/99‑1 (MQ=255)
     aCTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGt                                                                                          >  2:4731469/1‑99 (MQ=255)
       tctGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtc                                                                                        >  2:2680680/1‑99 (MQ=255)
       tctGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtc                                                                                        <  1:501180/99‑1 (MQ=255)
       tctGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtc                                                                                        >  2:3991478/1‑99 (MQ=255)
        ctGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtct                                                                                       <  1:4913649/99‑1 (MQ=255)
         tGGTCCGCGCTCAGCACAACAGAAAAATTTCAGCTTGAGCAAATCACGCGCAAATTGCTGTCGCGTCTCGACTTGATGGAACAAGACTGTGTGGTtctc                                                                                      <  2:1886337/99‑1 (MQ=255)
         tGGTCCGCGCTCAGCACAACAGAAAAAGATCATCTCGCGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtctc                                                                                      <  2:4692472/99‑1 (MQ=255)
         tGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTtctc                                                                                      >  1:407112/1‑99 (MQ=255)
           gTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGa                                                                                    >  1:3685770/1‑99 (MQ=255)
            tCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCACATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTCTGGTTCTCGaa                                                                                   <  2:4968208/99‑1 (MQ=255)
             ccGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAg                                                                                  <  2:3410149/99‑1 (MQ=255)
              cgcgCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGc                                                                                 <  1:116430/99‑1 (MQ=255)
                cgcTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGa                                                                               >  2:2762621/1‑99 (MQ=255)
                 gcTCAGAACAACAGAACAAGATCAGCTCGAGCAAATCCCCCGCCATTTGCTCTCCCGTCCCGACCAGATGGAACACGACTGTGTGGATCTCGAAgagac                                                                              <  2:2330496/99‑2 (MQ=255)
                   tCAGCACAACAGAAAAAGATCAGCTCGAGCATATTTCCTGCAAATTGTTTTTCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGa                                                                            <  1:137726/99‑1 (MQ=255)
                   tCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGa                                                                            >  2:1450691/1‑99 (MQ=255)
                    cAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAg                                                                           <  1:2901121/99‑1 (MQ=255)
                    cAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAg                                                                           <  1:4157410/99‑1 (MQ=255)
                    cAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAg                                                                           >  2:3625939/1‑99 (MQ=255)
                     aGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGc                                                                          >  2:424271/1‑99 (MQ=255)
                      gCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCt                                                                         <  2:4990763/99‑1 (MQ=255)
                       cacaACCGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGCAGCGATGAGCTa                                                                        >  2:1119109/1‑99 (MQ=255)
                          aacaGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAAcg                                                                     >  1:77715/1‑99 (MQ=255)
                          aacaGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAAcg                                                                     >  1:4411986/1‑99 (MQ=255)
                          aacaGAAAAAGATCAGCTCGAGCAAATCACCAGCAAATTGGTCTCCAGTCTCGACCAGATGGAAGAAGACTGTGTGGTTCTCGAAGCGATGAGCTAAcg                                                                     <  2:3178885/99‑1 (MQ=255)
                           acaGAAAAAGATCAGTTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAAcgc                                                                    <  1:5278962/99‑1 (MQ=255)
                           acaGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGAATGTGTGGTTCTCGAAGCGATGAGCTAActc                                                                    >  2:3093828/1‑97 (MQ=255)
                             aGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGt                                                                  >  1:4791940/1‑99 (MQ=255)
                             aGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGt                                                                  <  2:1547657/99‑1 (MQ=255)
                             aGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGt                                                                  <  1:3066206/99‑1 (MQ=255)
                              gAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTc                                                                 >  1:1987123/1‑99 (MQ=255)
                              gAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTc                                                                 >  1:1797769/1‑99 (MQ=255)
                               aaaaaGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCa                                                                <  1:3280331/99‑1 (MQ=255)
                               aaaaaGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCa                                                                >  1:3101447/1‑99 (MQ=255)
                               aaaaaGATCAGCTCGAGCAAATCACCCACAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCa                                                                >  2:1368409/1‑99 (MQ=255)
                                 aaaGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCAtc                                                              >  1:374885/1‑99 (MQ=255)
                                 aaaGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCAtc                                                              <  2:4443564/99‑1 (MQ=255)
                                  aaGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCAtct                                                             <  2:4621754/99‑1 (MQ=255)
                                     aTCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGc                                                          >  1:2792558/1‑99 (MQ=255)
                                      tCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCt                                                         <  2:76176/99‑1 (MQ=255)
                                      tCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCt                                                         <  1:4258046/99‑1 (MQ=255)
                                        aGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCa                                                       >  2:2863571/1‑99 (MQ=255)
                                          cTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaa                                                     >  1:802055/1‑99 (MQ=255)
                                          cTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaa                                                     >  1:894278/1‑99 (MQ=255)
                                              aGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATc                                                 <  2:3198899/99‑1 (MQ=255)
                                                cAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCa                                               >  2:566749/1‑99 (MQ=255)
                                                cAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCa                                               >  2:5120599/1‑99 (MQ=255)
                                                  aaTCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGa                                             <  2:3816871/99‑1 (MQ=255)
                                                   aTCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGAt                                            <  1:2403064/99‑1 (MQ=255)
                                                   aTCACCCGCAAATTGCCCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGAt                                            <  2:4326448/99‑1 (MQ=255)
                                                      aCCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTta                                         <  1:3039528/99‑1 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  1:2765455/1‑99 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  1:1311481/1‑99 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  2:2900421/1‑99 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  1:5184578/1‑99 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  1:166214/1‑99 (MQ=255)
                                                        ccGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtata                                       >  2:507987/1‑99 (MQ=255)
                                                         cGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATaa                                      <  1:1201285/99‑1 (MQ=255)
                                                          gCAAATTGCTCTCCCGGCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATaaa                                     <  2:3033416/99‑1 (MQ=255)
                                                           ctaaTTGCTCTTTCTTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATaaaa                                    <  1:967567/97‑1 (MQ=255)
                                                              aTTGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGaa                                 <  1:432750/99‑1 (MQ=255)
                                                               ttGCTCTCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGaaa                                <  1:3804057/99‑1 (MQ=255)
                                                                  ctctCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGaaaaaa                             <  1:2825886/99‑1 (MQ=255)
                                                                  ctctCCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGaaaaaa                             <  1:2956031/99‑1 (MQ=255)
                                                                   tctcCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAAt                            <  1:1443304/99‑1 (MQ=255)
                                                                   tctcCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAAt                            >  2:1483985/1‑99 (MQ=255)
                                                                   tctcCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAAt                            >  1:2140227/1‑99 (MQ=255)
                                                                     tcCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGa                          >  2:4159347/1‑99 (MQ=255)
                                                                     tcCCGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGa                          <  2:1992597/99‑1 (MQ=255)
                                                                       ccGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACt                        >  2:4598006/1‑99 (MQ=255)
                                                                       ccGTCTCGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACt                        >  1:385173/1‑99 (MQ=255)
                                                                          tctcGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGc                     <  2:3876418/99‑1 (MQ=255)
                                                                            tcGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGAGTTATAAAAGAAAAAATGACTGGCCa                   >  2:2077222/1‑99 (MQ=255)
                                                                             cGACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAg                  >  1:3924390/1‑99 (MQ=255)
                                                                              gACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGc                 >  1:2281874/1‑99 (MQ=255)
                                                                              gACCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGc                 >  2:2397868/1‑99 (MQ=255)
                                                                              gACCAGATGGAACAAGACTGTGTAGATCTCGATGCGATCCGCAAACGCGACATCACGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGc                 >  2:644054/1‑99 (MQ=255)
                                                                               aCCAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCa                <  1:5023111/99‑1 (MQ=255)
                                                                                ccAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGACCAGCAt               <  1:3146791/99‑1 (MQ=255)
                                                                                 cAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATc              >  2:3761065/1‑99 (MQ=255)
                                                                                 cAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATc              <  1:3594829/99‑1 (MQ=255)
                                                                                 cAGATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATc              <  1:3523685/99‑1 (MQ=255)
                                                                                   gATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGc            >  1:4980883/1‑99 (MQ=255)
                                                                                   gATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGc            >  1:636342/1‑99 (MQ=255)
                                                                                   gATGGAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTGATAAAAGAAAAAATGACTGGCCAGCATCGc            <  1:4375835/99‑1 (MQ=255)
                                                                                      ggAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAAc         >  1:4558887/1‑99 (MQ=255)
                                                                                       gAACACGACTGTTTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTAGCTCAAAAATCCAGTTTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACa        <  2:2999362/99‑1 (MQ=255)
                                                                                       gAACAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACa        <  2:403577/99‑1 (MQ=255)
                                                                                        aaCAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACAt       <  1:865683/99‑1 (MQ=255)
                                                                                        aaCAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACAt       >  2:106003/1‑99 (MQ=255)
                                                                                        aaCAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACAt       >  1:2060689/1‑99 (MQ=255)
                                                                                          cAAGACTGTGTGGTTCTCGAAGCGATGAGCTAATGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGc     <  1:3145295/99‑1 (MQ=255)
                                                                                          cAAGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGc     >  2:719683/1‑99 (MQ=255)
                                                                                           aaGACTGTGTGGTTCTAGAAGCGATGAGCTAACGCGTCATCTCGCTTAAAAATCCAGAAAAAAAAAAAAAAAAATGACTGGCCAGCATCGCAACATGCt    <  2:2376171/99‑1 (MQ=255)
                                                                                            aGACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTg   <  1:3253969/99‑1 (MQ=255)
                                                                                             gACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTgg  >  2:4559514/1‑99 (MQ=255)
                                                                                             gACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTgg  >  1:2610862/1‑99 (MQ=255)
                                                                                             gACTGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTgg  >  2:3151989/1‑99 (MQ=255)
                                                                                                |                                                                                               
CATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTATAAAAGAAAAAATGACTGGCCAGCATCGCAACATGCTGG  >  NC_000913/2811178‑2811369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: