Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 18.1% 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,3030AT18.1% 486.5 / 94.6 204intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (86/81);  new base T (19/18);  total (105/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 207 total aligned reads displayed.

GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402
                                                                                                   |                                                                                                   
gagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATa                                                                                                     <  2:3434318/100‑1 (MQ=255)
  gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   <  2:2946118/100‑1 (MQ=255)
  gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   <  1:1983850/100‑1 (MQ=255)
  gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   >  1:2762934/1‑100 (MQ=255)
    agcTATATAAGAGATCAGCACAATTACGAACGCTTTATTACACTGACAGTATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCGGGAATAcc                                                                                                   <  1:2371803/98‑1 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTTCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  1:2223802/1‑100 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  2:3187879/1‑100 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                >  2:75966/1‑100 (MQ=255)
      cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                                 >  1:336281/1‑98 (MQ=255)
         tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                            <  2:311277/100‑1 (MQ=255)
         tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGAGTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                            <  1:558184/100‑1 (MQ=255)
          ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATACAAGGGCTGGTATTCCAGCCCttt                                                                                           <  2:4092123/100‑2 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                          <  1:1150411/100‑3 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                          <  1:2226635/100‑3 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                         <  1:4146804/100‑4 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                         <  2:3343709/100‑1 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTGATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                         <  1:2518690/100‑4 (MQ=255)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctg                                                                                     <  2:1915966/100‑8 (MQ=38)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     <  2:28287/100‑1 (MQ=255)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     >  1:976006/1‑100 (MQ=255)
                gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATACTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     >  1:2497081/1‑100 (MQ=255)
                 aTCAGCACTATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctg                                                                                     <  1:3294758/99‑8 (MQ=38)
                   cAGCACAATTACGAACTCTCTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  1:3654495/98‑1 (MQ=255)
                    aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                   <  2:1394802/98‑10 (MQ=38)
                    aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                 >  2:783735/1‑100 (MQ=255)
                      cacaATTACGAACTCTTTATTACACTGACAGTATTGTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                               <  1:2099936/100‑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)
                                aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                     >  1:2543383/1‑100 (MQ=255)
                                aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                     <  2:2200183/100‑1 (MQ=255)
                                 aCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                    <  2:2771111/100‑1 (MQ=255)
                                  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   >  2:2006187/1‑100 (MQ=255)
                                   tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                  >  1:4288484/1‑100 (MQ=255)
                                    ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                 >  1:2781134/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)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCCGCCCTTGTTCGTGTTCATCCCAGCTCAGAATAACTTTCCCg                                                         >  1:193506/1‑100 (MQ=255)
                                             acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                        <  2:2105389/100‑1 (MQ=255)
                                                gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                     >  2:2534965/1‑100 (MQ=255)
                                                gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                     >  2:687245/1‑100 (MQ=255)
                                                  naGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                   >  2:2348027/2‑100 (MQ=255)
                                                    ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 >  1:2860926/1‑100 (MQ=255)
                                                    ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                 <  2:3792432/100‑1 (MQ=255)
                                                    ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCCCAGAATAACTTTCCCGGACTgg                                                   >  1:3431926/1‑98 (MQ=255)
                                                        tttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGaa                                             <  1:3286491/100‑1 (MQ=255)
                                                         ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                            <  2:1609563/100‑1 (MQ=255)
                                                          atcttatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                           <  2:3827611/96‑1 (MQ=255)
                                                           aCATAATTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                           <  1:4025918/99‑1 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:4089896/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:2640289/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         <  1:2913029/100‑1 (MQ=255)
                                                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACgagat                                       >  1:3827587/1‑96 (MQ=255)
                                                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                       >  2:1519134/1‑100 (MQ=255)
                                                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                       >  2:1936179/1‑100 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATan                                      <  2:2480488/100‑2 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                      <  2:2439077/100‑1 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                      >  2:2996091/1‑100 (MQ=255)
                                                                tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     >  2:4394677/1‑100 (MQ=255)
                                                                tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     <  2:1275495/100‑1 (MQ=255)
                                                                    aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 >  1:2121600/1‑100 (MQ=255)
                                                                     cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                >  1:933898/1‑100 (MQ=255)
                                                                     cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                <  2:169847/100‑1 (MQ=255)
                                                                      aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                               >  2:1106188/1‑100 (MQ=255)
                                                                           aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                          <  2:3247259/100‑1 (MQ=255)
                                                                              cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                       >  2:797567/1‑100 (MQ=255)
                                                                               tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      <  1:3039109/100‑1 (MQ=255)
                                                                               tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  2:4231705/1‑100 (MQ=255)
                                                                               tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCGCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  1:1018609/1‑100 (MQ=255)
                                                                               acacgaacaaGGGCTGGTATTCCAGCCCTTTTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  1:3502677/9‑100 (MQ=37)
                                                                               acacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  1:1726792/9‑100 (MQ=38)
                                                                               acacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  1:1678219/9‑100 (MQ=38)
                                                                               acacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                      >  1:1771015/9‑100 (MQ=38)
                                                                                cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                     >  2:2727707/1‑100 (MQ=255)
                                                                                cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                     <  1:3816349/100‑1 (MQ=255)
                                                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   <  2:2860942/100‑1 (MQ=255)
                                                                                  gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCCGGaa                   >  1:1888576/1‑100 (MQ=255)
                                                                                   gaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  >  2:1245940/5‑100 (MQ=39)
                                                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  <  1:3216962/100‑1 (MQ=255)
                                                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 >  2:611814/1‑100 (MQ=255)
                                                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 <  1:3329449/100‑1 (MQ=255)
                                                                                    aacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 >  2:93455/4‑100 (MQ=255)
                                                                                     aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                <  2:336886/100‑1 (MQ=255)
                                                                                      caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               >  2:184306/2‑100 (MQ=255)
                                                                                       aaGGGCTGGTATTCCAGCGATTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCTCACAGGTCCGAACGCATCGCGTCAAAGCCCTTCTGCAAatcat              >  1:1576134/1‑100 (MQ=255)
                                                                                       aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat              >  2:3865993/1‑100 (MQ=255)
                                                                                       aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGCACGCATAGCGTCAAAGCCCTTCTGGAAatcat              >  2:48306/1‑100 (MQ=255)
                                                                                       aaGGGCTGGAAGACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat              <  1:135083/100‑1 (MQ=255)
                                                                                        aGGGCTGGAATACCAGCCCTTGTTCGTGTGAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             >  1:3957011/1‑100 (MQ=255)
                                                                                         gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg            >  2:101734/1‑100 (MQ=255)
                                                                                         gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg            <  2:1975287/100‑1 (MQ=255)
                                                                                          ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa           >  1:356515/1‑100 (MQ=255)
                                                                                           gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt          <  2:2121608/100‑1 (MQ=255)
                                                                                            ntGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa         >  2:3710790/2‑100 (MQ=255)
                                                                                               gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag      <  2:396025/100‑1 (MQ=255)
                                                                                                 aTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaa    >  2:1433569/1‑100 (MQ=255)
                                                                                                  tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa   <  2:1642515/100‑1 (MQ=255)
                                                                                                   aCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGACCGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc  >  1:3138229/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402

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