Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6340TG32.5% 316.0 / 84.2 160G60G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (11/97);  new base G (52/0);  total (63/97)
Fisher's exact test for biased strand distribution p-value = 1.40e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 231 total aligned reads displayed.

GTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGC  >  NC_000913/433538‑433730
                                                                                                 |                                                                                                    
gTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                                                                                      <  1:10470064/97‑1 (MQ=255)
gTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                                                                                      <  1:11237496/97‑1 (MQ=255)
 tGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                                                                                     <  1:12460225/97‑1 (MQ=255)
   cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                   <  1:10913445/97‑1 (MQ=255)
   cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGCGag                                                                                                   >  1:6768365/1‑96 (MQ=255)
    gTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaa                                                                                                  <  1:1462950/97‑1 (MQ=255)
    gTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaa                                                                                                  <  1:18185599/97‑1 (MQ=255)
     tCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaa                                                                                                 <  1:8721698/97‑1 (MQ=255)
     tCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGgaaa                                                                                                 >  1:3943015/1‑97 (MQ=255)
      cATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaaa                                                                                                <  1:13508262/97‑1 (MQ=255)
      cATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGgaaaa                                                                                                >  1:5099380/1‑97 (MQ=255)
        ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGc                                                                                              <  1:20382937/97‑1 (MQ=255)
         tGTCAAAGATGGCGAAATTGTCGTTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                                                                                             <  1:5540987/97‑1 (MQ=255)
         tGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                                                                                             <  1:3209396/97‑1 (MQ=255)
         tGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGgaaagagc                                                                                             >  1:20503120/1‑93 (MQ=255)
            cAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACCaa                                                                                           >  1:1774678/1‑97 (MQ=255)
            cAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCaaa                                                                                          <  1:3101464/97‑1 (MQ=255)
            cAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCaaa                                                                                          <  1:19926246/97‑1 (MQ=255)
            cAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCaaa                                                                                          >  1:12561418/1‑97 (MQ=255)
              aaGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGACaaaag                                                                                        >  1:4610751/1‑95 (MQ=255)
               aGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                      <  1:10820773/97‑1 (MQ=255)
               aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                      <  1:20157588/97‑1 (MQ=255)
               aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                      <  1:14478800/97‑1 (MQ=255)
               aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                      <  1:6710910/97‑1 (MQ=255)
               aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                      <  1:16167124/97‑1 (MQ=255)
                 aTGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGc                                                                                   <  1:15199199/97‑1 (MQ=255)
                 aTGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGc                                                                                   <  1:15672877/97‑1 (MQ=255)
                 aTGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGGAAAAACAAAggagg                                                                                    >  1:9231178/1‑94 (MQ=255)
                    gCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                               <  1:5925308/97‑1 (MQ=255)
                    gCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                               <  1:6620585/97‑1 (MQ=255)
                    gCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                               <  1:15859749/97‑1 (MQ=255)
                     cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                                                                              <  1:14449524/97‑1 (MQ=255)
                       aaaTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGc                                                                            <  1:10152047/97‑1 (MQ=255)
                       aaaTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGc                                                                            <  1:16319839/97‑1 (MQ=255)
                       aaaTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAAGAGCGAcagg                                                                            >  1:2405018/1‑94 (MQ=255)
                        aaTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGcc                                                                           <  1:16481675/97‑1 (MQ=255)
                         aTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt                                                                          <  1:2755706/97‑1 (MQ=255)
                          ttGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa                                                                         <  1:3701994/97‑1 (MQ=255)
                            gTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGGCAAAAGGGGGACCGCCTCTg                                                                       >  1:4894110/1‑97 (MQ=255)
                             tCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                      <  1:18433038/97‑1 (MQ=255)
                             tCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                      <  1:16618648/97‑1 (MQ=255)
                             tCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                      <  1:18772611/97‑1 (MQ=255)
                             tCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                      <  1:7415750/97‑1 (MQ=255)
                              cGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTc                                                                     <  1:11879246/97‑1 (MQ=255)
                               ggTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca                                                                    <  1:17931026/97‑1 (MQ=255)
                                 tGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTcaca                                                                  <  1:1208388/97‑1 (MQ=255)
                                  gAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCacac                                                                 <  1:13774468/97‑1 (MQ=255)
                                  gAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCacac                                                                 <  1:7517509/97‑1 (MQ=255)
                                    aGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCCCACTc                                                               >  1:340825/1‑97 (MQ=255)
                                    aGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAAGGGGGGCCGCCCCTGTCACACTc                                                               >  1:15033418/1‑97 (MQ=255)
                                     ggTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                              <  1:17683377/97‑1 (MQ=255)
                                     ggTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                              <  1:5223396/97‑1 (MQ=255)
                                     ggTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCCCTGTCACACTcc                                                              >  1:20665981/1‑96 (MQ=255)
                                     ggTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGACCGCCTATGTCACACACg                                                              >  1:6260814/1‑97 (MQ=255)
                                      gTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGACAAAAGGGCGAACGCCCATGTCAAACTCgc                                                             >  1:9205994/1‑96 (MQ=255)
                                         accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACAAAGGGGGGACCGCCTCTGTTACACTCCAAcc                                                         >  1:18453909/1‑97 (MQ=255)
                                         accaccAACGGGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAAcc                                                         <  1:2002570/97‑1 (MQ=255)
                                           ccccaAGGTGGGGGTGAACCACATGCCGAAGTACGCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                       <  1:7828703/95‑1 (MQ=255)
                                           caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                       <  1:19569482/97‑1 (MQ=255)
                                           caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCGCCCCTGTCACACTCGCAAccc                                                        >  1:218782/1‑97 (MQ=255)
                                            accaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAGGAGCGACCGCCCATGTCACACTCGAACCCtc                                                      >  1:17345279/1‑96 (MQ=255)
                                             ccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGACAAAGGGGCGACCGCCCATGTTACACTCGAACCCTGt                                                     >  1:13255026/1‑97 (MQ=255)
                                               aaCGTGGGGGTGAACCCCATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAg                                                   <  1:15433924/97‑1 (MQ=255)
                                               aaCGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAg                                                   <  1:17795241/97‑1 (MQ=255)
                                               aaCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCTATGTCACACCCGAAACCCGTAg                                                   >  1:18576149/1‑97 (MQ=255)
                                                aCGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCTTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGc                                                  <  1:5450724/97‑1 (MQ=255)
                                                aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGCGCCCCCATGTCACACCCCAACCCCGTAgg                                                  >  1:19464067/1‑96 (MQ=255)
                                                aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACAAAAGAGGGAGCCCCCATGTCACACCCCAAACCCGTAGc                                                  >  1:17149704/1‑97 (MQ=255)
                                                 cGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGcc                                                 <  1:11363177/97‑1 (MQ=255)
                                                  gTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCTATGTCCCACTCGCACCCCGTAGCca                                                >  1:10355028/1‑97 (MQ=255)
                                                     cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCcatcc                                             >  1:3536655/1‑96 (MQ=255)
                                                      gggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGACAAAAGTGGGACCGCCTATGTCACACTCGAAACCTGTAGCcatcct                                            >  1:15046673/1‑95 (MQ=255)
                                                      gggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGAGGGACCGCCCATGTCACACCCGAACCCTGTAGGcatcat                                            >  1:4183375/1‑97 (MQ=255)
                                                          gAACCACATGCCGAAGTACACGCGTTGGGTATGGCGGGTAAAAAACCAAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                        <  1:3380307/97‑1 (MQ=255)
                                                          gAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGGAAAACCCAAAGGGGGGGCGCCCCTGTCACAACCGGAGCCCGTAGCCCACATGGtg                                        >  1:9289599/1‑96 (MQ=255)
                                                          gAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGCGACCGCCTATATCACACCCCAACCCCGTAGCCATCATGGTc                                        >  1:9597232/1‑97 (MQ=255)
                                                           acccACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCg                                       <  1:17284826/95‑1 (MQ=255)
                                                           aaCCACATGCCGAAGTCCACGCGTTGCGTATGGGGGGTAAAAAACCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCg                                       <  1:1761672/97‑1 (MQ=255)
                                                            cccAAATGCCGAAGTACCCGGGTTGCTTAGGGCGGGTAAAAAACCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGt                                      <  1:10782627/96‑1 (MQ=255)
                                                              cacaTGCCGAAGTACACGCGTTGCGTATGGGGGGTAAAAAAGCAAAAGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTAc                                    <  1:9017216/97‑1 (MQ=255)
                                                              cacaTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCATGTCACACTCGCACCCCGTAGCCCTCACGGTCGTAc                                    >  1:12711256/1‑97 (MQ=255)
                                                              cacaGGCCGAAGTACCCGCGTTGCTTAGGGCGGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTAc                                    <  1:6830096/97‑1 (MQ=255)
                                                               gattGCCGAAGTCCCCCAGTTGCATATGGGGGGTGAAAAAGAAAAAGTGGGCCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACg                                   <  1:7715471/94‑1 (MQ=255)
                                                               acaTGCCGAAGTCCCCGCGTGGCGTAGGGCGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACg                                   <  1:13194788/97‑1 (MQ=255)
                                                                caTGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGGGACCGCCTATGTTACCCTCGAAACCCGTAGCCATCATTGTCGCACGc                                  >  1:11625594/1‑97 (MQ=255)
                                                                 aTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGGGACCGCCTATGTCACACTCGAAACCTGTAGCCACCATGGTCGTAcccc                                 >  1:19800183/1‑97 (MQ=255)
                                                                  tcccGAGGTCCACGCGTTGCGTAGGGCGGGTAAAAAAGCCAAAGGTCCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                <  1:12578283/95‑1 (MQ=255)
                                                                   gCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCCCACTCGCACCCCGGAGCCATCATGGTCGTACGCCAc                               >  1:4104986/1‑97 (MQ=255)
                                                                     cGAAGTACACGCGTTGCGTAGGGGGGGTAAAAAACCAAAAGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCg                             <  1:2202179/97‑1 (MQ=255)
                                                                               ggtttGCGTAGGGGGGGGGAAAAACCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                   <  1:5724880/94‑1 (MQ=255)
                                                                                 tttGCGTGTGGCGGGTAAAAAACCAAAAGGTGCCCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCa                 <  1:8248560/96‑1 (MQ=255)
                                                                                       ttgggCGGGAAAAAAACCAAAAGGTGCCCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATc           <  1:16887432/94‑1 (MQ=255)
                                                                                        aTGGGGGGTAAAAAAACCAAAGGTGCCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATc           <  1:11297971/97‑1 (MQ=255)
                                                                                           gCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccg       >  1:18054969/1‑97 (MQ=255)
                                                                                            gggggTTAAAAAACCAAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgcc        <  1:10250495/95‑1 (MQ=255)
                                                                                            cGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccgc      >  1:14058522/1‑97 (MQ=255)
                                                                                               ggTTAAAAAACCAAAGGGGGGACCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCt     <  1:8526273/97‑1 (MQ=255)
                                                                                               ggTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTg    >  1:11916751/1‑97 (MQ=255)
                                                                                                 tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGc  >  1:15525056/1‑97 (MQ=255)
                                                                                                 tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGc  >  1:16742498/1‑97 (MQ=255)
                                                                                                 |                                                                                                    
GTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGC  >  NC_000913/433538‑433730

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