Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 15.5% 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,3030AT15.5% 289.5 / 43.2 116intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (56/42);  new base T (5/13);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 3.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 118 total aligned reads displayed.

AGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAA  >  NC_000913/3790205‑3790401
                                                                                                  |                                                                                                  
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:404785/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)
       atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                            <  1:1580309/100‑1 (MQ=255)
        tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                           <  1:1201637/100‑1 (MQ=255)
         ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                          >  1:723045/1‑100 (MQ=255)
         ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                          >  1:600271/1‑100 (MQ=255)
          taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                         <  1:2020349/100‑3 (MQ=255)
           aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                        <  1:2231717/100‑4 (MQ=255)
           aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttta                                                                                        <  2:1404113/100‑4 (MQ=255)
           aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                        >  2:356375/1‑100 (MQ=255)
           aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                        >  1:459619/1‑100 (MQ=255)
           aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCATTGtt                                                                                        >  2:1890174/1‑100 (MQ=255)
           aaGACATCAGCACAATTACAAACTCTTTATTACACTTACATGATTTTACATATTTAACATATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGtt                                                                                        <  2:1207391/100‑1 (MQ=255)
             gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                      <  2:167738/100‑1 (MQ=255)
              agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatct                                                                                     <  1:1996785/100‑7 (MQ=38)
              agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                     <  2:261067/100‑1 (MQ=255)
               gaTCAGCTCAATTACAACCTCTTTATTACCCTTACATGATTTTACTTATTTAACAGATTATCCTCAGATAAAAGTGCTGGAATACCAGCCCTTGTTCgtg                                                                                    <  2:408668/100‑1 (MQ=255)
               gaTCAGCACAATTACGAACTTTTTATTACACTGAACGGATTTTACATATTTAAAAGATTATCCTAAGGTAAAATGGCTGGAATACAAGcgattattc                                                                                       >  2:2014266/1‑89 (MQ=255)
               gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                    >  1:904165/1‑100 (MQ=255)
               gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                    <  1:1453971/100‑1 (MQ=255)
               gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                    >  2:1184910/1‑100 (MQ=255)
                aTCAGCACAATTTCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                      <  2:100523/97‑6 (MQ=38)
                aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctga                                                                                   <  1:976432/100‑9 (MQ=38)
                 tCAGCATAATTACGAACTCTTTATTACACTGACAGAATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                                  <  2:2262705/100‑10 (MQ=37)
                 tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                  <  2:673569/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)
                      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)
                                 ctctTTATTACACTGCCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTTTTCGTGTTACTCCCAGCTCa                                                                      >  2:2047662/1‑96 (MQ=255)
                                   ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTAAGAATaa                                                                >  2:652674/1‑100 (MQ=255)
                                          acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                         >  1:1345084/1‑100 (MQ=255)
                                          acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                         >  1:430263/1‑100 (MQ=255)
                                           cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                        <  2:1886522/100‑1 (MQ=255)
                                           cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                        <  2:1886767/100‑1 (MQ=255)
                                           cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                        <  1:1590347/100‑1 (MQ=255)
                                             cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                      >  2:20105/1‑100 (MQ=255)
                                              tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                     <  1:2223233/100‑1 (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:2441719/1‑100 (MQ=255)
                                                aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTg                                                   <  2:503401/100‑1 (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)
                                                                ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                   >  2:2296035/1‑100 (MQ=255)
                                                                ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTACCGGACTGGCCCGAACGCATAGCg                                   >  1:1190212/1‑100 (MQ=255)
                                                                 tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                   <  1:2022542/99‑1 (MQ=255)
                                                                  aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                 >  2:690083/1‑100 (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)
                                                                       aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                            <  2:1572797/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)
                                                                               cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                    >  1:2495708/1‑100 (MQ=255)
                                                                                 gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                  <  1:699223/100‑1 (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)
                                                                                               aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga    >  2:61409/1‑100 (MQ=255)
                                                                                               aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga    <  1:565872/100‑1 (MQ=255)
                                                                                               aaTAACACCCCTTCTTAGTGTTACTAACCGATCAGAATAAAATACCAGGAATGGACCGAACGCATAGAGTCAAAGCACTTCTGGAAATCATCGATAgaga    <  2:2137253/100‑1 (MQ=255)
                                                                                                 tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa  >  1:1382835/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAA  >  NC_000913/3790205‑3790401

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: