Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 17.0% intergenic (‑222/+17) waaH ← / ← tdh LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,3030AT17.0% 430.8 / 76.9 170intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (76/65);  new base T (16/13);  total (92/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 175 total aligned reads displayed.

GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402
                                                                                                   |                                                                                                   
  gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTCTTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATACAAGGGCTGGAATAcc                                                                                                   >  2:1597130/1‑100 (MQ=255)
  gcagcTATATAAGAGAGCAGCACAATTAGCAACTCTTTATTACACTGACAGGCAGTTACATATTTAACAGAGTATTCTCAGATCAAAGGGCTGGAATa                                                                                                     <  1:3855626/98‑1 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  2:2892460/1‑100 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  2:1649776/1‑100 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCACAGATAAAAGGGCTGGAATACCAg                                                                                                 >  1:3084191/1‑99 (MQ=255)
      cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                               <  2:1759010/100‑1 (MQ=255)
        atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                             <  2:233652/100‑1 (MQ=255)
        atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                             >  1:3303193/1‑100 (MQ=255)
         tataAGAGATCAGCACAATTCCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCATCAGATAAAAGGACTGGAATACCAGCCCt                                                                                             >  1:2430504/1‑99 (MQ=255)
         tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                            <  2:371454/100‑1 (MQ=255)
         tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCtt                                                                                            >  2:795797/1‑100 (MQ=255)
          ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                           <  1:3391033/100‑2 (MQ=255)
          ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                           <  1:3015104/100‑1 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                          >  2:508323/1‑100 (MQ=255)
           taAAAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttn                                                                                          <  1:2292467/100‑3 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAGCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         >  2:137077/1‑100 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                         <  2:1785552/100‑4 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         >  2:923733/1‑100 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACAGATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         <  1:1846574/100‑1 (MQ=255)
             agagATCAGCATAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                        <  1:813148/100‑1 (MQ=255)
             agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                        <  2:2700715/100‑5 (MQ=39)
              gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         <  1:2310423/98‑1 (MQ=255)
              gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                       <  1:371905/100‑1 (MQ=255)
               agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatct                                                                                      <  2:201653/100‑7 (MQ=38)
               agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatct                                                                                      <  2:1823055/100‑7 (MQ=38)
               agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                      <  1:109605/100‑1 (MQ=255)
               agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                      >  1:1955697/1‑100 (MQ=255)
               agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                      <  2:1298532/100‑1 (MQ=255)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtn                                                                                     <  2:387423/100‑2 (MQ=255)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     >  2:1468564/1‑100 (MQ=255)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  1:463383/100‑1 (MQ=255)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  1:2926130/100‑1 (MQ=255)
                  tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                   <  2:264870/100‑10 (MQ=38)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgagg                                                                                  <  2:3378852/100‑11 (MQ=37)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  >  1:2115816/1‑100 (MQ=255)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  >  2:1273289/1‑100 (MQ=255)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  >  2:1175576/1‑100 (MQ=255)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  >  2:2566364/1‑100 (MQ=255)
                        caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                             <  2:2492102/100‑1 (MQ=255)
                        caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                             <  2:3789447/100‑1 (MQ=255)
                            tACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                         >  1:1585252/1‑100 (MQ=255)
                              cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                       >  2:407543/1‑100 (MQ=255)
                              cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                       >  2:384896/1‑100 (MQ=255)
                                     tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                >  2:646272/1‑100 (MQ=255)
                                      ttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACt                                                               >  2:720/1‑100 (MQ=255)
                                        attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                             >  1:1856969/1‑100 (MQ=255)
                                        attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                             <  2:2934633/100‑1 (MQ=255)
                                        attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                             >  1:479738/1‑100 (MQ=255)
                                        attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                             >  1:3570542/1‑100 (MQ=255)
                                           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGACCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                          >  1:650036/1‑100 (MQ=255)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         >  2:3032556/1‑100 (MQ=255)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         >  2:2057673/1‑100 (MQ=255)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         >  2:816309/1‑100 (MQ=255)
                                               tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                      <  2:3177479/100‑1 (MQ=255)
                                                  cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                   >  2:101900/1‑100 (MQ=255)
                                                  cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                   <  1:1668097/100‑1 (MQ=255)
                                                   aagATTTTACATATTTAACATATTCTGCTCAGTTAAAAGGGCTGTAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCGGAATAACTTTCCCGGACTGGc                                                  <  1:680996/98‑1 (MQ=255)
                                                    ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 >  2:2946532/1‑100 (MQ=255)
                                                    ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 <  1:1067196/100‑1 (MQ=255)
                                                     gaatgtACATATTTAACACATTATCCTCAGATAAAAGGGCTGGAAGACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCTGACTg                                                    <  1:2248310/91‑1 (MQ=255)
                                                        tttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGaa                                             >  1:1366321/1‑100 (MQ=255)
                                                          tACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                           >  1:975009/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGcc                                         >  1:2230840/1‑99 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         <  1:386998/100‑1 (MQ=255)
                                                             atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTAAGAATAACTTTCCCGGACTGGCCAGAACGCAt                                        >  1:1147419/1‑100 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                      >  2:2251721/1‑100 (MQ=255)
                                                                 ttAATAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                    <  2:3265907/100‑1 (MQ=255)
                                                                 ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                    <  2:3543976/100‑1 (MQ=255)
                                                                 ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                    >  2:1521320/1‑100 (MQ=255)
                                                                 ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGATAATCCAAGCTCAGAATAACTTTCACGGACTGGGAGGAACGCATacag                                    >  2:3469522/1‑97 (MQ=255)
                                                                  tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                   >  1:2189552/1‑100 (MQ=255)
                                                                   aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                  <  2:3190267/100‑1 (MQ=255)
                                                                   aaCAGATTATCCTCAGATAAAAGGGCAGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                  >  2:2360425/1‑100 (MQ=255)
                                                                    aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 <  1:1577884/100‑1 (MQ=255)
                                                                    aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 >  1:3000070/1‑100 (MQ=255)
                                                                      aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                               <  2:2281117/100‑1 (MQ=255)
                                                                       gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                              <  2:2541216/100‑1 (MQ=255)
                                                                           aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                          >  2:2834540/1‑100 (MQ=255)
                                                                               acacgaacaaGGGCTGGTATTCCAGCCCTTTTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  2:1810489/9‑100 (MQ=37)
                                                                                cacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                     >  2:2170263/8‑100 (MQ=38)
                                                                                cAGAGAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      <  1:752006/99‑1 (MQ=255)
                                                                                 aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    <  2:1952991/100‑1 (MQ=255)
                                                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCATCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   <  1:2481542/100‑1 (MQ=255)
                                                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   <  1:2485615/100‑1 (MQ=255)
                                                                                   gaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCTGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  >  2:1684643/5‑100 (MQ=38)
                                                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 >  2:3344372/1‑100 (MQ=255)
                                                                                     acaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                >  1:2857103/3‑100 (MQ=255)
                                                                                     aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                <  2:1047267/100‑1 (MQ=255)
                                                                                      caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCGCAGAATAACTTTACCGGACTGGCCCGAACGAATAGCGTCAAAGCCCTGCTGGaaa                  >  1:2485655/2‑97 (MQ=255)
                                                                                       aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat              >  2:3821992/1‑100 (MQ=255)
                                                                                        ngggCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             >  2:460461/2‑100 (MQ=255)
                                                                                        aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAACCCTTCTGGAAatcatc             >  2:3685319/1‑100 (MQ=255)
                                                                                        aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  2:1671313/100‑1 (MQ=255)
                                                                                          ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGCATAACTTTCCCGGACTGGCCCGACCGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa           >  2:1548509/1‑100 (MQ=255)
                                                                                          ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa           >  2:2399905/1‑100 (MQ=255)
                                                                                            cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCATTCTGGAAATCATCGATa         >  2:3668666/1‑100 (MQ=255)
                                                                                               gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACCTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag      <  1:2028200/100‑1 (MQ=255)
                                                                                                aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga     <  1:3175357/100‑1 (MQ=255)
                                                                                                   aCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc  <  2:3823935/100‑1 (MQ=255)
                                                                                                   aCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc  >  1:1719861/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402

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