Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,300 A→T 15.7% intergenic (‑219/+20) waaH ← / ← tdh LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,3000AT15.7% 273.4 / 45.0 115intergenic (‑219/+20)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (56/41);  new base T (5/13);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-02

Only 100 of 119 total aligned reads displayed.

TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGA  >  NC_000913/3790201‑3790399
                                                                                                   |                                                                                                   
tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                     <  1:2046491/100‑1 (MQ=255)
  agagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAAt                                                                                                   <  2:323105/100‑1 (MQ=255)
    agcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAc                                                                                                 >  2:904161/1‑100 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCCGATAAAAGGGCTGGAATa                                                                                                  >  2:25645/1‑98 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                <  2:120423/100‑1 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                >  1:1572808/1‑100 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                >  1:825919/1‑100 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGAATTTACATATTTAAAAGATTATCATCAGCTAAAAGGGCTGGAATaac                                                                                                >  1:198123/1‑98 (MQ=255)
       agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                              <  2:825912/100‑1 (MQ=255)
        gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                             >  1:963138/1‑100 (MQ=255)
         cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                            <  2:848770/100‑1 (MQ=255)
          tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                            <  2:850786/99‑1 (MQ=255)
           atataAGAGATCAGCACAATTACGAACTCTTTATTATACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                          >  1:824646/1‑100 (MQ=255)
           atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                          <  1:1580309/100‑1 (MQ=255)
            tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                         <  1:1201637/100‑1 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                        >  1:600271/1‑100 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                       <  2:1056904/100‑3 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                       <  1:2020349/100‑3 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                      <  2:1404113/100‑4 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      >  2:356375/1‑100 (MQ=255)
               aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                      >  1:459619/1‑100 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  2:1701876/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                    <  2:167738/100‑1 (MQ=255)
                  agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatct                                                                                   <  1:1996785/100‑7 (MQ=38)
                   gaTCAGCTCAATTACAACCTCTTTATTACCCTTACATGATTTTACTTATTTAACAGATTATCCTCAGATAAAAGTGCTGGAATACCAGCCCTTGTTCgtg                                                                                  <  2:408668/100‑1 (MQ=255)
                   gaTCAGCACAATTACGAACTTTTTATTACACTGAACGGATTTTACATATTTAAAAGATTATCCTAAGGTAAAATGGCTGGAATACAAGcgattattc                                                                                     >  2:2014266/1‑89 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  <  1:1453971/100‑1 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  >  2:1184910/1‑100 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                  >  1:904165/1‑100 (MQ=255)
                    aTCAGCACAATTTCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                    <  2:100523/97‑6 (MQ=38)
                     tCAGCATAATTACGAACTCTTTATTACACTGACAGAATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                <  2:2262705/100‑10 (MQ=37)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                <  2:190830/100‑1 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                <  2:673569/100‑1 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                <  2:1322685/100‑1 (MQ=255)
                       ngCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                              >  1:1996587/2‑100 (MQ=255)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                              <  1:407767/100‑1 (MQ=255)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                              >  2:785284/1‑100 (MQ=255)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                              >  1:2281124/1‑100 (MQ=255)
                        gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAAt                                                                             >  2:87922/1‑100 (MQ=255)
                        gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAAt                                                                             >  1:597375/1‑100 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                           <  1:2398604/100‑1 (MQ=255)
                              ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAg                                                                        >  2:742601/1‑99 (MQ=255)
                              ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGc                                                                       <  2:786170/100‑1 (MQ=255)
                               tACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                      >  1:807858/1‑100 (MQ=255)
                                 cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                    >  1:1295968/1‑100 (MQ=255)
                                    aCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                 <  1:1532862/100‑1 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTAAGAATaa                                                              >  2:652674/1‑100 (MQ=255)
                                              acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                       >  1:430263/1‑100 (MQ=255)
                                              acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                       >  1:1345084/1‑100 (MQ=255)
                                               cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                      <  1:1590347/100‑1 (MQ=255)
                                               cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                      <  2:1886522/100‑1 (MQ=255)
                                               cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                      >  2:2223219/1‑100 (MQ=255)
                                               cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                      <  2:1886767/100‑1 (MQ=255)
                                                 cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                    >  2:20105/1‑100 (MQ=255)
                                                  tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   <  1:2223233/100‑1 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                   >  2:684648/1‑99 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                  <  1:1184918/100‑1 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                  >  1:1583752/1‑100 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                  >  1:161291/1‑100 (MQ=255)
                                                    aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTg                                                 <  2:503401/100‑1 (MQ=255)
                                                    aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTg                                                 >  2:2441719/1‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                              <  1:1339254/100‑1 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                    >  2:497420/1‑100 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                    >  1:59054/1‑100 (MQ=255)
                                                                   tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                  >  1:613626/1‑100 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCGTAGCg                                 <  1:1985924/100‑1 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTACCGGACTGGCCCGAACGCATAGCg                                 >  1:1190212/1‑100 (MQ=255)
                                                                     tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                 <  1:2022542/99‑1 (MQ=255)
                                                                      aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCTTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCTGGACTGGCCCGAACGCATAGCg                                 >  2:605529/1‑98 (MQ=255)
                                                                       aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                              <  2:751043/100‑1 (MQ=255)
                                                                         aGATTATCCTCAGATAACAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAACTCAAAATAACTTTCCCTGACTGGCCCGAACGCATAGCGTCaaa                            >  2:1098234/1‑100 (MQ=255)
                                                                         aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGATCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                            >  2:1642774/1‑100 (MQ=255)
                                                                          gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                           >  1:47989/1‑100 (MQ=255)
                                                                           aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                          <  1:2011099/100‑1 (MQ=255)
                                                                             tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                        >  1:101085/1‑100 (MQ=255)
                                                                              aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCt                       <  1:244177/100‑1 (MQ=255)
                                                                               tCCTCTGATAAAAGGGCTGGAATAACAGTCCTTGTTCGGGTTAAGCCCAGCTCAGAATAACCTTCCCGGACTGGCCCGAACGCATAGGGTCAAAGCCCtt                      <  2:618043/100‑1 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                     >  2:1891726/1‑99 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTg                   <  2:668231/100‑1 (MQ=255)
                                                                                   cacgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                  >  1:224410/8‑100 (MQ=38)
                                                                                   cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                  >  1:2495708/1‑100 (MQ=255)
                                                                                     gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                <  1:699223/100‑1 (MQ=255)
                                                                                     cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                >  1:427458/6‑100 (MQ=39)
                                                                                       tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat              >  2:2532457/1‑100 (MQ=255)
                                                                                        acaaGGGCAGGTATTCCAGCCCTTGTTGGTGTTAATCACAGCGCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGACCTTCTGGAAatc             >  2:897022/3‑100 (MQ=255)
                                                                                        aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc             <  2:163204/100‑1 (MQ=255)
                                                                                        aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc             <  2:1905780/100‑1 (MQ=255)
                                                                                         caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca            >  2:1312621/2‑100 (MQ=255)
                                                                                         aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca            <  1:1474997/100‑1 (MQ=255)
                                                                                         aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca            <  2:1211249/100‑1 (MQ=255)
                                                                                         aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca            <  2:2429153/100‑1 (MQ=255)
                                                                                         aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATACCAGCTCAGAATAACTTTCCCGGACTGGCACGCACGCATAGCGTCAAAGCCCTTCTGGCActca            >  2:2124524/1‑100 (MQ=255)
                                                                                          aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  1:1302849/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  1:1888575/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat           >  1:644285/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTCCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCAGAGCGTCACAGCCCTTCCgg                  >  2:2392373/1‑93 (MQ=255)
                                                                                          aaGGGATGGAATACCACCCCTTGTTCGTGTTAAACCAAACACAGAATAACTCGCCATGACTGGACCGAACGTATAGCGTCAAAGCGGTTCTgcgacgca            >  2:1503884/1‑92 (MQ=255)
                                                                                              gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt       <  2:2213762/100‑1 (MQ=255)
                                                                                                tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATag     >  1:174437/1‑100 (MQ=255)
                                                                                                   aaTAACACCCCTTCTTAGTGTTACTAACCGATCAGAATAAAATACCAGGAATGGACCGAACGCATAGAGTCAAAGCACTTCTGGAAATCATCGATAgaga  <  2:2137253/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: