Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 16.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,3030AT16.0% 323.2 / 48.7 125intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (52/53);  new base T (7/13);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 125 total aligned reads displayed.

CAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAA  >  NC_000913/3790207‑3790401
                                                                                                |                                                                                                  
cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAccn                                                                                                 <  1:1396707/100‑2 (MQ=255)
cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                 <  2:2206899/100‑1 (MQ=255)
 agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                                <  2:1637283/100‑1 (MQ=255)
    tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGccc                                                                                             <  1:1017990/100‑1 (MQ=255)
    tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGccc                                                                                             <  2:2445640/100‑1 (MQ=255)
     atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTTCATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                            <  2:684562/100‑1 (MQ=255)
     atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                            <  1:287257/100‑1 (MQ=255)
       ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                          <  1:306924/100‑2 (MQ=255)
       ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                          <  1:566211/100‑2 (MQ=255)
       ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                          <  2:2653839/100‑2 (MQ=255)
       ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                          >  2:1538064/1‑100 (MQ=255)
        taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                         <  1:2163855/100‑3 (MQ=255)
         aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                        <  2:2794847/100‑1 (MQ=255)
          agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                       <  1:100977/100‑5 (MQ=39)
           gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                      <  2:2743608/100‑6 (MQ=39)
           gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                      <  1:1883292/100‑6 (MQ=39)
             gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                    >  1:2150260/1‑100 (MQ=255)
             gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                    >  1:752958/1‑100 (MQ=255)
              atcatcACAATTACAAACTCTTGATTACACTGACATGATTTTACCTATTTAACAGATTATCCTCAGATAAAAGGGCGGGTATTCCAGCCCttttatct                                                                                     <  1:85085/98‑7 (MQ=25)
              aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTACCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                   <  1:1885501/100‑1 (MQ=255)
              aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                   <  1:1389186/100‑1 (MQ=255)
              aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                   <  1:2573562/100‑1 (MQ=255)
              aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                   <  2:2204714/100‑1 (MQ=255)
               tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                  <  2:902195/100‑1 (MQ=255)
               tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                  <  2:1492978/100‑1 (MQ=255)
               tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                  <  2:1817609/100‑1 (MQ=255)
               tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                  <  1:1401499/100‑1 (MQ=255)
                 aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                >  1:2117827/1‑100 (MQ=255)
                  gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAAt                                                                               >  1:2514511/1‑100 (MQ=255)
                   cacaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                              <  2:517265/100‑1 (MQ=255)
                    acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                             >  1:2398258/1‑100 (MQ=255)
                    acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                              <  1:827240/99‑1 (MQ=255)
                     caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcca                                                                            >  1:1449097/1‑99 (MQ=255)
                      aaTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCa                                                                           <  2:2151597/100‑1 (MQ=255)
                       aTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                            >  1:1176004/1‑98 (MQ=255)
                           cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                      >  1:181830/1‑100 (MQ=255)
                           cGAACTCTTTATTAAAATAAAAGGATTTTAAATATTTAACATATTATCCTCAGATAAAAGAGATAGAATACAAGTCATTATTCGTGTTAATccaaaat                                                                        >  1:769314/1‑93 (MQ=255)
                            gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCan                                                                     <  2:2859965/100‑2 (MQ=255)
                             aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                    >  1:1066910/1‑100 (MQ=255)
                               ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                  <  2:1643324/100‑1 (MQ=255)
                                tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                 <  1:1952440/100‑1 (MQ=255)
                                  tgtattaCACTGACAGGATTTTACATATTTCACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                               <  1:1770103/98‑1 (MQ=255)
                                   ttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACt                                                              <  1:2197149/100‑1 (MQ=255)
                                     attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                            >  2:1627247/1‑100 (MQ=255)
                                        acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                         >  2:2242915/1‑100 (MQ=255)
                                        acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCACTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                         >  1:2790009/1‑100 (MQ=255)
                                          acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                       <  2:1141155/100‑1 (MQ=255)
                                          acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                       >  1:1065402/1‑100 (MQ=255)
                                               cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                  >  2:186745/1‑100 (MQ=255)
                                               cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                  >  2:1654737/1‑100 (MQ=255)
                                               cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                  <  2:2636905/100‑1 (MQ=255)
                                                 ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                >  2:350489/1‑100 (MQ=255)
                                                 ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                <  1:442685/100‑1 (MQ=255)
                                                 ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                <  2:2086558/100‑1 (MQ=255)
                                                 ggATTTTACATATGTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                <  2:2719400/100‑1 (MQ=255)
                                                  gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                               <  2:2035575/100‑1 (MQ=255)
                                                         cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                        <  2:161298/100‑1 (MQ=255)
                                                         cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                        <  2:457480/100‑1 (MQ=255)
                                                            atTTAACAGATTATCCTCAGATAAAATGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                     >  2:1771938/1‑100 (MQ=255)
                                                            atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                     >  2:1594317/1‑100 (MQ=255)
                                                            atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                     >  1:480129/1‑100 (MQ=255)
                                                               tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                  >  2:278026/1‑100 (MQ=255)
                                                               tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                  >  2:1869119/1‑100 (MQ=255)
                                                               tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCATTGTTCGTGTTAATCCCAGCTCAGAATAGCTTTCCCGGACTGGCCCGAACGCATAGCg                                   <  1:1609276/99‑1 (MQ=255)
                                                                   aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCGTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                              >  1:1552443/1‑100 (MQ=255)
                                                                      ttATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGcc                           <  2:374578/100‑1 (MQ=255)
                                                                         tcctccGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATAACAGCTAAGAATAACTTTCCCAAACTGGCCCGAACGCATAGCGTCAACGCCCtt                        >  1:1281041/1‑100 (MQ=255)
                                                                         tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                        >  1:1554041/1‑100 (MQ=255)
                                                                          taacacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                       >  2:2491971/11‑100 (MQ=37)
                                                                          taacacgaacaaGGGCTGGTATTCCAGCCCTTGTGCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                       <  2:2177513/90‑1 (MQ=37)
                                                                             cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                    <  1:157061/100‑1 (MQ=255)
                                                                               gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                  >  2:1477421/1‑100 (MQ=255)
                                                                                aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                 <  2:2257582/100‑1 (MQ=255)
                                                                                aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                 <  1:1744495/100‑1 (MQ=255)
                                                                                aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                 <  2:549406/100‑1 (MQ=255)
                                                                                 tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                >  2:517970/1‑100 (MQ=255)
                                                                                 tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                <  1:417944/100‑1 (MQ=255)
                                                                                  acaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCACGAACGCATAGCGTCAAAGCCCTTCTGGaaaca               >  1:1318248/3‑98 (MQ=255)
                                                                                   aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca              >  1:2570942/1‑100 (MQ=255)
                                                                                   aaaGGGCTGGAATACCAGCCCTGGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAAGCCTAGCGTCAAAGCCCTTCTGGAAatca              >  1:2637401/1‑100 (MQ=255)
                                                                                    aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat             >  1:2211820/1‑100 (MQ=255)
                                                                                    aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat             >  2:408837/1‑100 (MQ=255)
                                                                                    aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat             >  2:1074484/1‑100 (MQ=255)
                                                                                     aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc            >  2:2401139/1‑100 (MQ=255)
                                                                                     aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc            >  2:1550468/1‑100 (MQ=255)
                                                                                     aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc            >  1:2459528/1‑100 (MQ=255)
                                                                                     aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc            <  2:1615693/100‑1 (MQ=255)
                                                                                     aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc            <  1:1920248/100‑1 (MQ=255)
                                                                                      gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca              >  1:1768488/1‑97 (MQ=255)
                                                                                      gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg           <  1:512368/100‑1 (MQ=255)
                                                                                       ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa          >  1:738608/1‑100 (MQ=255)
                                                                                       ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATAATCg           >  2:1686862/1‑99 (MQ=255)
                                                                                        gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa          <  1:227767/99‑1 (MQ=255)
                                                                                        gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt         >  2:902855/1‑100 (MQ=255)
                                                                                        gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt         >  2:1085875/1‑100 (MQ=255)
                                                                                        gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt         >  1:1316869/1‑100 (MQ=255)
                                                                                         cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa        >  2:2864629/1‑100 (MQ=255)
                                                                                          tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATag       <  2:2398269/100‑1 (MQ=255)
                                                                                             aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga    >  2:852364/1‑100 (MQ=255)
                                                                                               tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa  >  2:561680/1‑100 (MQ=255)
                                                                                                |                                                                                                  
CAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAA  >  NC_000913/3790207‑3790401

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: