Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,714,643 T→G I89M (ATT→ATG gstA → glutathionine S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,714,6430TG95.5% 588.8 / 3.4 177I89M (ATT→ATGgstAglutathionine S‑transferase
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (72/97);  total (74/103)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

TGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTACA  >  NC_000913/1714546‑1714740
                                                                                                 |                                                                                                 
tgGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTc                                                                                                 <  1:1077127/100‑1 (MQ=255)
 ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTcc                                                                                                >  1:1201494/1‑100 (MQ=255)
 ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTcc                                                                                                <  1:5223342/100‑1 (MQ=255)
  gATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTccc                                                                                               <  2:4573172/100‑1 (MQ=255)
    tGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGc                                                                                             <  2:2381176/100‑1 (MQ=255)
     gACGGTACTGTGCGGACGGAAGGCGTAGCGATGATGCAGTATCTGGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCt                                                                                            <  2:4148467/100‑1 (MQ=255)
      aCGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCta                                                                                           >  2:3714201/1‑100 (MQ=255)
       cGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCtat                                                                                          <  2:1164139/100‑1 (MQ=255)
       cGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCtat                                                                                          >  2:2070975/1‑100 (MQ=255)
        ggTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCtata                                                                                         <  2:2847012/100‑1 (MQ=255)
        ggTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCtata                                                                                         <  1:3738782/100‑1 (MQ=255)
         gTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATaa                                                                                        <  2:2749369/100‑1 (MQ=255)
         gTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGGATGTCCCGCTATaa                                                                                        <  2:4017066/100‑1 (MQ=255)
           acgttGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGGGTCCCCGACCGCCAGTTGCTGGCACCGGTAAGCAGTATGTCCCGCTATaaaa                                                                                      <  2:961507/97‑1 (MQ=255)
           aCTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATaaaa                                                                                      <  2:5021134/100‑1 (MQ=255)
            cTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCTACAGCGTCCCCGACCGCGAGTTGCTGGCACCGGTAAACAGTATGTCCAGCTATAAAAc                                                                                     <  2:2694694/100‑1 (MQ=255)
              ttGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCa                                                                                   <  1:3662310/100‑1 (MQ=255)
                gCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATc                                                                                 <  2:508785/100‑1 (MQ=255)
                 cTGACGGAAGGCGTAGCGATTATGCTGTATCTTGCCTACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCg                                                                                <  2:5214164/100‑1 (MQ=255)
                 cTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCg                                                                                >  1:3226061/1‑100 (MQ=255)
                 cTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCg                                                                                >  2:3125273/1‑100 (MQ=255)
                     cGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATg                                                                            >  1:2132201/1‑100 (MQ=255)
                     cGGAAGACGTAGCGATTATGCAGTATCTTACCGACAGCGTCACCGACCGCCAGTTACGAGCACCGGTAACAAGTATGTCCCGGTGTAAAACCGTCGCAtc                                                                            >  2:5240320/1‑99 (MQ=255)
                      ggAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATgg                                                                           <  2:1504099/100‑1 (MQ=255)
                       gAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATTTCCCGCTATAAAACCATCGAATGGc                                                                          <  1:3966690/100‑1 (MQ=255)
                          ggCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGa                                                                       <  2:4379998/100‑1 (MQ=255)
                          ggCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGa                                                                       >  2:2785946/1‑100 (MQ=255)
                            cGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAAt                                                                     <  2:1023727/100‑1 (MQ=255)
                             gTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAAtt                                                                    >  2:263635/1‑100 (MQ=255)
                             gTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAAtt                                                                    >  1:940315/1‑100 (MQ=255)
                             gTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCCATAAAACCATCGAATGGCTGAAtt                                                                    >  2:3693878/1‑100 (MQ=255)
                               aGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTAc                                                                  <  1:124033/100‑1 (MQ=255)
                                gCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACa                                                                 <  1:1061767/100‑1 (MQ=255)
                                gCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACa                                                                 <  2:4071943/100‑1 (MQ=255)
                                 cGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACAt                                                                <  1:1224510/100‑1 (MQ=255)
                                 cGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACAt                                                                >  1:2144951/1‑100 (MQ=255)
                                 cGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACAt                                                                <  1:2513117/100‑1 (MQ=255)
                                    ttATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGc                                                             <  2:2687728/100‑1 (MQ=255)
                                    ttATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGc                                                             >  1:2717073/1‑100 (MQ=255)
                                       tGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCAc                                                          >  2:3575108/1‑100 (MQ=255)
                                        gCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCAcc                                                         >  1:2878666/1‑100 (MQ=255)
                                        gCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCAcc                                                         >  2:604626/1‑100 (MQ=255)
                                        gCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCAcc                                                         >  1:1791313/1‑100 (MQ=255)
                                              tCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTCTGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTg                                                   >  1:1424171/1‑100 (MQ=255)
                                              tCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTg                                                   <  2:1377469/100‑1 (MQ=255)
                                                ttGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAACACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCa                                                 <  2:551351/100‑1 (MQ=255)
                                                ttGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCa                                                 <  2:215493/100‑1 (MQ=255)
                                                   ccGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATaa                                              >  1:3793806/1‑100 (MQ=255)
                                                     gACAGCGTCCCCGACCGCTTGTTGCTGGCACTGGTTTACAGTGTGTTCCGCTATAAATCTATCGAATGGCTGATTTACATCGCCACCGCTCTGCATACAg                                            >  2:2635228/1‑100 (MQ=255)
                                                       cAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGt                                          <  1:5341846/100‑1 (MQ=255)
                                                       cAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGt                                          >  1:2266917/1‑100 (MQ=255)
                                                       cAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGt                                          >  1:3669107/1‑100 (MQ=255)
                                                        aGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGtt                                         <  2:4370073/100‑1 (MQ=255)
                                                         gCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGttt                                        <  1:975309/100‑1 (MQ=255)
                                                          cctCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTc                                       <  2:830627/98‑1 (MQ=255)
                                                          cctCCACGACCGCTAGTTGCTGGAACCGGTAAAGAGTATGTCCCGCTATAAAACCATCGATTGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTc                                       <  2:4362701/98‑1 (MQ=255)
                                                          cGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTc                                       <  1:3021628/100‑1 (MQ=255)
                                                            tCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTcac                                     <  2:2571006/100‑1 (MQ=255)
                                                            tCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTcac                                     <  1:63767/100‑1 (MQ=255)
                                                            tCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTcac                                     <  1:1077988/100‑1 (MQ=255)
                                                             ccccGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTcaca                                    <  1:1706344/100‑1 (MQ=255)
                                                              cccGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGTCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCacac                                   >  1:4922036/1‑100 (MQ=255)
                                                              cccGACCGCAAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCacac                                   >  1:2412524/1‑100 (MQ=255)
                                                               ccGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACAcc                                  <  1:1518443/100‑1 (MQ=255)
                                                               ccGACCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAAGGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACAcc                                  <  1:1189973/100‑1 (MQ=255)
                                                                  aCCGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACctct                               >  2:3326636/1‑100 (MQ=255)
                                                                   ccGCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTg                              >  2:2479740/1‑100 (MQ=255)
                                                                     gCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGtt                            >  1:3016623/1‑100 (MQ=255)
                                                                     gCCAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGtt                            >  2:1286612/1‑100 (MQ=255)
                                                                      ccAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGttt                           <  2:3635305/100‑1 (MQ=255)
                                                                       cAGTTGCTGGCACCGGTAAACAGTATGTCTCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATACAGGTTTCACACCTCTGTTTc                          <  2:1537990/100‑1 (MQ=255)
                                                                       cAGTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTc                          <  1:4175593/100‑1 (MQ=255)
                                                                         gTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGc                        <  2:3195993/100‑1 (MQ=255)
                                                                         gTTGCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGc                        <  2:2292391/100‑1 (MQ=255)
                                                                           tGCTGGCACCGGTAAACAGTATTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGccc                      <  1:103256/100‑1 (MQ=255)
                                                                            gCTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCt                     <  2:471915/100‑1 (MQ=255)
                                                                             cTGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTg                    <  2:4904439/100‑1 (MQ=255)
                                                                              tGGCACCGGTAAACAGTATTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGa                   <  1:4486691/100‑1 (MQ=255)
                                                                              tGGCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGa                   >  2:3477046/1‑100 (MQ=255)
                                                                               ggCACCGGTAAACAGTATGTCCCGCTATGAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGAt                  >  1:4387029/1‑100 (MQ=255)
                                                                               ggCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGCGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGAt                  >  2:2206430/1‑100 (MQ=255)
                                                                               ggCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGAt                  >  2:3920329/1‑100 (MQ=255)
                                                                               ggCACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGAt                  <  1:4143410/100‑1 (MQ=255)
                                                                                 naCCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATac                >  1:5486300/2‑100 (MQ=255)
                                                                                 cACCGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATac                <  2:581591/100‑1 (MQ=255)
                                                                                  aCCGGTAAACAGTATTTCCTGCTATAAAACCATCGTATGGCTGAATTACATCGCCACCTAGCTGCCTAAAGGTTTCACAACTCTGTTTGGCCCTGATaca               >  2:2828290/1‑100 (MQ=255)
                                                                                    cGGTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACAcc             <  2:1185431/100‑1 (MQ=255)
                                                                                     ggTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCg            >  2:4028024/1‑100 (MQ=255)
                                                                                     ggTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCg            >  1:2138479/1‑100 (MQ=255)
                                                                                      gTAAACAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCgg           >  1:4549718/1‑100 (MQ=255)
                                                                                        aaaCAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGaa         <  2:1405426/100‑1 (MQ=255)
                                                                                        aaaCAGTAGGTCCCGCTATAAAACCATCGAATGGCGGAATTACATCGCCACCGAGCGGCATAAAGGTTTCACACGTCTCTTTCGCCCTGATACAACGGaa         >  2:3026373/1‑100 (MQ=255)
                                                                                         aaCAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAag        <  2:647932/100‑1 (MQ=255)
                                                                                           cAGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAagag      <  1:3386214/100‑1 (MQ=255)
                                                                                            aGTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAagagc     >  1:2286940/1‑99 (MQ=255)
                                                                                             gTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTa    <  2:4862791/100‑1 (MQ=255)
                                                                                             gTATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTa    >  1:5459265/1‑100 (MQ=255)
                                                                                              tATTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTAc   <  2:5321746/100‑1 (MQ=255)
                                                                                              tATGTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTAc   <  1:427844/100‑1 (MQ=255)
                                                                                               aTTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGGTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGGACa  >  2:1913823/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
TGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCCAGTTGCTGGCACCGGTAAACAGTATTTCCCGCTATAAAACCATCGAATGGCTGAATTACATCGCCACCGAGCTGCATAAAGGTTTCACACCTCTGTTTCGCCCTGATACACCGGAAGAGTACA  >  NC_000913/1714546‑1714740

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