Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,466 G→A 21.1% S25S (AGC→AGT)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,4660GA21.1% 254.3 / 56.3 118S25S (AGC→AGT) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base G (45/48);  new base A (16/9);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

AGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTC  >  NC_000913/2147367‑2147565
                                                                                                   |                                                                                                   
agTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTg                                                                                                     <  1:744650/100‑1 (MQ=255)
   aTGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTc                                                                                                  >  1:2908912/1‑100 (MQ=255)
   aTGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTc                                                                                                  <  1:645337/100‑1 (MQ=255)
    tGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCa                                                                                                 >  2:291821/1‑100 (MQ=255)
      gCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCAcc                                                                                               >  1:1840327/1‑100 (MQ=255)
              tatCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggc                                                                                       >  1:57163/1‑100 (MQ=255)
              tatCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggc                                                                                       >  2:2492236/1‑100 (MQ=255)
              tatCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggc                                                                                       >  1:2611856/1‑100 (MQ=255)
                tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcgg                                                                                     >  2:947175/1‑100 (MQ=255)
                tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcgg                                                                                     >  1:424378/1‑100 (MQ=255)
                tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcgg                                                                                     <  1:3213056/100‑1 (MQ=255)
                tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcgg                                                                                     >  2:1404074/1‑100 (MQ=255)
                tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcgg                                                                                     >  2:757857/1‑100 (MQ=255)
                 ccGGAATAGCAGGCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcggc                                                                                    >  2:2278264/1‑100 (MQ=255)
                  cGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGgcggcg                                                                                   <  1:2910360/100‑1 (MQ=255)
                  cGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGgcggcg                                                                                   <  1:1423032/100‑1 (MQ=255)
                   ggAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGAGACCGTTTCCACACTGCTCACCGCACGGGCGGCga                                                                                  >  2:924861/1‑100 (MQ=255)
                    gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgag                                                                                 <  2:3208802/100‑1 (MQ=255)
                    gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgag                                                                                 >  2:1684404/1‑100 (MQ=255)
                     aaTAGCGGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAAGAACTGACCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgaga                                                                                <  1:651576/100‑1 (MQ=255)
                     aaTAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgaga                                                                                <  2:3182899/100‑1 (MQ=255)
                        aGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAAt                                                                             >  2:3335281/1‑100 (MQ=255)
                        aGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAAt                                                                             >  2:1032450/1‑100 (MQ=255)
                        aGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAAt                                                                             <  2:1845746/100‑1 (MQ=255)
                         gtaGTCACCACGCGGCGATATTCGCCCACCGACAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATc                                                                            >  2:2334599/3‑100 (MQ=255)
                          cAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATcc                                                                           <  2:2882326/100‑1 (MQ=255)
                            gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                         <  2:3100434/100‑1 (MQ=255)
                            gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                         <  2:779877/100‑1 (MQ=255)
                            gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCa                                                                         >  2:1291247/1‑100 (MQ=255)
                             tCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATcc                                                                           <  2:26250/97‑1 (MQ=255)
                             tCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCaa                                                                        <  1:2639580/100‑1 (MQ=255)
                              caccacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAc                                                                       >  2:2235580/1‑100 (MQ=255)
                                ccacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACAt                                                                     >  1:1237942/1‑100 (MQ=255)
                                 cacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATc                                                                    >  2:3228716/1‑100 (MQ=255)
                                   cgccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCa                                                                  >  2:1053220/1‑100 (MQ=255)
                                   cgccgcGATATTCGCCCACCGCCAGACCGCGGGCATAATAACTGGCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATcc                                                                   <  1:2649763/99‑1 (MQ=255)
                                     ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCAc                                                                 >  1:1116927/1‑99 (MQ=255)
                                     ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg                                                                >  2:1176992/1‑100 (MQ=255)
                                     ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGCTTCCACACTGCTCACCGCCCGGGCGGCGAGAAATCCCACCATCCACg                                                                >  1:1283763/1‑100 (MQ=255)
                                       gcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                              >  1:1809148/1‑100 (MQ=255)
                                           tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGt                                                          <  1:2699317/100‑1 (MQ=255)
                                           tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGt                                                          <  2:645144/100‑1 (MQ=255)
                                             ttCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCa                                                        >  2:2705182/1‑100 (MQ=255)
                                               cGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCCCGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATa                                                      >  2:965403/1‑100 (MQ=255)
                                                gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAc                                                     <  1:86423/100‑1 (MQ=255)
                                                   cACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggc                                                  >  1:2923487/1‑100 (MQ=255)
                                                    aCCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggc                                                  <  1:2930561/99‑1 (MQ=255)
                                                      cGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACggcggc                                               <  1:1176981/100‑1 (MQ=255)
                                                        ccAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                             <  1:1543852/100‑1 (MQ=255)
                                                         cAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGc                                            <  1:1045326/100‑1 (MQ=255)
                                                         cAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGACAAATCCCAACATCCACGACCAGTCATACGGCGGCTGc                                            <  1:1854771/100‑1 (MQ=255)
                                                          aGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGcc                                           <  1:1236399/100‑1 (MQ=255)
                                                            aCTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAg                                         >  1:1152030/1‑100 (MQ=255)
                                                              tACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGtt                                       >  2:956077/1‑100 (MQ=255)
                                                              tACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGtt                                       >  2:1987320/1‑100 (MQ=255)
                                                               aCGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTc                                      <  2:1844828/100‑1 (MQ=255)
                                                               aCGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGGTc                                      >  2:2758515/1‑100 (MQ=255)
                                                               aCGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTc                                      <  1:757851/100‑1 (MQ=255)
                                                                  ggCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAgg                                   <  1:2212736/100‑1 (MQ=255)
                                                                  ggCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAgg                                   <  2:18999/100‑1 (MQ=255)
                                                                   gCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                  >  1:1629726/1‑100 (MQ=255)
                                                                   gCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                  <  1:1973014/100‑1 (MQ=255)
                                                                   gCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                  <  1:1797168/100‑1 (MQ=255)
                                                                   gCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                                  <  2:1836075/100‑1 (MQ=255)
                                                                    cATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGt                                 >  2:2556300/1‑100 (MQ=255)
                                                                      taataaCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtat                               <  2:354994/100‑1 (MQ=255)
                                                                      taataaCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtat                               <  2:2746978/100‑1 (MQ=255)
                                                                        ataaCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATAc                             >  2:1988100/1‑100 (MQ=255)
                                                                          aaCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACAt                           <  2:451345/100‑1 (MQ=255)
                                                                           aCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                          >  2:2433802/1‑100 (MQ=255)
                                                                           aCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                          >  1:734958/1‑100 (MQ=255)
                                                                           aCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                          <  2:1237953/100‑1 (MQ=255)
                                                                            cTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCg                         >  1:612498/1‑100 (MQ=255)
                                                                            cTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCg                         >  2:1380013/1‑100 (MQ=255)
                                                                            cTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCg                         >  2:2666942/1‑100 (MQ=255)
                                                                              gTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                       <  1:1032441/100‑1 (MQ=255)
                                                                              gTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                       >  1:1740773/1‑100 (MQ=255)
                                                                              gTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                       <  1:1421353/100‑1 (MQ=255)
                                                                                ccGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtc                     >  1:838204/1‑100 (MQ=255)
                                                                                   cgACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTcc                  <  2:1261112/100‑1 (MQ=255)
                                                                                    gACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCt                 >  2:3394025/1‑100 (MQ=255)
                                                                                    gACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCt                 <  2:2611875/100‑1 (MQ=255)
                                                                                    gACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCt                 <  2:2905659/100‑1 (MQ=255)
                                                                                      ccGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTa               >  2:320703/1‑100 (MQ=255)
                                                                                       cGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              >  2:2689858/1‑100 (MQ=255)
                                                                                       cGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              <  2:879172/100‑1 (MQ=255)
                                                                                       cGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              <  2:2157228/100‑1 (MQ=255)
                                                                                        gTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt             >  1:322385/1‑100 (MQ=255)
                                                                                        gTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt             >  1:896003/1‑100 (MQ=255)
                                                                                        gTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt              >  1:61798/1‑99 (MQ=255)
                                                                                        gTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt             <  1:898570/100‑1 (MQ=255)
                                                                                         tttCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTc            >  1:2557156/1‑100 (MQ=255)
                                                                                          ttCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa           <  1:2620524/100‑1 (MQ=255)
                                                                                          ttCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa           <  2:945248/100‑1 (MQ=255)
                                                                                          ttCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa           <  2:424382/100‑1 (MQ=255)
                                                                                           tCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt          >  2:2469907/1‑100 (MQ=255)
                                                                                            ccACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCCTATTCATc         >  1:302713/1‑100 (MQ=255)
                                                                                             cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCATTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc        <  1:358502/100‑1 (MQ=255)
                                                                                             cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc        >  1:444270/1‑100 (MQ=255)
                                                                                               cacTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGc      <  2:612502/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCTTTC  >  NC_000913/2147367‑2147565

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