Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,320 A→T 16.8% intergenic (‑147/+20) ygfI ← / ← yggE putative DNA‑binding transcriptional regulator/oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3200AT16.8% 519.3 / 91.5 214intergenic (‑147/+20)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (83/95);  new base T (17/19);  total (100/114)
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.09e-01

Only 100 of 228 total aligned reads displayed.

CCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACT  >  NC_000913/3067221‑3067419
                                                                                                   |                                                                                                   
ccTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCcgacga                                                                                                     >  2:1342707/1‑100 (MQ=255)
ccTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCcgacga                                                                                                     <  1:4201497/100‑1 (MQ=255)
   tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                  <  1:2839887/100‑1 (MQ=255)
    gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAACATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATAc                                                                                                 >  1:2803031/1‑100 (MQ=255)
      tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                               >  2:860001/1‑100 (MQ=255)
      tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                               >  2:3664576/1‑100 (MQ=255)
       gaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                                              <  1:4243642/100‑1 (MQ=255)
        aaaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATtcgtcg                                                                                             <  2:1828755/100‑1 (MQ=255)
          aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCggg                                                                                           <  2:2753158/100‑1 (MQ=255)
          aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGCCggg                                                                                           >  1:1042201/1‑100 (MQ=255)
           ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                           <  2:3004312/99‑1 (MQ=255)
           ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                          <  2:2065475/100‑1 (MQ=255)
           ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                          <  2:2428427/100‑1 (MQ=255)
           ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                          <  2:781972/100‑1 (MQ=255)
           ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATAGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                          >  1:2682011/1‑100 (MQ=255)
            cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                          <  1:2926444/99‑1 (MQ=255)
            cgcgAGCAAAGATAGGTGATTACGTCAGGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                          <  1:1512263/99‑1 (MQ=255)
              cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                       <  2:1082898/100‑1 (MQ=255)
                aGCAAGGATAGGTGAGTACGTCATGGGTTTATAGACAATTTCAGAAGAAGGAGGCAATATCGGGTAATGGCATTAGCCCTACGAATACGTCTGGc                                                                                          <  1:2384731/95‑1 (MQ=255)
                aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaatg                                                                                     <  2:3977330/100‑4 (MQ=39)
                    aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAatat                                                                                 <  2:2810796/100‑1 (MQ=255)
                    aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAatat                                                                                 >  2:3379312/1‑100 (MQ=255)
                       aTAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtattat                                                                              >  2:2482237/1‑100 (MQ=255)
                        tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                             >  1:1759651/1‑100 (MQ=255)
                         aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                            >  2:2834445/1‑100 (MQ=255)
                           gTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgtg                                                                          <  2:1635649/100‑1 (MQ=255)
                            tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                         >  2:190065/1‑100 (MQ=255)
                              aTTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctg                                                                       >  2:767516/1‑100 (MQ=255)
                                tACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCa                                                                     <  2:2577421/100‑1 (MQ=255)
                                  cGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAgg                                                                   >  1:2172884/1‑100 (MQ=255)
                                   gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCGATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                                                  >  2:3099443/1‑100 (MQ=255)
                                    tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                                                 <  1:3722808/100‑1 (MQ=255)
                                    tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCCGCAGgtg                                                                 <  2:3280829/100‑1 (MQ=255)
                                      attgTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                               <  2:241054/97‑1 (MQ=255)
                                         gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                                            >  2:3520333/1‑100 (MQ=255)
                                         gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGCTTTa                                                            >  1:2010943/1‑100 (MQ=255)
                                            ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                         <  2:4056628/100‑1 (MQ=255)
                                              acagaaaattacagaacaAGGAGGCAATATCGGGTAAAGTCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                       <  1:52111/100‑1 (MQ=255)
                                               cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                      >  1:2469810/1‑100 (MQ=255)
                                                agaaacttacagaaaaAGGAGGCAATATCGGGTAAAGTCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgggn                                                      <  1:1181522/99‑2 (MQ=255)
                                                 gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                    >  2:456144/1‑100 (MQ=255)
                                                 gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                    >  2:394833/1‑100 (MQ=255)
                                                 gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                    >  2:1686178/1‑100 (MQ=255)
                                                 gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                    >  2:146977/1‑100 (MQ=255)
                                                  aaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                   <  2:2635542/100‑1 (MQ=255)
                                                    aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                 <  1:762497/100‑1 (MQ=255)
                                                    aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                 <  1:739225/100‑1 (MQ=255)
                                                    aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                 <  2:1330036/100‑1 (MQ=255)
                                                      ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGGTGCAGGGGGTTTAGGGGGTTGTTTAt                                               >  1:3848137/1‑100 (MQ=255)
                                                      ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGGTTTAGCGGGTTGTTGAt                                               >  2:1965268/1‑100 (MQ=255)
                                                       tacagaaaaAGGAGGCACTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                                              <  1:3507324/100‑1 (MQ=255)
                                                        acagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                             <  2:3609337/100‑1 (MQ=255)
                                                          agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac                                           >  2:2381670/1‑100 (MQ=255)
                                                           gaaaaAGGAGGCAATATCGGGTAAATGCATTAGCCCGACGAATACGTCGGGCTACAATTATTATTGTGCTGCAGGTGTTTTGGCGGGTTGTTGATcc                                             <  1:947338/97‑1 (MQ=255)
                                                           gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                          <  1:1262342/100‑1 (MQ=255)
                                                            aaaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                          <  2:1759657/99‑1 (MQ=255)
                                                              aaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                        >  1:3105250/1‑99 (MQ=255)
                                                              aaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGt                                       <  2:554446/100‑1 (MQ=255)
                                                               aaGGAGGCAATATCGGGTAAAGGCATTAGTCCGACGAATACGTCGGGCTACAAATATTATTGTGCTCCAGGTGTTTTAGCGGGTTGTTGATCCACAGGtt                                      <  1:2643272/100‑1 (MQ=255)
                                                               aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGt                                       <  2:1663/99‑1 (MQ=255)
                                                                 ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                    <  2:3756212/100‑1 (MQ=255)
                                                                 ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                    <  2:4308689/100‑1 (MQ=255)
                                                                  gaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTa                                   <  1:2833262/100‑1 (MQ=255)
                                                                     gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                <  1:860002/100‑1 (MQ=255)
                                                                     gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                >  2:3808259/1‑100 (MQ=255)
                                                                     gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                <  1:1619754/100‑1 (MQ=255)
                                                                      cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                               <  2:4219772/100‑1 (MQ=255)
                                                                          atCGGGTAAAGGCATTAGCCCGCCGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                           >  1:2554019/1‑100 (MQ=255)
                                                                          atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACACATATCATCGCGCTGCACGTGTATTAGCGGGATGTTCACCCACAGGTTCTAACTAGac                            <  1:1742545/99‑2 (MQ=255)
                                                                          atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                           <  2:602223/100‑1 (MQ=255)
                                                                          atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                           >  1:1468251/1‑100 (MQ=255)
                                                                             gatctaaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                        <  1:3639472/95‑1 (MQ=255)
                                                                              ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACca                       <  2:1497003/100‑1 (MQ=255)
                                                                              ggTAAAAGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGCAGAcc                        >  1:2088846/1‑99 (MQ=255)
                                                                               gTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                      <  2:2005163/100‑1 (MQ=255)
                                                                                 aaaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACAt                    >  2:1020269/1‑100 (MQ=255)
                                                                                  aaGGCATTAGCCCGACGAATACGTCGGGCTACAAATCTTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATc                   <  1:49056/100‑1 (MQ=255)
                                                                                  aaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATc                   <  2:4089003/100‑1 (MQ=255)
                                                                                  aaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATc                   >  2:2957187/1‑100 (MQ=255)
                                                                                  aaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGATTTGCCGGGTTGTTGATACTCAGGTTCTAACTGCAAGAACACATc                   >  1:4383756/1‑100 (MQ=255)
                                                                                   aGGCATTAGCCCGACGAATACGTCGGGCTACCAATATTATTGTGCTGCAGGTGTCTTAGCGGGCTGTTGATCCCCAGTTTCTCACTGGCAGACcaca                     >  1:1344731/1‑97 (MQ=255)
                                                                                   aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTCGATGCACAGGTTCTAACTGGGAGACcacagg                   >  1:3592511/1‑97 (MQ=255)
                                                                                    ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                 >  2:3461193/1‑100 (MQ=255)
                                                                                     tttgtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAc                >  2:2227464/5‑100 (MQ=39)
                                                                                          aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat           <  1:295472/100‑1 (MQ=255)
                                                                                          aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat           >  1:3528111/1‑100 (MQ=255)
                                                                                           gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTGTAACTGGAAGACCACATCGACCTGatc          >  2:2022235/1‑100 (MQ=255)
                                                                                           gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  1:3588163/1‑100 (MQ=255)
                                                                                           gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  1:2213057/1‑100 (MQ=255)
                                                                                           gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATTCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  2:566630/1‑100 (MQ=255)
                                                                                           gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  2:2301718/1‑100 (MQ=255)
                                                                                           gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  1:3535872/1‑100 (MQ=255)
                                                                                           gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc          >  1:984185/1‑100 (MQ=255)
                                                                                            cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca         <  2:4237414/100‑1 (MQ=255)
                                                                                            cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca         <  2:1440732/100‑1 (MQ=255)
                                                                                              cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatcatc       >  1:615921/1‑100 (MQ=255)
                                                                                              cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatcatc       >  2:1346439/1‑100 (MQ=255)
                                                                                                  gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc   <  2:1378602/100‑1 (MQ=255)
                                                                                                   aaTACGTCGGGGTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCAAACTGGAAGACCACATCGACCTGATCATCAAACt  >  1:2313748/1‑100 (MQ=255)
                                                                                                   aaTACGTCAGGCTACAAATATTATTGTGCTGCAGGTGTGTTAGCGGCTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc   <  1:156903/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACT  >  NC_000913/3067221‑3067419

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: