Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 19.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,3030AT19.0% 459.2 / 93.1 194intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (79/78);  new base T (17/20);  total (96/98)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790206‑3790402
                                                                                                 |                                                                                                   
gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   <  2:2699023/100‑1 (MQ=255)
gcacgtataaaAGCGATCAGAAGAATTACGATCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAGAGGGCTGGAATAc                                                                                                    <  1:2599990/94‑1 (MQ=255)
 cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                  >  2:1467025/1‑100 (MQ=255)
  agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                                 <  2:3054266/100‑1 (MQ=255)
   gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  1:63744/1‑100 (MQ=255)
   gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  1:3517983/1‑100 (MQ=255)
   gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  2:3564762/1‑100 (MQ=255)
    cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGcc                                                                                               <  1:514591/100‑1 (MQ=255)
    cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                               <  1:1888806/100‑1 (MQ=255)
      atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                             <  1:3201774/100‑1 (MQ=255)
       tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                            <  2:321062/100‑1 (MQ=255)
        ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                           <  2:514197/100‑1 (MQ=255)
        ataAGAGAGCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATGTTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                           <  1:294044/100‑2 (MQ=255)
         taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                          >  2:1968278/1‑100 (MQ=255)
          aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                         <  2:1188888/100‑4 (MQ=255)
          aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         <  1:2983931/100‑1 (MQ=255)
           agagATCAGCACAATTACGAACTTTTTATTACACTGACAGTATTTTACATATTTAACAGATTATCCTCATATATTTTGGCTAGTATTCCAGCCCttt                                                                                           <  1:1006462/97‑2 (MQ=255)
            gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                       <  2:2403023/100‑6 (MQ=39)
            gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                       <  2:1933987/100‑6 (MQ=39)
              gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctg                                                                                     <  1:1815282/100‑8 (MQ=38)
              gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     <  2:3741347/100‑1 (MQ=255)
                tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatccgag                                                                                   <  2:3632818/100‑10 (MQ=38)
                 cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  <  2:2743647/100‑1 (MQ=255)
                  aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                 >  2:2225130/1‑100 (MQ=255)
                    cacaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                               >  2:2870400/1‑100 (MQ=255)
                         ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGc                                                                          <  2:1581067/100‑1 (MQ=255)
                          tACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                         <  2:3587940/100‑1 (MQ=255)
                           aCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTc                                                                        <  2:981787/100‑1 (MQ=255)
                           aCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTc                                                                        <  1:3563880/100‑1 (MQ=255)
                            cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                       >  2:647705/1‑100 (MQ=255)
                                ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCCCAGAAt                                                                   >  2:2031036/1‑100 (MQ=255)
                                ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATACCAGCTCAGAAt                                                                   >  1:838977/1‑100 (MQ=255)
                                 tctTTATTACACTGACAGGATTTTACGTATTTAACAGATTATCCTCGGATAAAAGGGCTGAAATACCAGCCCTTGTTCGTGTAAAACGCAGCTCAg                                                                      >  1:2876585/1‑96 (MQ=255)
                                 tctTTATTACACTGACAGGATTTTACATATTTAACAGTTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                  <  1:2574970/100‑1 (MQ=255)
                                 tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                  <  2:3144129/100‑1 (MQ=255)
                                  ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCTTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                 >  1:766372/1‑100 (MQ=255)
                                  ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                 >  2:2331548/1‑100 (MQ=255)
                                  ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                 >  2:3597081/1‑100 (MQ=255)
                                       ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTc                                                            >  1:1498916/1‑100 (MQ=255)
                                         acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                          <  1:3708267/100‑1 (MQ=255)
                                          cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTTAGAATAACTTTCCCg                                                         >  2:3387472/1‑100 (MQ=255)
                                          cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         >  2:862659/1‑100 (MQ=255)
                                           acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                        <  2:3128827/100‑1 (MQ=255)
                                            cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                       >  2:2472897/1‑100 (MQ=255)
                                             tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                      <  2:1153021/100‑1 (MQ=255)
                                             tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                      >  1:3787692/1‑100 (MQ=255)
                                              gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                     >  2:508839/1‑100 (MQ=255)
                                               aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTg                                                    <  2:1531273/100‑1 (MQ=255)
                                                ctggATTTTACATATTTAACAGATTATGCTCAGATAAAAGGGCTGGAATTCCAGCCCTTGTTCGTGTTAATCACATCTAAGAATAACttt                                                             >  1:3444616/3‑90 (MQ=255)
                                                cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                   >  1:1841959/1‑100 (MQ=255)
                                                  ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 <  2:695268/100‑1 (MQ=255)
                                                  ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 <  2:2358308/100‑1 (MQ=255)
                                                    aTTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACAAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                               >  1:3298844/1‑100 (MQ=255)
                                                          cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  2:3025632/1‑100 (MQ=255)
                                                          cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCCTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:790243/1‑100 (MQ=255)
                                                            tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                       >  2:1684841/1‑100 (MQ=255)
                                                             atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCCGCCCTTGTTCGTGTTAATCCCAGCTCAGCATAACTTTCCCGGACTGGCCCGAACGCATAg                                      <  2:3817620/100‑1 (MQ=255)
                                                              tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     >  1:1597693/1‑100 (MQ=255)
                                                              tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     <  2:1007190/100‑1 (MQ=255)
                                                              tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     <  2:2103120/100‑1 (MQ=255)
                                                              tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     >  2:416700/1‑100 (MQ=255)
                                                                  aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 <  1:2297628/100‑1 (MQ=255)
                                                                  aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 >  2:2306753/1‑100 (MQ=255)
                                                                     gATTATCCTCCGATAAACGGGCTGCAATACCAGGCCTTATTCGTCTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCCAACGCCTAGCGGCaatg                              >  1:3097791/1‑98 (MQ=255)
                                                                     gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                              >  2:3107680/1‑100 (MQ=255)
                                                                      aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                             <  2:78694/100‑1 (MQ=255)
                                                                          tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTGTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                           <  1:1104532/98‑1 (MQ=255)
                                                                            cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                       >  2:1731942/1‑100 (MQ=255)
                                                                            cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                       >  2:1069717/1‑100 (MQ=255)
                                                                              cacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                     >  2:3263192/8‑100 (MQ=38)
                                                                              cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                     <  2:1498925/100‑1 (MQ=255)
                                                                               aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    >  1:538456/1‑100 (MQ=255)
                                                                               aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    <  1:279449/100‑1 (MQ=255)
                                                                               aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    <  2:2824695/100‑1 (MQ=255)
                                                                                gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   >  2:2516102/1‑100 (MQ=255)
                                                                                gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   >  2:405738/1‑100 (MQ=255)
                                                                                gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   <  2:1761912/100‑1 (MQ=255)
                                                                                cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                       >  1:1646913/6‑96 (MQ=39)
                                                                                 gaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  >  1:3387837/5‑100 (MQ=39)
                                                                                 atctaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  <  1:2726179/96‑1 (MQ=255)
                                                                                  tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAACTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 <  1:3494241/100‑1 (MQ=255)
                                                                                   acaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                >  1:2052614/3‑100 (MQ=255)
                                                                                    caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCGCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               >  1:242983/2‑100 (MQ=255)
                                                                                    aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               <  1:67498/100‑1 (MQ=255)
                                                                                     aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGTCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               >  2:1250600/1‑99 (MQ=255)
                                                                                     aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat              >  2:1626747/1‑100 (MQ=255)
                                                                                     aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTAAAAGCCCTTCTgca                    >  1:2814725/1‑92 (MQ=255)
                                                                                      aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             >  1:899965/1‑100 (MQ=255)
                                                                                      aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             >  1:2641022/1‑100 (MQ=255)
                                                                                      aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  1:1770991/100‑1 (MQ=255)
                                                                                      aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  1:543364/100‑1 (MQ=255)
                                                                                       gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATCGCCTCAAAGCCCTTCTCGAAATCATCg            >  1:3327419/1‑100 (MQ=255)
                                                                                         gatctaaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt          <  1:1191631/95‑1 (MQ=255)
                                                                                          cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa         >  1:3462290/1‑100 (MQ=255)
                                                                                          cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa         <  2:1370056/100‑1 (MQ=255)
                                                                                            ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga       >  1:3275658/1‑100 (MQ=255)
                                                                                            ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAACCATCGATaga       >  1:2345766/1‑100 (MQ=255)
                                                                                            ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTCCTGGAAATCATCGATaga       >  1:2462647/1‑100 (MQ=255)
                                                                                               aTACCAGTCTTTGTTCGTGTTAATCACAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATTGATAGCGaa    <  1:139114/100‑1 (MQ=255)
                                                                                                tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCACAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa   <  1:2028923/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790206‑3790402

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: