Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,300 A→T 15.3% intergenic (‑219/+20) 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,3000AT15.3% 387.3 / 53.1 148intergenic (‑219/+20)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (60/65);  new base T (10/13);  total (70/78)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 152 total aligned reads displayed.

TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGA  >  NC_000913/3790201‑3790399
                                                                                                   |                                                                                                   
tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     >  2:2334619/1‑100 (MQ=255)
    agcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAc                                                                                                 >  1:2311291/1‑100 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                <  1:2486042/100‑1 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                <  2:368362/100‑1 (MQ=255)
        gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                             >  1:690246/1‑100 (MQ=255)
        gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                             <  1:3555598/100‑1 (MQ=255)
         cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                            >  1:2536404/1‑100 (MQ=255)
           atatcAGAGATCAGCACAATTACGATCTCTGTATTACTCTGACAGGATTTTACAGATTTAACAGATTATCCTCATATAAAAGGGCCGGTATTCCAGGCCt                                                                                          <  2:261448/100‑1 (MQ=255)
           atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                          <  1:1337619/100‑1 (MQ=255)
             ataAGAGATCAGCACTATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                         <  2:1522192/99‑1 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                        <  1:3092871/100‑2 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                       <  2:568314/100‑3 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTACCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                       <  1:1284168/99‑1 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                      <  2:1001525/100‑4 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      <  2:2996791/100‑1 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      <  2:314212/100‑1 (MQ=255)
                agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                     >  2:82494/1‑100 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  2:1368782/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  2:1934479/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  1:1786533/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  1:186717/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    >  1:3537802/1‑100 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    <  1:2577345/100‑1 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    <  1:3154221/100‑1 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    >  2:968655/1‑100 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACGCTGACAGAATTTTACATATTTGACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctg                                                                                  <  2:1197713/100‑8 (MQ=37)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  <  2:2055505/100‑1 (MQ=255)
                    aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                 <  1:1806152/100‑1 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAGAAGGGCTGGAATACCAGTCCTTGTTCGTGtt                                                                                <  1:2455333/100‑1 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                <  1:2152968/100‑10 (MQ=38)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                <  2:3504924/100‑1 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                           >  1:416288/1‑100 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGAGTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                            >  1:1630040/1‑99 (MQ=255)
                           caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                          >  1:2773100/1‑100 (MQ=255)
                             aTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAg                                                                        >  2:600811/1‑100 (MQ=255)
                                 cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                    >  1:3231119/1‑100 (MQ=255)
                                   aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                  >  1:3432655/1‑100 (MQ=255)
                                     ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                <  1:587486/100‑1 (MQ=255)
                                     ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                <  2:1621558/100‑1 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                              >  2:3409298/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                              >  1:2081281/1‑100 (MQ=255)
                                        tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                             <  1:426134/100‑1 (MQ=255)
                                           attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACt                                                            <  2:690243/98‑1 (MQ=255)
                                              acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                       <  1:642966/100‑1 (MQ=255)
                                                 cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                    >  1:2333384/1‑100 (MQ=255)
                                                 cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                    <  1:2663259/100‑1 (MQ=255)
                                                  tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   <  2:354123/100‑1 (MQ=255)
                                                  tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   <  2:1032432/100‑1 (MQ=255)
                                                  tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   <  1:19830/100‑1 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                  >  1:3233847/1‑100 (MQ=255)
                                                     cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                <  2:852678/100‑1 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              >  1:3578111/1‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              >  1:490628/1‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAACCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                >  2:2049370/1‑98 (MQ=255)
                                                       ggATTTTAAATAATTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              <  2:2730150/100‑1 (MQ=255)
                                                         aTTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                            <  2:3002478/100‑1 (MQ=255)
                                                           tttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGaa                                          >  2:1560576/1‑100 (MQ=255)
                                                            ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                         <  1:2275628/100‑1 (MQ=255)
                                                             tACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                        <  2:2674439/100‑1 (MQ=255)
                                                             tACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                        >  2:1179913/1‑100 (MQ=255)
                                                             tACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                        <  1:1287003/100‑1 (MQ=255)
                                                               cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                      >  1:3164127/1‑100 (MQ=255)
                                                                  atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCAAACGCATAg                                   >  2:1205362/1‑100 (MQ=255)
                                                                  atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCACgcccacaacg                                   >  2:2872205/1‑91 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                 >  1:2122101/1‑100 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                 >  1:1118691/1‑100 (MQ=255)
                                                                     tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                  <  2:723064/98‑1 (MQ=255)
                                                                     tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                >  2:2621993/1‑100 (MQ=255)
                                                                     tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                <  1:1127528/100‑1 (MQ=255)
                                                                      naCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                               >  2:1883921/2‑100 (MQ=255)
                                                                       aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                              <  1:3033538/100‑1 (MQ=255)
                                                                        cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                             >  1:2954898/1‑100 (MQ=255)
                                                                         aGATTGTCCTCATATAAAAGGGCTGGAATACCTGCCCTTGTTCGTGTTAATCCCAGCTCAGAACAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                              <  2:554280/98‑1 (MQ=255)
                                                                         aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                            <  1:2872228/100‑1 (MQ=255)
                                                                          gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                           <  2:1258844/100‑1 (MQ=255)
                                                                            ttATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGcc                         <  1:3640443/100‑1 (MQ=255)
                                                                              aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTTAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                       <  2:177567/100‑1 (MQ=255)
                                                                                 cTCAGATAAAATTGCTGGAATACCAGCCATTGTTCGTATTAATCCCAGCTCAGAATAACTTTCCCGGACCGGCCAGAACTCATAGCGTTAAAGCCCTTCt                    >  2:1270587/1‑100 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                     <  2:1414024/99‑1 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                    >  2:1169685/1‑100 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   <  1:1351909/100‑1 (MQ=255)
                                                                                  acacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   >  1:3492644/9‑100 (MQ=38)
                                                                                   nacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                  >  1:2255718/8‑100 (MQ=38)
                                                                                    acgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                 >  1:1742281/7‑100 (MQ=38)
                                                                                    aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAAGCCCAGCCCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                 <  2:2542445/100‑1 (MQ=255)
                                                                                     cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                >  2:421489/6‑100 (MQ=39)
                                                                                        acaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat              >  2:3127434/3‑99 (MQ=255)
                                                                                        aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc             <  2:1746784/100‑1 (MQ=255)
                                                                                        aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc             <  2:2649697/100‑1 (MQ=255)
                                                                                          aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  2:516903/1‑100 (MQ=255)
                                                                                           aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc          >  1:294448/1‑100 (MQ=255)
                                                                                           aGGGCTGGAATACCAGTCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc          <  2:2122088/100‑1 (MQ=255)
                                                                                            gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg         <  1:1691753/100‑1 (MQ=255)
                                                                                             ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa        >  2:3419747/1‑100 (MQ=255)
                                                                                              gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt       <  1:1537949/100‑1 (MQ=255)
                                                                                               cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa      <  1:2583955/100‑1 (MQ=255)
                                                                                                tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATag     <  1:2687280/100‑1 (MQ=255)
                                                                                                 gtaaTACCAGCCCTTGTTCGTGTTAATACCCGCTCAGAATAACTTTCCCGGACTGGTACGAACGCATAGCGTCAAAGCCCTTCTAGCAATAATCGATag     >  2:1659899/3‑99 (MQ=255)
                                                                                                  gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag   >  1:1537408/1‑100 (MQ=255)
                                                                                                   aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga  <  1:680917/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGA  >  NC_000913/3790201‑3790399

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