Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,300 A→T 18.5% 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,3000AT18.5% 457.1 / 97.7 201intergenic (‑219/+20)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (84/79);  new base T (19/18);  total (103/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

Only 100 of 207 total aligned reads displayed.

TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGA  >  NC_000913/3790201‑3790399
                                                                                                   |                                                                                                   
tgaatgtatcTATATAAGAGATCTGTACAATTATGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     <  1:4165771/91‑1 (MQ=255)
tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     >  2:2516377/1‑100 (MQ=255)
tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     >  2:3387494/1‑100 (MQ=255)
tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     >  1:2446101/1‑100 (MQ=255)
   gagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATa                                                                                                  <  2:3434318/100‑1 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                >  1:2762934/1‑100 (MQ=255)
      cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                               >  2:1686630/1‑100 (MQ=255)
        gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                             >  2:75966/1‑100 (MQ=255)
         cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                            <  1:4104100/100‑1 (MQ=255)
          tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATGTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                             <  1:4039801/98‑1 (MQ=255)
          tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGAGTTTACATATTTAGCAGATTATCCGCAGATAAAAGGGCTGGAATACCAGcc                                                                                            <  1:2037684/99‑1 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATACAAGGGCTGGTATTCCAGCCCttt                                                                                        <  2:4092123/100‑2 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                        <  1:3542017/100‑2 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                        >  1:1941986/1‑100 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTTTTCCAGCCCtttt                                                                                       <  2:2957709/100‑3 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                       <  1:1150411/100‑3 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                       >  2:2926607/1‑100 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                       >  2:1770774/1‑100 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                      <  2:3103122/100‑4 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      <  2:2583497/100‑1 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      <  2:1907942/100‑1 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTGATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                      <  1:2518690/100‑4 (MQ=255)
                agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                     <  2:3222599/100‑5 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACTCTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  2:282264/100‑6 (MQ=38)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    >  2:3296591/1‑100 (MQ=255)
                  agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGAGCTGGAATACCAGCCCTTCTTTgt                                                                                   >  1:3360146/1‑100 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  >  1:976006/1‑100 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  <  1:1770768/100‑1 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  <  1:2221859/100‑1 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATACTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  >  1:2497081/1‑100 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                <  1:2348856/100‑1 (MQ=255)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                <  2:1394802/98‑10 (MQ=38)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                              >  2:783735/1‑100 (MQ=255)
                          acctttACGAACCCTTTATTACACTGACAGGATTGTACATATTTAACAGCTTATCCTCAGATAAAACGGCTGGACTACCGGCCCTTGTTCGTGTTAATc                                                                            <  1:3625928/95‑1 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                           <  2:4351076/100‑1 (MQ=255)
                           cgaTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                          >  2:108926/3‑100 (MQ=255)
                           caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                          >  2:1076818/1‑100 (MQ=255)
                              ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGc                                                                       <  2:2228533/100‑1 (MQ=255)
                               tatattcgttttATTACTCTGACACGAGTTTATATTTTTAACAGGTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCa                                                                         <  1:4001937/88‑1 (MQ=255)
                               tACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                      >  1:856269/1‑100 (MQ=255)
                                 cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                    >  1:2224910/1‑100 (MQ=255)
                                  gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAg                                                                   <  1:1454436/100‑1 (MQ=255)
                                   aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                  >  1:2543383/1‑100 (MQ=255)
                                   aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                  <  2:2200183/100‑1 (MQ=255)
                                      tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                               >  1:4288484/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                              >  1:2781134/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                              >  2:3363600/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                              >  1:19253/1‑100 (MQ=255)
                                        tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                             <  2:2280362/100‑1 (MQ=255)
                                           attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                          >  1:495326/1‑100 (MQ=255)
                                           attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAACAGCCTTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                          >  1:893817/1‑100 (MQ=255)
                                                acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                     <  2:2105389/100‑1 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                  >  2:687245/1‑100 (MQ=255)
                                                    aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACTAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   >  2:2851290/1‑98 (MQ=255)
                                                     naGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                >  2:2348027/2‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              >  1:2860926/1‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              <  1:3808244/100‑1 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              <  2:3792432/100‑1 (MQ=255)
                                                            ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                         <  2:1609563/100‑1 (MQ=255)
                                                             atcttatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                        <  2:3827611/96‑1 (MQ=255)
                                                               cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                      <  1:2913029/100‑1 (MQ=255)
                                                               cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                      >  1:2640289/1‑100 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                    >  1:3498489/1‑100 (MQ=255)
                                                                  atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATan                                   <  2:2480488/100‑2 (MQ=255)
                                                                  atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                   <  2:2439077/100‑1 (MQ=255)
                                                                   tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                  >  2:4394677/1‑100 (MQ=255)
                                                                   tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                  <  2:1275495/100‑1 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTAGCCCGAACGCATAGc                                  >  1:3288970/1‑99 (MQ=255)
                                                                     tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGt                                <  1:3296226/100‑1 (MQ=255)
                                                                       aCAGATTATCCTCAGATAAAAGCGCTGGCATACCCTCCCTTGTTCGTGTCAATCCCAGCTCAGAATAACTTTCCCGGACTGTCCCGAACGCATATCGTc                               <  1:3218414/99‑1 (MQ=255)
                                                                        cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                             >  1:933898/1‑100 (MQ=255)
                                                                         aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                            >  2:1106188/1‑100 (MQ=255)
                                                                         aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                            >  1:1669694/1‑100 (MQ=255)
                                                                                ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTCCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                     <  1:3415463/100‑1 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                    >  2:797567/1‑100 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   <  2:1988140/100‑1 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   >  2:4231705/1‑100 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   <  1:3039109/100‑1 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   <  1:1702985/100‑1 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCGCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   >  1:1018609/1‑100 (MQ=255)
                                                                                    aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                 <  2:1932511/100‑1 (MQ=255)
                                                                                    aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                 <  1:2108404/100‑1 (MQ=255)
                                                                                     gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                >  1:1784756/1‑100 (MQ=255)
                                                                                      gaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa               >  2:1245940/5‑100 (MQ=39)
                                                                                      aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa               <  1:3876437/100‑1 (MQ=255)
                                                                                       tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat              >  2:611814/1‑100 (MQ=255)
                                                                                       aacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat              >  2:93455/4‑100 (MQ=255)
                                                                                         caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAAGCATAGCGTAAAAGCCCTTCTGGAAat              >  1:3827046/2‑98 (MQ=255)
                                                                                          aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  2:3865993/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGCACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  2:48306/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAAGACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           <  1:135083/100‑1 (MQ=255)
                                                                                           aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc          >  2:2742564/1‑100 (MQ=255)
                                                                                            gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg         >  2:101734/1‑100 (MQ=255)
                                                                                              gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGCATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt       >  2:1517169/1‑100 (MQ=255)
                                                                                              gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt       <  2:2121608/100‑1 (MQ=255)
                                                                                               ntGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa      >  2:3710790/2‑100 (MQ=255)
                                                                                               cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa      >  1:3179333/1‑100 (MQ=255)
                                                                                                tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATag     <  2:31863/100‑1 (MQ=255)
                                                                                                 ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga    >  2:194085/1‑100 (MQ=255)
                                                                                                  gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag   <  2:396025/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: