Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,532,150 A→G T35A (ACC→GCC)  b1463 → putative N‑hydroxyarylamine O‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,532,1500AG99.2% 427.0 / ‑4.9 133T35A (ACC→GCC) b1463putative N‑hydroxyarylamine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (58/73);  total (60/73)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 138 total aligned reads displayed.

CTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGT  >  NC_000913/1532060‑1532239
                                                                                          |                                                                                         
cTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGcc                                                                                         <  2:4623168/93‑1 (MQ=255)
cTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGcc                                                                                         >  1:614433/1‑93 (MQ=255)
 tGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCa                                                                                        >  2:3373297/1‑93 (MQ=255)
  gAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCAt                                                                                       <  2:4590576/93‑1 (MQ=255)
    aTCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTc                                                                                     <  2:650449/93‑1 (MQ=255)
     tCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTcc                                                                                    <  1:1257263/93‑1 (MQ=255)
      cACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCg                                                                                   >  1:3132273/1‑93 (MQ=255)
       aCTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGt                                                                                  <  2:3973955/93‑1 (MQ=255)
        cTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGtt                                                                                 <  1:3634155/93‑1 (MQ=255)
         tATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGttt                                                                                >  1:3348373/1‑93 (MQ=255)
         tATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGttt                                                                                >  2:1154562/1‑93 (MQ=255)
             ttGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGaaa                                                                            <  2:1903272/93‑1 (MQ=255)
             ttGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGaaa                                                                            <  1:451109/93‑1 (MQ=255)
             ttGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGaaa                                                                            <  2:1415242/93‑1 (MQ=255)
              tGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGaaaa                                                                           <  2:2732305/93‑1 (MQ=255)
               gCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAAc                                                                          >  1:1375907/1‑93 (MQ=255)
               gCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAAc                                                                          <  2:3701364/93‑1 (MQ=255)
                cccGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAAcc                                                                         <  2:2020970/93‑1 (MQ=255)
                cccGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAAcc                                                                         >  2:2940721/1‑93 (MQ=255)
                 ccGTATTAACTGGTCGTGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCt                                                                        >  2:457791/1‑93 (MQ=255)
                  cGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTc                                                                       <  2:1884214/93‑1 (MQ=255)
                  cGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTc                                                                       >  2:1094342/1‑93 (MQ=255)
                   gTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCg                                                                      <  1:4314534/93‑1 (MQ=255)
                   gTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCg                                                                      >  2:1677750/1‑93 (MQ=255)
                    tATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGa                                                                     >  2:1638900/1‑93 (MQ=255)
                    tATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGa                                                                     >  1:4458253/1‑93 (MQ=255)
                    tATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGa                                                                     <  2:3861563/93‑1 (MQ=255)
                     aTTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGAc                                                                    >  2:2068893/1‑93 (MQ=255)
                         aCTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGttt                                                                >  2:1716201/1‑93 (MQ=255)
                         aCTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGttt                                                                >  2:2914909/1‑93 (MQ=255)
                          cTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGtttt                                                               <  1:3412131/93‑1 (MQ=255)
                          cTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGtttt                                                               <  1:1176085/93‑1 (MQ=255)
                           tGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTtg                                                              >  2:4741189/1‑93 (MQ=255)
                            ggTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTtgc                                                             >  2:1209105/1‑93 (MQ=255)
                                gggAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACTTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGcc                                                         >  1:4221850/1‑93 (MQ=255)
                                gggAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGcc                                                         <  1:3275711/93‑1 (MQ=255)
                                 ggAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCg                                                        >  1:3118360/1‑93 (MQ=255)
                                  gAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGa                                                       <  2:1946688/93‑1 (MQ=255)
                                  gAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGa                                                       <  2:1946689/93‑1 (MQ=255)
                                  gAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGa                                                       >  2:2722665/1‑93 (MQ=255)
                                   aGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAg                                                      <  2:4488255/93‑1 (MQ=255)
                                    gctgctGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAgg                                                     <  2:1512119/93‑1 (MQ=255)
                                    gctgctGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAgg                                                     >  1:1009264/1‑93 (MQ=255)
                                    gctgctGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAgg                                                     >  1:1482769/1‑93 (MQ=255)
                                      tgctgcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGa                                                   >  2:3669615/1‑93 (MQ=255)
                                      tgctgcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGa                                                   >  1:760579/1‑93 (MQ=255)
                                       tctgcGGTCAATATTGATACGATTCGTGCATTGCACCTGAAACACAATTGCGCACTTCCGTTTGAAAACCTCGACGTTTTTCTGCCGCggggg                                                  >  2:2501968/2‑91 (MQ=255)
                                       gctgcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGaa                                                  >  1:2384435/1‑93 (MQ=255)
                                       gctgcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGaa                                                  >  1:2383356/1‑93 (MQ=255)
                                        ctgcGGTCAATATTGATTCGCGTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAATCTCGACGGTTTGCGGCCGAGGGaaa                                                 <  2:1663493/93‑1 (MQ=255)
                                         tgtggTCAATATTGATACGCTTCGTGCATTGCACCTGTAACACAATTGCGCCATTCCGTCTGAAAACCTCGACTTTTTGCTGCCGAGGGatat                                                >  2:4218905/4‑93 (MQ=255)
                                          gcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATa                                               <  2:1830814/93‑1 (MQ=255)
                                          gcGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATa                                               <  1:3928905/93‑1 (MQ=255)
                                            ggTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACa                                             >  1:3133585/1‑93 (MQ=255)
                                             gTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAg                                            <  1:3081751/93‑1 (MQ=255)
                                             gTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCCCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAg                                            >  2:70177/1‑93 (MQ=255)
                                                aaTATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCtt                                         <  2:2903920/93‑1 (MQ=255)
                                                aaTATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCtt                                         <  1:784832/93‑1 (MQ=255)
                                                 atatTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTg                                        <  1:4236514/93‑1 (MQ=255)
                                                    ttGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATa                                     >  2:65128/1‑93 (MQ=255)
                                                    ttGATACGCTACGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAACCCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATa                                     <  2:2447998/93‑1 (MQ=255)
                                                     tGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGCTaa                                    >  1:3891570/1‑93 (MQ=255)
                                                     tGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaa                                    <  1:2102497/93‑1 (MQ=255)
                                                     tGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaa                                    <  2:17189/93‑1 (MQ=255)
                                                     tGATACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaa                                    >  1:3891569/1‑93 (MQ=255)
                                                       aTACGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaatc                                  <  1:1799608/93‑1 (MQ=255)
                                                          cGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaatcaat                               <  2:1130492/93‑1 (MQ=255)
                                                          cGCTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATaatcaat                               <  1:1848562/93‑1 (MQ=255)
                                                            cTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCg                             <  1:4095152/93‑1 (MQ=255)
                                                            cTTCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCg                             <  1:1892717/93‑1 (MQ=255)
                                                              tCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAACCCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGcc                           <  2:907130/93‑1 (MQ=255)
                                                              tCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGcc                           <  1:3541872/93‑1 (MQ=255)
                                                              tCGTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGcc                           <  2:3933149/93‑1 (MQ=255)
                                                                gTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCgg                         >  1:52213/1‑93 (MQ=255)
                                                                gTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCgg                         >  1:3660827/1‑93 (MQ=255)
                                                                gTGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCgg                         >  1:248472/1‑93 (MQ=255)
                                                                 tGCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGa                        >  2:3052654/1‑93 (MQ=255)
                                                                  gCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGaa                       >  2:2376613/1‑93 (MQ=255)
                                                                  gCATTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGaa                       >  1:2584744/1‑93 (MQ=255)
                                                                    aTTGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAaga                     <  1:63584/93‑1 (MQ=255)
                                                                      tGCACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAgaga                   <  2:76149/93‑1 (MQ=255)
                                                                       gccccTGAAACACAACTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGaa                  <  2:4336445/90‑1 (MQ=255)
                                                                        cACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGaaa                 <  1:2371190/93‑1 (MQ=255)
                                                                        cACCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGaaa                 <  1:1451950/93‑1 (MQ=255)
                                                                         cccTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAAc                <  1:1980622/92‑1 (MQ=255)
                                                                         aCCTGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAAc                <  2:2784162/93‑1 (MQ=255)
                                                                          ccTGATACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACt               <  2:1251356/93‑1 (MQ=255)
                                                                            tGAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTgg             <  1:2111357/93‑1 (MQ=255)
                                                                             gAAACACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGt            <  2:4707744/93‑1 (MQ=255)
                                                                              aaaCACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTg           >  1:960288/1‑93 (MQ=255)
                                                                              aaaCACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTg           >  1:2932392/1‑93 (MQ=255)
                                                                              aaaCACAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTg           >  1:1464968/1‑93 (MQ=255)
                                                                               aaCACAATGTCGCCAGTCCGTGTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGa          <  1:1275827/93‑1 (MQ=255)
                                                                                acacAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATCCAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGaa         >  1:1104078/1‑92 (MQ=255)
                                                                                acacAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGAt         <  2:2182887/93‑1 (MQ=255)
                                                                                acacAATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGAt         >  1:1104079/1‑93 (MQ=255)
                                                                                 cacaATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATa        <  1:4024446/93‑1 (MQ=255)
                                                                                 cacaATTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATa        <  1:3033088/93‑1 (MQ=255)
                                                                                    aaTTGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGcc     >  1:2741167/1‑93 (MQ=255)
                                                                                       tGCGCCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCcgt  <  1:2634972/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTGAATCACTATTTTGCCCGTATTAACTGGTCGGGAGCTGCTGCGGTCAATATTGATACGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCCGTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGT  >  NC_000913/1532060‑1532239

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