Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,147,466 G→A 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,4660GA97.3% 349.4 / ‑0.4 112S25S (AGC→AGT) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (52/56);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 114 total aligned reads displayed.

GGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCT  >  NC_000913/2147372‑2147562
                                                                                              |                                                                                                
ggCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCAc                                                                                             >  1:2322294/1‑100 (MQ=255)
 gCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCAcc                                                                                            <  2:3005287/100‑1 (MQ=255)
 gCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCAcc                                                                                            <  2:1346418/100‑1 (MQ=255)
   gggCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGc                                                                                          <  2:1904732/100‑1 (MQ=255)
          atCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcg                                                                                   >  2:980787/1‑100 (MQ=255)
           tCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcgg                                                                                  >  1:781081/1‑100 (MQ=255)
            ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcggc                                                                                 <  2:657302/100‑1 (MQ=255)
            ccGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGggcggc                                                                                 >  1:2116807/1‑100 (MQ=255)
               gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCgag                                                                              <  2:2892523/100‑1 (MQ=255)
               gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCgag                                                                              <  2:1479897/100‑1 (MQ=255)
               gAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCgag                                                                              <  1:95950/100‑1 (MQ=255)
               cttccgatctCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCgag                                                                              <  2:2967327/91‑1 (MQ=255)
                aaTAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCgaga                                                                             <  1:92831/100‑1 (MQ=255)
                 aTAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGaa                                                                            >  1:489180/1‑100 (MQ=255)
                  tAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGaaa                                                                           <  1:2386797/100‑1 (MQ=255)
                    gCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATc                                                                         <  2:728623/100‑1 (MQ=255)
                     cAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATcc                                                                        <  1:2622772/100‑1 (MQ=255)
                      attCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATccc                                                                       >  1:2038604/3‑100 (MQ=255)
                      aGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATccc                                                                       <  1:2882159/100‑1 (MQ=255)
                      aGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATccc                                                                       <  2:2209055/100‑1 (MQ=255)
                       gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCa                                                                      >  1:1079894/1‑100 (MQ=255)
                       gTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCa                                                                      >  2:482478/1‑100 (MQ=255)
                         caccacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAGATcggaag                                                                    >  1:2967353/1‑95 (MQ=255)
                         caccacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAAc                                                                    >  1:2098318/1‑100 (MQ=255)
                         caccacGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAAc                                                                    >  2:1432845/1‑100 (MQ=255)
                              cgccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCa                                                               <  2:2634008/100‑1 (MQ=255)
                               gccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCAc                                                              >  1:2900102/1‑100 (MQ=255)
                                ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTCCTCACCGCCCGGGCGGCGAGAAATCCCAACCTCCAc                                                              >  1:557020/1‑99 (MQ=255)
                                ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg                                                             <  1:1146358/100‑1 (MQ=255)
                                ccgcGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg                                                             <  1:435721/100‑1 (MQ=255)
                                  gcgATATTCGCCCACCGCCAGACTACGGGCATAGTAACTGTCCGCGACCGTTTCCACACTACTCACCGCCCGGGCGGCGAGAAATCCCCGCATCCACGAc                                                           >  2:1997266/1‑100 (MQ=255)
                                  gcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                           >  1:2744726/1‑100 (MQ=255)
                                  gcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                           >  1:402404/1‑100 (MQ=255)
                                  gcgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc                                                           <  1:1193502/100‑1 (MQ=255)
                                   cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                          <  1:1531491/100‑1 (MQ=255)
                                   cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                          <  1:96241/100‑1 (MQ=255)
                                   cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc                                                          <  1:3029978/100‑1 (MQ=255)
                                   cgATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCAAGAAATCCCAACATCCACGAcc                                                          >  2:1017490/1‑100 (MQ=255)
                                      tatTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGt                                                       <  1:2061257/100‑1 (MQ=255)
                                       atTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTc                                                      <  1:1744265/100‑1 (MQ=255)
                                       atTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTc                                                      <  2:1507079/100‑1 (MQ=255)
                                         tCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGATAAATCCCAACATCCACGACCAGTCAt                                                    <  1:1851597/100‑1 (MQ=255)
                                           gCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAc                                                  <  1:1705880/100‑1 (MQ=255)
                                             ccACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGATAAATCCCAACATCCACGACCAGTCATAcgg                                                <  2:1042325/100‑1 (MQ=255)
                                             ccACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCAAGGGCGGCGCAAAATCCCAACATCCACGACCAGTCATAcg                                                 >  1:2531245/1‑99 (MQ=255)
                                              cACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggc                                               <  1:1774426/100‑1 (MQ=255)
                                                ccGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATAcggcgg                                             <  2:1242451/100‑1 (MQ=255)
                                                 cGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACggcggc                                            >  2:2414347/1‑100 (MQ=255)
                                                   ccAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                          >  2:761279/1‑100 (MQ=255)
                                                   ccAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTg                                          <  2:1022595/100‑1 (MQ=255)
                                                    cAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGc                                         >  2:2760692/1‑100 (MQ=255)
                                                     aGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGcc                                        <  1:552052/100‑1 (MQ=255)
                                                     aGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGcc                                        <  2:1401966/100‑1 (MQ=255)
                                                      gACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCa                                       >  1:1089830/1‑100 (MQ=255)
                                                           cGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCa                                  >  2:420619/1‑100 (MQ=255)
                                                           cGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCa                                  >  2:1461211/1‑100 (MQ=255)
                                                           cGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCa                                  >  1:2785064/1‑100 (MQ=255)
                                                            gggCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAg                                 <  1:2218706/100‑1 (MQ=255)
                                                             ggCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAgg                                >  2:918515/1‑100 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCCCACTACTCACCGCACGGTCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  2:2716722/100‑1 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  2:85700/100‑1 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  2:1646168/100‑1 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  1:2813627/100‑1 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  1:2013490/100‑1 (MQ=255)
                                                              gCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAggg                               <  2:312411/100‑1 (MQ=255)
                                                                 taataaCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtat                            <  1:725454/100‑1 (MQ=255)
                                                                 taataaCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtat                            >  1:2565952/1‑100 (MQ=255)
                                                                  aataaCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGtata                           <  2:3033069/100‑1 (MQ=255)
                                                                   ataaCTGTCCGCGACCGTTTCCACACTACTCACCGCATGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATAc                          <  1:771611/100‑1 (MQ=255)
                                                                    taaCTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACa                         >  1:1844330/1‑100 (MQ=255)
                                                                      aCTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATc                       <  1:2692196/100‑1 (MQ=255)
                                                                       cTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCg                      >  2:1944954/1‑100 (MQ=255)
                                                                         gTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGc                     >  2:645311/1‑99 (MQ=255)
                                                                         gTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                    >  1:1281667/1‑100 (MQ=255)
                                                                         gTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCa                    <  1:2414369/100‑1 (MQ=255)
                                                                          tCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAt                   >  1:1717948/1‑100 (MQ=255)
                                                                           ccGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtc                  >  1:3015018/1‑100 (MQ=255)
                                                                           ccGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtc                  >  2:1148881/1‑100 (MQ=255)
                                                                            cgcgACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCCTACGGCGGCTGCCAGTTCAGGGTATACATCGCAt                   <  2:27658/98‑1 (MQ=255)
                                                                            cgcgACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtct                 >  2:1034958/1‑100 (MQ=255)
                                                                             gcgACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCAtctc                <  2:736018/100‑1 (MQ=255)
                                                                                  cGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt           <  2:489175/100‑1 (MQ=255)
                                                                                  cGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAt           >  2:1037458/1‑100 (MQ=255)
                                                                                   gTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt          <  1:1019451/100‑1 (MQ=255)
                                                                                   gTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTAtt          <  1:833492/100‑1 (MQ=255)
                                                                                     ttCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa        <  1:761286/100‑1 (MQ=255)
                                                                                     ttCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa        <  2:1531052/100‑1 (MQ=255)
                                                                                     ntCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCa        >  2:2280677/2‑100 (MQ=255)
                                                                                      tCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt       >  2:1984077/1‑100 (MQ=255)
                                                                                      tCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt       >  1:2066397/1‑100 (MQ=255)
                                                                                      tCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCAt       <  2:2821702/100‑1 (MQ=255)
                                                                                       ccACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATc      >  2:1291860/1‑100 (MQ=255)
                                                                                        cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc     >  2:888721/1‑100 (MQ=255)
                                                                                        cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc     >  1:703209/1‑100 (MQ=255)
                                                                                        cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc     <  1:352736/100‑1 (MQ=255)
                                                                                        cacaCTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATcc     <  1:1432857/100‑1 (MQ=255)
                                                                                          cacTACTCACTGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGc   <  2:1996006/100‑1 (MQ=255)
                                                                                          cacTACTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGc   >  1:259376/1‑100 (MQ=255)
                                                                                           actactCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCGGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCt  <  2:1628668/100‑1 (MQ=255)
                                                                                           actactCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCt  >  1:117696/1‑100 (MQ=255)
                                                                                              |                                                                                                
GGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCATACGGCGGCTGCCAGTTCAGGGTATACATCGCATCTCCTTATTCATCCGCT  >  NC_000913/2147372‑2147562

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