Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,585,150 Δ1 bp 18.8% coding (1337/2280 nt) ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,585,1490G.18.8% 423.1 / 101.8 149coding (1338/2280 nt)ydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (48/72);  new base . (20/8);  total (68/81)
Fisher's exact test for biased strand distribution p-value = 3.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

Only 100 of 153 total aligned reads displayed.

CGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACC  >  NC_000913/1585050‑1585244
                                                                                                   |                                                                                                
cGGTACAGCGCTCGCTTCCCGAACTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTg                                                                                                  >  2:2431143/1‑100 (MQ=255)
 ggTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATACAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTgg                                                                                                 <  1:1125122/100‑1 (MQ=255)
    aCAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTG‑CCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                             <  1:3519370/100‑1 (MQ=255)
     cAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                             <  1:2371340/100‑1 (MQ=255)
      aGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAtt                                                                                            >  1:2321510/1‑100 (MQ=255)
           tcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                      <  1:3618022/100‑1 (MQ=255)
           tcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                      <  1:1515893/100‑1 (MQ=255)
            cgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagca                                                                                     <  1:1687238/100‑1 (MQ=255)
            cgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagca                                                                                     >  1:3294802/1‑100 (MQ=255)
             gctTCCCGATCGGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cag                                                                                       <  1:3994406/97‑1 (MQ=255)
              ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                      >  2:3706650/1‑97 (MQ=255)
              ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATg                                                                                   <  1:2717858/100‑1 (MQ=255)
              ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATg                                                                                   <  2:1711704/100‑1 (MQ=255)
               ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTTCCCCGTACATATTGCTTCCATGCTGGCAATTg‑cagc                                                                                      >  2:4173354/1‑96 (MQ=255)
               ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGt                                                                                  <  2:906092/100‑1 (MQ=255)
               ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTg‑ccagca                                                                                    <  2:1927359/97‑6 (MQ=255)
                tCCCCATCCGGCAGTGCCAGCGCAAAATTGCCCCCCACGGAGATCAATTCTCGAGCCCGCCCCGTACATATTGCTTGCTTGCT‑GCAAGTGCCAGCATgc                                                                                  <  2:2565539/99‑2 (MQ=255)
                  ccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                               <  1:415416/100‑1 (MQ=255)
                  ccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                               >  2:4157488/1‑100 (MQ=255)
                   cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaagc                                                                              <  1:1290209/100‑6 (MQ=38)
                   cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCAAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaan                                                                               <  2:3974558/99‑5 (MQ=255)
                    gATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAgg                                                                             <  1:1566716/100‑1 (MQ=255)
                       cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTg                                                                           <  2:1265170/99‑1 (MQ=255)
                         ggCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAt                                                                        <  1:3581439/100‑1 (MQ=255)
                              tGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGt                                                                   <  2:1170107/100‑1 (MQ=255)
                               gCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTg                                                                  >  2:1541918/1‑100 (MQ=255)
                                ccAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgg                                                                 >  1:499067/1‑100 (MQ=255)
                                 cAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                <  1:982133/100‑1 (MQ=255)
                                  aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                               <  1:638952/100‑1 (MQ=255)
                                  aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                               <  2:3755707/100‑1 (MQ=255)
                                   gcgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGt                                                              >  1:2892091/1‑100 (MQ=255)
                                    cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                             <  2:1354850/100‑1 (MQ=255)
                                    cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                             <  1:3330300/100‑1 (MQ=255)
                                    cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                             <  1:1019055/100‑1 (MQ=255)
                                     gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGa                                                            >  1:3169206/1‑100 (MQ=255)
                                       aaaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTGGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                          <  2:753977/100‑1 (MQ=255)
                                        aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                         >  1:3434822/1‑100 (MQ=255)
                                        aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                         <  1:306105/100‑1 (MQ=255)
                                         aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGcc                                                        <  2:800982/100‑1 (MQ=255)
                                         aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGcc                                                        <  1:1036251/100‑1 (MQ=255)
                                            tGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                     <  2:209718/100‑1 (MQ=255)
                                              ccccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgc                                                   <  2:1554089/100‑1 (MQ=255)
                                               cccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                                  >  2:3295622/1‑100 (MQ=255)
                                               cccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGAATGAGGTGGAGTGAAGCCATAgc                                                   >  1:326094/1‑99 (MQ=255)
                                                ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                                 >  2:1060245/1‑100 (MQ=255)
                                                 cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCt                                                >  2:2768113/1‑100 (MQ=255)
                                                  ccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTc                                               >  1:535361/1‑100 (MQ=255)
                                                   cATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCa                                              >  2:2633386/1‑100 (MQ=255)
                                                    aTGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGCCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                             <  2:2005750/100‑1 (MQ=255)
                                                    aTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                             <  2:499066/100‑1 (MQ=255)
                                                     tGCAGATCAATGCTCGAACCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                            <  1:2081506/100‑1 (MQ=255)
                                                      gCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAccc                                           <  2:2766624/100‑1 (MQ=255)
                                                         gATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                        <  2:3192413/100‑1 (MQ=255)
                                                         gATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                        >  2:1157713/1‑100 (MQ=255)
                                                           tCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACg                                      >  2:3204562/1‑100 (MQ=255)
                                                                gCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                 <  2:2734044/100‑1 (MQ=255)
                                                                gCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                 <  1:754074/100‑1 (MQ=255)
                                                                  tCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                               <  2:2892082/100‑1 (MQ=255)
                                                                     aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc                            <  1:2260534/100‑1 (MQ=255)
                                                                     aGCCTGCCCCGTACATAGTGCTGTCATGCGGGCAATGG‑CAGCACGTCCCGGTGCCTGAGGTGGCGTGAAGTCATAGCGCCCACCCAGACGAGCCAGAAAc                            <  2:3809109/100‑1 (MQ=255)
                                                                        cTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                         <  1:3956017/100‑1 (MQ=255)
                                                                        cTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                         <  1:4067034/100‑1 (MQ=255)
                                                                         tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTg                        <  2:179633/100‑1 (MQ=255)
                                                                          gCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCCTGTCCAGGTGCAGGAGGCGGGGTCAAGCCATAGCGCTCACCCAGACGAGCCAGCAACTCCGc                       >  2:2045366/1‑100 (MQ=255)
                                                                          gCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGc                       <  2:292041/100‑1 (MQ=255)
                                                                          gCCCCGTACATATTACTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGc                       <  1:1452656/100‑1 (MQ=255)
                                                                              cGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGAc                   >  1:3540518/1‑100 (MQ=255)
                                                                               gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGGCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg                  <  1:2139319/100‑1 (MQ=255)
                                                                                   atatTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGttt              <  2:2321505/100‑1 (MQ=255)
                                                                                    tatTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtc             >  2:3271805/1‑100 (MQ=255)
                                                                                    tatTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtc             >  2:1048006/1‑100 (MQ=255)
                                                                                     atTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct            >  1:3041364/1‑100 (MQ=255)
                                                                                      ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtctc           <  2:3540501/100‑1 (MQ=255)
                                                                                       cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg          >  2:2784473/7‑100 (MQ=38)
                                                                                       cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg          >  1:321461/7‑100 (MQ=38)
                                                                                       cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg          >  1:3413958/7‑100 (MQ=38)
                                                                                       cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg          >  1:1526733/7‑100 (MQ=38)
                                                                                        gcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg         <  1:2221528/100‑1 (MQ=255)
                                                                                        ctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg         >  1:1858305/6‑100 (MQ=38)
                                                                                        ctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg         >  1:1991210/6‑100 (MQ=38)
                                                                                        ctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg         >  1:1857986/6‑100 (MQ=38)
                                                                                         tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg         >  2:2465744/5‑99 (MQ=38)
                                                                                         tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  2:4118002/5‑100 (MQ=255)
                                                                                         tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  1:965711/5‑100 (MQ=255)
                                                                                         tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  1:373805/5‑100 (MQ=255)
                                                                                         cttgcGTGCTGGCAACTG‑CAGCGTGTCCAGGTGCATCAGGTGGGGTGACGCCATAGCGCTTGCCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        <  2:473167/100‑1 (MQ=255)
                                                                                         cttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  2:2398128/1‑100 (MQ=255)
                                                                                         cttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  2:3669647/1‑100 (MQ=255)
                                                                                         cttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt        >  1:1551464/1‑100 (MQ=255)
                                                                                          ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg       >  1:383874/1‑100 (MQ=255)
                                                                                          ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg       <  1:2405164/100‑1 (MQ=255)
                                                                                          ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg       <  1:2465382/100‑1 (MQ=255)
                                                                                           tgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa      <  1:507116/100‑1 (MQ=255)
                                                                                           gacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa      >  2:2456730/3‑100 (MQ=255)
                                                                                            gcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGAt     <  1:1667148/100‑1 (MQ=255)
                                                                                             cATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATa    <  1:2419064/100‑1 (MQ=255)
                                                                                             cATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTAACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATa    >  2:693535/1‑100 (MQ=255)
                                                                                             cATGATGTCAATTT‑CAGCAGGTCCAAGTGCATGAGGGTGGGTGAAGACATAGTGCTCTCCCAGACGAGCCAGAACCCCTGCAGACGGTTTCTCGGTg       <  2:681693/97‑1 (MQ=255)
                                                                                              aTGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATAc   <  1:3272343/100‑1 (MQ=255)
                                                                                              aTGCGGGCAATTG‑CAGCACGTCCAGGTGCAGGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg       <  1:1727204/96‑1 (MQ=255)
                                                                                                   |                                                                                                
CGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACC  >  NC_000913/1585050‑1585244

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: ‑