Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 15.9% intergenic (‑150/+17) 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,3230AT15.9% 441.9 / 73.8 182intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (66/87);  new base T (9/20);  total (75/107)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

Only 100 of 187 total aligned reads displayed.

TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGA  >  NC_000913/3067224‑3067421
                                                                                                   |                                                                                                  
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                    <  2:407194/100‑1 (MQ=255)
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                    >  1:2500549/1‑100 (MQ=255)
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                    >  2:59928/1‑100 (MQ=255)
tgatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                    >  1:678343/1‑100 (MQ=255)
 gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATAc                                                                                                   >  1:1273695/1‑100 (MQ=255)
     aaaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                               <  1:3303095/100‑1 (MQ=255)
     aaaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                               <  2:3512732/100‑1 (MQ=255)
      aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCTATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                              <  1:2358863/100‑1 (MQ=255)
      aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                              <  1:1419017/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                             <  2:4100317/99‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                            <  2:2370084/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                            <  1:309477/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                            <  2:1082130/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                            <  2:720889/100‑1 (MQ=255)
        ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCggg                                                                                             >  2:1909312/1‑99 (MQ=255)
        ccGCGAGCAAAGATAGGTGATGACGTCATGCTTTTATAGAAAATTAGAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                             <  1:2973247/99‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                           <  2:3021790/100‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                           <  2:2218906/100‑1 (MQ=255)
         cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                           <  1:622250/100‑1 (MQ=255)
          gcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                            <  2:4169035/98‑1 (MQ=255)
           cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                         >  2:4374431/1‑100 (MQ=255)
            gAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaat                                                                                        <  2:785088/100‑3 (MQ=255)
            gAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACa                                                                                        <  2:311044/100‑1 (MQ=255)
             aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaatg                                                                                       <  1:229460/100‑4 (MQ=39)
             aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaatg                                                                                       <  2:3172281/100‑4 (MQ=39)
               cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCACGACGCATACGTAGGCATACaca                                                                                      >  2:1907832/1‑97 (MQ=255)
                aaaGATAGGTGATTACGTCATGTTTTTACATCAAATTACATATAAAGGAGGCAATATCGGGTAAAGGCAGTAGCCCGACGTATTCGTCGGGCTaatgcc                                                                                     <  1:1130738/99‑6 (MQ=25)
                aaaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAata                                                                                    <  2:873012/100‑1 (MQ=255)
                    atgggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtattat                                                                                <  1:1541037/97‑1 (MQ=255)
                     tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                               >  2:2527532/1‑100 (MQ=255)
                     tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATattatt                                                                               <  2:2811123/100‑1 (MQ=255)
                      aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                              <  1:2745480/100‑1 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                           >  2:3893170/1‑100 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                           <  1:3785168/100‑1 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                           >  1:2066699/1‑100 (MQ=255)
                         tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                           <  2:1853400/100‑1 (MQ=255)
                             tACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCa                                                                       <  1:786280/100‑1 (MQ=255)
                              aCGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                                                      <  1:3148238/100‑1 (MQ=255)
                                   aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                                 <  1:2543762/100‑1 (MQ=255)
                                   aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                                 <  1:2134423/100‑1 (MQ=255)
                                    tGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtgtag                                                                >  1:1589317/1‑98 (MQ=255)
                                      gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                                              >  2:147899/1‑100 (MQ=255)
                                      gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                                              >  2:4017289/1‑100 (MQ=255)
                                         ttacagaaaatttcagaaaaAGGAGGGAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCTGGCTAGAGATATTATTGTTCTGGAGGTGGGTTAGCg                                                           >  1:1300889/1‑100 (MQ=255)
                                         ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                           <  1:4097180/100‑1 (MQ=255)
                                         ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                           <  1:1694155/100‑1 (MQ=255)
                                         ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                           <  1:1207599/100‑1 (MQ=255)
                                         gatctgaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                           <  1:3827322/95‑1 (MQ=255)
                                           acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAGATTATTGTGCTGCAGGTGTTTTAGCggg                                                         <  1:1002463/100‑1 (MQ=255)
                                            cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                        >  2:1269143/1‑100 (MQ=255)
                                              gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttg                                                      >  2:2544088/1‑100 (MQ=255)
                                                aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                    <  1:1827372/100‑1 (MQ=255)
                                                aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATCGTGCTGCAGGTGTTTTAGCGGgttgtt                                                    <  2:1569276/100‑1 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                   <  1:3974266/100‑1 (MQ=255)
                                                 aattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGATACAAATATTATTGTGATGCAGGTGTTTTAGCGAGttattg                                                   >  1:1012882/1‑100 (MQ=255)
                                                   ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAATTATTATTGTGCTGCAGGTGTTTTAGAGGGTGGTTGAt                                                 >  1:3555967/1‑100 (MQ=255)
                                                   ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                 >  2:1392244/1‑100 (MQ=255)
                                                   ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                 <  1:584473/100‑1 (MQ=255)
                                                   ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                 <  2:4055657/100‑1 (MQ=255)
                                                     acagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                               <  2:1641066/100‑1 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                               <  1:461183/99‑1 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                              >  2:757290/1‑100 (MQ=255)
                                                      cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                              >  2:4234087/1‑100 (MQ=255)
                                                       agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac                                             <  2:3279677/100‑1 (MQ=255)
                                                       agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac                                             >  1:2695008/1‑100 (MQ=255)
                                                        gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                            >  1:3342839/1‑100 (MQ=255)
                                                           aaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGt                                         >  1:358699/1‑100 (MQ=255)
                                                             aggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTc                                       >  1:265846/1‑100 (MQ=255)
                                                              ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                      <  1:537338/100‑1 (MQ=255)
                                                               gaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTa                                     >  1:3404451/1‑100 (MQ=255)
                                                                  gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGTTTCTAAc                                   <  1:1854135/99‑1 (MQ=255)
                                                                         cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                           <  2:601490/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                          <  2:2560736/100‑1 (MQ=255)
                                                                          gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                          <  2:2852209/100‑1 (MQ=255)
                                                                           ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAcc                          >  2:1713341/1‑99 (MQ=255)
                                                                              aatatttgtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcacct                      >  1:2679909/9‑98 (MQ=37)
                                                                                aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCg                    >  2:895127/1‑100 (MQ=255)
                                                                                 ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                   >  1:1731568/1‑100 (MQ=255)
                                                                                 ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                   >  1:3006759/1‑100 (MQ=255)
                                                                                 atttgtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATcta                   >  1:4252444/6‑98 (MQ=38)
                                                                                  gCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGa                   <  1:3898151/99‑1 (MQ=255)
                                                                                   ttgtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                 >  2:1174107/4‑100 (MQ=39)
                                                                                   cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                 >  2:3546637/1‑100 (MQ=255)
                                                                                     ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg               <  1:3449240/100‑1 (MQ=255)
                                                                                      tAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa              >  1:4331697/1‑100 (MQ=255)
                                                                                       aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat             >  2:4028123/1‑100 (MQ=255)
                                                                                       aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat             >  2:743851/1‑100 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTGTTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc            >  2:4248895/1‑100 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCGCATCGACCTGatc            >  2:3740395/1‑100 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc            <  2:2083612/100‑1 (MQ=255)
                                                                                        gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc            <  2:2233797/100‑1 (MQ=255)
                                                                                         cccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatca           >  1:3272053/1‑100 (MQ=255)
                                                                                          ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatcat          <  1:2086527/100‑1 (MQ=255)
                                                                                            gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGAtcatca        <  2:2425762/100‑1 (MQ=255)
                                                                                              cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa      >  2:3792372/1‑100 (MQ=255)
                                                                                              cgaATACGTCGGGCTACAAATATTATTATGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa      >  1:3891418/1‑100 (MQ=255)
                                                                                               gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGTGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc     <  2:2109209/100‑1 (MQ=255)
                                                                                               gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc     <  2:1731575/100‑1 (MQ=255)
                                                                                               gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAAc     <  2:1492373/100‑1 (MQ=255)
                                                                                                 aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACACGTTCTAACTGGAAGACCACATCGCCCTGATCATCaacc     >  1:1969180/1‑96 (MQ=255)
                                                                                                   |                                                                                                  
TGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGA  >  NC_000913/3067224‑3067421

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: