Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,824 T→G L250R (CTG→CGG)  yigM → putative inner membrane EamA‑like transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,8240TG92.3% 550.3 / 24.3 183L250R (CTG→CGG) yigMputative inner membrane EamA‑like transporter
Reads supporting (aligned to +/- strand):  ref base T (10/4);  new base G (97/72);  total (107/76)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 186 total aligned reads displayed.

GGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGC  >  NC_000913/4011727‑4011922
                                                                                                 |                                                                                                  
ggCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCggg                                                                                                  >  1:3747759/1‑100 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGc                                                                                                 <  1:3993150/100‑1 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGc                                                                                                 >  1:1579404/1‑100 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGc                                                                                                 >  1:1647289/1‑100 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGc                                                                                                 >  1:3574891/1‑100 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGc                                                                                                 >  1:2903539/1‑100 (MQ=255)
 gCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGc                                                                                                 >  1:1895592/1‑100 (MQ=255)
  cTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCt                                                                                                >  1:739641/1‑100 (MQ=255)
     cTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATc                                                                                             <  2:3547784/100‑1 (MQ=255)
      tttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCt                                                                                            >  1:2252279/1‑100 (MQ=255)
       ttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTg                                                                                           <  1:1390419/100‑1 (MQ=255)
       ttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTg                                                                                           <  2:3676364/100‑1 (MQ=255)
         aTGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGc                                                                                         <  1:1162726/100‑1 (MQ=255)
          tgtgGAACTACGGCGCGACACAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCa                                                                                        >  2:3655358/1‑100 (MQ=255)
           gtgGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCAc                                                                                       <  1:3078034/100‑1 (MQ=255)
           gtgGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCAc                                                                                       >  2:1882738/1‑100 (MQ=255)
           gtgGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCAc                                                                                       >  1:2509017/1‑100 (MQ=255)
            tgGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCAcc                                                                                      <  2:586901/100‑1 (MQ=255)
             ggAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCa                                                                                     >  1:3226270/1‑100 (MQ=255)
                aCTATGGCGCGACGTAGGTGGACGCCTTAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACa                                                                                  <  2:1997305/100‑1 (MQ=255)
                 gtACTGCGCGACGCAGGTGGTCGCCGGAACGCTGGGCATTATGAATAATAGGCACGTTCCAGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACac                                                                                 >  2:1759133/2‑99 (MQ=255)
                 cTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAg                                                                                 >  2:2613108/1‑100 (MQ=255)
                 cTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAg                                                                                 >  2:2243067/1‑100 (MQ=255)
                  tACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGc                                                                                <  1:516503/100‑1 (MQ=255)
                    cGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCg                                                                              >  1:2411788/1‑100 (MQ=255)
                    cGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCg                                                                              <  2:2976898/100‑1 (MQ=255)
                     ggCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGc                                                                             >  1:4046335/1‑100 (MQ=255)
                          gACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTgg                                                                        >  1:2313266/1‑100 (MQ=255)
                           aCGCAGGTGGAGGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGc                                                                       <  2:478932/100‑1 (MQ=255)
                           aCGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGc                                                                       >  2:3950256/1‑100 (MQ=255)
                           aCGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGc                                                                       <  1:3160893/100‑1 (MQ=255)
                              cAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCaa                                                                    >  2:1088378/1‑100 (MQ=255)
                              cAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCaa                                                                    >  2:1355299/1‑100 (MQ=255)
                                ggTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACg                                                                  >  2:1282702/1‑100 (MQ=255)
                                ggTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACg                                                                  >  2:2953274/1‑100 (MQ=255)
                                ggTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACg                                                                  <  1:608946/100‑1 (MQ=255)
                                    gACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTtta                                                              >  1:3822856/1‑100 (MQ=255)
                                    gACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTtta                                                              >  2:758659/1‑100 (MQ=255)
                                     aCGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttat                                                             <  2:430982/100‑1 (MQ=255)
                                     aCGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttat                                                             >  1:2772205/1‑100 (MQ=255)
                                      cGCCGGAACGCTGTGCATTATGAATAAGATGCACGTTCCGGCAGGGCTGCTTGGAAACCGGGCTATCTGGCACCAACAGCCGCACGGGCCAACGTttatt                                                            >  2:2981882/1‑100 (MQ=255)
                                      cGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatt                                                            >  1:1004252/1‑100 (MQ=255)
                                      cGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatt                                                            >  2:349500/1‑100 (MQ=255)
                                      cGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatt                                                            <  1:4014926/100‑1 (MQ=255)
                                        ccGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                          >  2:1063976/1‑100 (MQ=255)
                                        ccGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                          <  1:784045/100‑1 (MQ=255)
                                        ccGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                          >  1:2164989/1‑100 (MQ=255)
                                        ccGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                          >  1:2419073/1‑100 (MQ=255)
                                        ccGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTAc                                                          <  1:2846241/100‑1 (MQ=255)
                                             aCGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGgcg                                                     >  1:1450019/1‑100 (MQ=255)
                                               gCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGgcgcg                                                   <  2:3376502/100‑1 (MQ=255)
                                                cTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGgcgcgc                                                  >  2:2604183/1‑100 (MQ=255)
                                                  gggCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTg                                                >  2:1511678/1‑100 (MQ=255)
                                                   ggCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTgg                                               <  2:61361/100‑1 (MQ=255)
                                                      aTTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGa                                            >  2:2617908/1‑100 (MQ=255)
                                                       ttATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTgct                                           >  2:3602890/1‑98 (MQ=255)
                                                       ttATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGAt                                           <  1:2960185/100‑1 (MQ=255)
                                                        tATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATc                                          <  1:3526471/100‑1 (MQ=255)
                                                         aTGAATACTATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTCGCACCAACCGCCGCACTGGCCAACGTTTATTACAGGCGAGCTCGTGATcc                                         >  2:2380766/1‑100 (MQ=255)
                                                         aTGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATcc                                         <  2:1006784/100‑1 (MQ=255)
                                                         aTGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATcc                                         <  2:2819790/100‑1 (MQ=255)
                                                          tGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCt                                        >  2:2535557/1‑100 (MQ=255)
                                                          tGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCt                                        >  1:3688859/1‑100 (MQ=255)
                                                            aataatATGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTgg                                      <  2:1684510/100‑1 (MQ=255)
                                                             ataataTGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCCCTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGc                                     >  1:1842223/1‑100 (MQ=255)
                                                             ataataTGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGc                                     >  2:3254045/1‑100 (MQ=255)
                                                             ataataTGCACGTTCCGACAGGGCTGCTGGTAAACAGGGCTATCTGGCACCAACCGCCGCACTGGCCAACGTTTATTACACGCGCGCAGCTAATccagac                                     >  2:694085/1‑96 (MQ=255)
                                                              taataTGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGcc                                    <  1:3783712/100‑1 (MQ=255)
                                                              taataTGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGcc                                    <  1:2066453/100‑1 (MQ=255)
                                                              taataTGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGcc                                    <  1:2066294/100‑1 (MQ=255)
                                                                atatGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTc                                  <  2:2413906/100‑1 (MQ=255)
                                                                   tGCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACt                               >  1:1455122/1‑100 (MQ=255)
                                                                    gCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtg                              <  2:1117299/100‑1 (MQ=255)
                                                                    gCACGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCAAACGTTTATTACAGGCGCGCTGGTGATCCTGGCGTCACta                              >  1:3693136/1‑99 (MQ=255)
                                                                      aCGTTCCGGAAGGGCGGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtgtg                            <  2:1608137/100‑1 (MQ=255)
                                                                       cGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTgg                           >  2:812361/1‑100 (MQ=255)
                                                                       cGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTgg                           >  2:1448483/1‑100 (MQ=255)
                                                                       cGTTCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTgg                           >  1:2815225/1‑100 (MQ=255)
                                                                          tCCGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCATGGCCTCACTGTGGGTg                        >  2:3142350/1‑100 (MQ=255)
                                                                           ccGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGc                       >  2:1334589/1‑100 (MQ=255)
                                                                            cGGCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCa                      >  1:1699549/1‑100 (MQ=255)
                                                                             ggCAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCAt                     >  1:2119602/1‑100 (MQ=255)
                                                                               cAGGGCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCg                   >  2:1973474/1‑100 (MQ=255)
                                                                                 gggCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTa                 >  2:3236148/1‑100 (MQ=255)
                                                                                 gggCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTa                 <  2:1230059/100‑1 (MQ=255)
                                                                                  ggCTGCTGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTaa                >  2:1664752/1‑100 (MQ=255)
                                                                                   gctgctGGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAg               >  2:1030715/1‑100 (MQ=255)
                                                                                       ctgGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTggg           >  2:449772/1‑100 (MQ=255)
                                                                                       ctgGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTggg           <  1:1788676/100‑1 (MQ=255)
                                                                                       ctgGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTggg           >  2:2515789/1‑100 (MQ=255)
                                                                                       ctgGTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTggg           >  1:848842/1‑100 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         >  1:1240549/1‑100 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         <  2:833278/100‑1 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         <  1:1538006/100‑1 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         >  1:1946571/1‑100 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         <  2:2705060/100‑1 (MQ=255)
                                                                                         ggTAAACCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc         >  2:2655921/1‑100 (MQ=255)
                                                                                          gTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCg        <  1:2653544/100‑1 (MQ=255)
                                                                                              aCCGGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTcc    >  1:2418035/1‑100 (MQ=255)
                                                                                                cTGGCTATCTGGCACCAACAGCCGCACTGGCCAACTTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCcgc  >  2:3900203/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGC  >  NC_000913/4011727‑4011922

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