Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 18.5% 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,3230AT18.5% 271.3 / 60.7 119intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (49/48);  new base T (9/13);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 120 total aligned reads displayed.

TGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAA  >  NC_000913/3067227‑3067416
                                                                                                |                                                                                             
tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATagat                                                                                            >  2:2794431/1‑97 (MQ=255)
tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                            >  1:2046093/1‑100 (MQ=255)
   aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                                           >  1:1976243/1‑98 (MQ=255)
   aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                          >  1:2739446/1‑99 (MQ=255)
    aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAGTATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                        <  2:1963904/100‑1 (MQ=255)
    aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                        <  2:948362/100‑1 (MQ=255)
    aCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                        <  2:1272602/100‑1 (MQ=255)
     ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                       <  2:1018267/100‑1 (MQ=255)
     ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                       <  2:3418507/100‑1 (MQ=255)
     ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                         >  1:2719860/1‑98 (MQ=255)
     ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                       <  1:118343/100‑1 (MQ=255)
     ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGGAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                       <  1:2293162/100‑1 (MQ=255)
      cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                      <  1:435964/100‑1 (MQ=255)
      cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                      <  1:2518354/100‑1 (MQ=255)
       gcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGTCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                         <  1:1643991/96‑1 (MQ=255)
       gcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTa                                                                                     <  2:2754572/100‑1 (MQ=255)
        cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaa                                                                                    <  2:3166647/100‑2 (MQ=255)
        cgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTAc                                                                                    <  2:2387338/100‑1 (MQ=255)
            cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                >  1:3275052/1‑100 (MQ=255)
            cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                >  1:3194162/1‑100 (MQ=255)
            cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                >  2:1228683/1‑100 (MQ=255)
            cAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                >  1:1439500/1‑100 (MQ=255)
              aggATAGGTGGTTACGTCATGGTTTTAGAGATGATTAAAGAGAACGGAGGTGATATCGGGTATAGGCATTAGCGCGAAGAATAGGTCggg                                                                                        >  1:3243647/3‑90 (MQ=255)
              aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAatat                                                                              >  2:2831406/1‑100 (MQ=255)
               aGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtatt                                                                             <  2:2895375/100‑1 (MQ=255)
                   aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                         >  1:1324103/1‑100 (MQ=255)
                    ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                                                        >  1:1440680/1‑100 (MQ=255)
                    ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                                                        <  2:621074/100‑1 (MQ=255)
                     gTGATTGCGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATAGCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgtg                                                                       <  2:1611967/100‑1 (MQ=255)
                      tGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                      <  1:1754026/100‑1 (MQ=255)
                       gATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgct                                                                     >  2:2752952/1‑100 (MQ=255)
                             gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                                               >  2:149016/1‑100 (MQ=255)
                              tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                                              >  2:2633480/1‑100 (MQ=255)
                              tCATGGTTTTAAAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                                              <  1:378905/100‑1 (MQ=255)
                                 tGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGttt                                                           >  2:1653773/1‑100 (MQ=255)
                                  ggTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                                          <  2:2011838/100‑1 (MQ=255)
                                  ggTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                                          <  2:2765999/100‑1 (MQ=255)
                                   gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGCCGAATACGTCGGGCTACAAATATTATTGTGCTGCCGGTGTTTTa                                                         >  2:2037422/1‑100 (MQ=255)
                                   gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                                         >  1:2146682/1‑100 (MQ=255)
                                    ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                                                        <  2:1294484/100‑1 (MQ=255)
                                     tttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGTTACAAATATTATTGTGCTGCAGGTGTTTTAGc                                                       >  1:1375687/1‑100 (MQ=255)
                                      ttacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                      <  1:43458/100‑1 (MQ=255)
                                        acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                    >  2:1877572/1‑100 (MQ=255)
                                        acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                    >  2:2294833/1‑100 (MQ=255)
                                         cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgt                                                   >  2:1293457/1‑100 (MQ=255)
                                          agaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgtt                                                  >  1:2959530/1‑100 (MQ=255)
                                            aaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                <  1:1445971/100‑1 (MQ=255)
                                            aaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                <  2:234192/100‑1 (MQ=255)
                                             aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                               <  1:2037411/100‑1 (MQ=255)
                                             aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                               <  2:1152617/100‑1 (MQ=255)
                                                ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                            <  1:1450641/100‑1 (MQ=255)
                                                 tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTCGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTCTTCGCGGGTTGTTGATc                                           <  1:1536911/100‑1 (MQ=255)
                                                  acagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                          <  2:2870111/100‑1 (MQ=255)
                                                   cagaaaaAGTATGCAATATCTGGTAAATGCATTAGCCCGACGAATACGTCTTGTTACTATTATTATTTTGTTGCAGGTGTTTTAGCGTGTTGTTGATCca                                         <  1:1230519/100‑1 (MQ=255)
                                                    agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGCCGAATACGTCGGGCTACAAGTATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGACCcac                                        <  1:1864427/100‑1 (MQ=255)
                                                     gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                       <  2:2764786/100‑1 (MQ=255)
                                                       aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCATGTGTTTTAGCGGGTTGTTGATCtac                                        <  2:3132868/97‑1 (MQ=255)
                                                       aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                     <  1:369733/100‑1 (MQ=255)
                                                          aggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTc                                  <  1:1402385/100‑1 (MQ=255)
                                                           ggagtCAATATCGTTCAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATAGTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTc                                  <  1:1995409/99‑1 (MQ=255)
                                                             aggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTaa                               >  2:3390875/1‑100 (MQ=255)
                                                              ggCCATATCGGGAAAAGCCATTAGCTGCACGAATACGTCGGGCTACAAATATTATTCTGCAGCAGGTGTTGTAGCGGGTTGTTGTTCcacat                                      <  1:1377740/92‑2 (MQ=255)
                                                              ggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                              >  2:234041/1‑100 (MQ=255)
                                                              ggCAATATAGGGTAACGGCATTAGCCCGACGAATAGGTCGGACTACAAATAATATTGGGCGGCAGGTGCTTGAGCGGGTTGTTGATCCGCTGGTCCTa                                >  1:3120852/1‑98 (MQ=255)
                                                               gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATCATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc                              <  1:2474462/99‑1 (MQ=255)
                                                                cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                            >  2:326179/1‑100 (MQ=255)
                                                                cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                            >  2:2508850/1‑100 (MQ=255)
                                                                 aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTgg                           <  1:3236905/100‑1 (MQ=255)
                                                                 aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTgg                           >  2:2104518/1‑100 (MQ=255)
                                                                  atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                          >  1:1769735/1‑100 (MQ=255)
                                                                  atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                          <  2:743075/100‑1 (MQ=255)
                                                                   tatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                         <  2:168404/100‑1 (MQ=255)
                                                                    atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                        <  1:406225/100‑1 (MQ=255)
                                                                    atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                        <  2:3011102/100‑1 (MQ=255)
                                                                     tCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                       <  1:1750890/100‑1 (MQ=255)
                                                                     tCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                       <  1:1365361/100‑1 (MQ=255)
                                                                      cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                      <  2:1040048/100‑1 (MQ=255)
                                                                         gTAAAGGCATTAGTCCGACGAATACGTCGGGCTACAAATATTATTGTACTGCAGGTCTTTAAGCGGGTTGTTGATCCACAGGTTCTAGCTAGAAGACcac                   >  2:3406065/1‑100 (MQ=255)
                                                                         gTAAAGGCATTAGCCCGACGAATACGTCGGGCTAAAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                   >  1:1133690/1‑100 (MQ=255)
                                                                           aaaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACAt                 <  2:1397880/100‑1 (MQ=255)
                                                                           aaaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACAt                 <  2:1324109/100‑1 (MQ=255)
                                                                            aaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATc                >  1:2918364/1‑100 (MQ=255)
                                                                                ttgtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc            >  2:1940158/4‑100 (MQ=39)
                                                                                cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc            <  2:1440686/100‑1 (MQ=255)
                                                                                cATTAGCCCGACGAATACGTCGGGCCACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGACGACCACATCGAca            >  1:2182794/1‑99 (MQ=255)
                                                                                  ttAGCCTGACGAATACGTCGGGCTACAAATATTATTTTGCTGCAGGTGTCTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCtt          <  1:1111459/100‑2 (MQ=255)
                                                                                  ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg          >  1:1223219/1‑100 (MQ=255)
                                                                                  ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg          <  2:3224742/100‑1 (MQ=255)
                                                                                  gtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg          >  2:658547/2‑100 (MQ=255)
                                                                                  gtAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg          >  1:2954057/2‑100 (MQ=255)
                                                                                   tAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa         >  2:1536921/1‑100 (MQ=255)
                                                                                    aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGat        >  2:57239/1‑100 (MQ=255)
                                                                                     gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc       >  2:1618832/1‑100 (MQ=255)
                                                                                     gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc       >  2:466588/1‑100 (MQ=255)
                                                                                     gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc       >  2:2792764/1‑100 (MQ=255)
                                                                                     gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc       <  2:3292859/100‑1 (MQ=255)
                                                                                     gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc       >  2:1744587/1‑100 (MQ=255)
                                                                                        cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatcatc    >  1:1785644/1‑100 (MQ=255)
                                                                                         gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGAtcatca   <  1:632327/100‑1 (MQ=255)
                                                                                          acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaa  <  1:769340/100‑1 (MQ=255)
                                                                                                |                                                                                             
TGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAA  >  NC_000913/3067227‑3067416

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: