Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 16.8% intergenic (‑176/+171) yjhU ← / ← yjhF putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,520,5000A.16.8% 343.8 / 69.9 125intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (59/45);  new base . (8/13);  total (67/58)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-02

Only 100 of 130 total aligned reads displayed.

TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520401‑4520599
                                                                                                   |                                                                                                    
tGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGgaga                                                                                                      <  1:2220584/100‑1 (MQ=255)
 ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  2:2627477/100‑1 (MQ=255)
 ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  2:493114/100‑1 (MQ=255)
  gCATGATTTTATACAGATATTTTATTCTTTTATATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  1:949600/99‑1 (MQ=255)
  gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    >  1:1271599/1‑100 (MQ=255)
   cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   <  1:344783/100‑1 (MQ=255)
       aTTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACAATTTGAGCGATAGCAAGCTGTTACCCCATACTGTAAAGCCGGAGACATg                                                                                                  >  2:2447217/1‑97 (MQ=255)
          ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCgg                                                                                           <  2:2911167/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                          <  2:1856768/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGc                                                                                          <  2:794034/100‑1 (MQ=39)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                        <  2:795995/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                        <  2:104223/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                        <  2:796265/100‑1 (MQ=38)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                        <  1:2215314/100‑1 (MQ=38)
              aCATATATTTTATTCTTTTACATGAAATAAACCATCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                       <  2:929328/100‑1 (MQ=37)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                       <  1:1704532/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                       <  1:1520644/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                       <  2:325829/100‑1 (MQ=38)
                 tcttccgatctCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                    <  1:1950119/90‑1 (MQ=25)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                    <  1:567151/100‑1 (MQ=38)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                    <  1:1737504/100‑1 (MQ=38)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATAACATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                     <  1:34933/99‑1 (MQ=255)
                  atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                   <  1:910657/100‑1 (MQ=38)
                  atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                   <  1:1457707/100‑1 (MQ=38)
                    atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                 >  1:908986/1‑100 (MQ=255)
                    actttATTCTTTTACATGAAATAAACCAGCATACCATTTGGTCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCTGATACATGCTCCGGCTTTACAGt                                                                                  <  2:1481318/97‑1 (MQ=255)
                     ttttATTCTTTTACATGAAATAAACCAGCATGCCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGTCGTAGACATGCTCCGGCTTTACAGTa                                                                                 <  2:2682777/99‑1 (MQ=255)
                          ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACt                                                                            >  2:2184237/1‑99 (MQ=255)
                          ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                           >  2:1529746/1‑100 (MQ=255)
                           tCTTTTACATGAAATCAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTga                                                                          <  1:3000453/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAGATATGCTCCGGCTTTACAGTATTACTgag                                                                         <  2:702473/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                         >  1:122035/1‑100 (MQ=255)
                             ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                        <  1:2824730/100‑1 (MQ=255)
                             ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                        <  1:2840008/100‑1 (MQ=255)
                             ttttACATGAAATAAACCAGAATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                        >  2:161558/1‑100 (MQ=255)
                              tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                       <  2:352656/100‑1 (MQ=255)
                              tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAAACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                       >  1:502247/1‑100 (MQ=255)
                                taactGAACTAAACCAGCATACCATTTGTGCGATCGCCCGCTGTTGCCCCATACTGTAAAGCGTTAGACATGCTCCGGTTTTACAGTATGACTgag                                                                         <  2:2904867/92‑1 (MQ=255)
                                 aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATac                                                                    >  1:1812591/1‑100 (MQ=255)
                                  cATGAAATAAACCAGTATACCATTTGAGCGATCGCTCACTGTTCCCCCATACTGTAAATCCTGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  2:2660536/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACGGTATTACTGAGAATaca                                                                   >  1:2170524/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  2:959972/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCACCATACTGTAAAGCCGGAGACCTGCTACGGCTTTACAGTATTATTGAGAATcca                                                                   >  2:245667/1‑100 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAATCCGGAGACATGCTCCAGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1136347/100‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                  <  2:407344/99‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCGCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1996617/100‑1 (MQ=255)
                                     gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                >  2:2073817/1‑100 (MQ=255)
                                      aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               >  1:2225646/1‑100 (MQ=255)
                                      aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               >  2:2254686/1‑100 (MQ=255)
                                       aaTACACCAGCATACCATTTGAGCGATCGCCCGCTGTTCACCCATACTGTAAAGCAGGAGACATACTCCGCCTTTACAGTATTACAGAAAAAACACCCAt                                                              >  2:2154147/1‑100 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGGTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          <  1:2219203/100‑1 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  1:3033402/1‑100 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  1:2331508/1‑100 (MQ=255)
                                            aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCAGACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                         <  2:2331481/100‑1 (MQ=255)
                                             ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                        <  1:245670/100‑1 (MQ=255)
                                                  aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                   >  1:1496197/1‑100 (MQ=255)
                                                   tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  >  1:1082136/1‑100 (MQ=255)
                                                   tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  <  1:2195176/100‑1 (MQ=255)
                                                    aCCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  <  2:2225623/99‑1 (MQ=255)
                                                    aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                 >  2:673321/1‑100 (MQ=255)
                                                     ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                >  1:2789301/1‑100 (MQ=255)
                                                       aTTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                              <  2:533105/100‑1 (MQ=255)
                                                        tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                >  2:1306386/1‑97 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:2706918/100‑1 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:2618685/100‑1 (MQ=255)
                                                            aGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                         <  2:230406/100‑1 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   >  2:1217750/1‑100 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAGTACACTCATCAAGCATAATCGAAGGTTCAATATCGc                                   <  2:2114975/99‑1 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  <  1:625068/100‑1 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  2:943851/1‑100 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  2:2905147/1‑100 (MQ=255)
                                                                    cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 <  1:902725/100‑1 (MQ=255)
                                                                      cGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                               <  1:1325952/100‑1 (MQ=255)
                                                                         ttctcagtaaTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                            <  1:1573612/91‑1 (MQ=21)
                                                                              ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       >  1:2245500/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  1:2689461/100‑1 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  1:3019787/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  2:1705731/1‑100 (MQ=255)
                                                                                ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGACCCTTTTTGATCGc                     >  2:2067240/1‑100 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    >  1:2907344/1‑100 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    <  2:2265585/100‑1 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    >  1:81511/1‑100 (MQ=255)
                                                                                   tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  1:2169233/1‑100 (MQ=255)
                                                                                    aCTGTAAAGCCGGAGCGTGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGAATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  1:586548/1‑99 (MQ=37)
                                                                                     cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  1:23834/1‑100 (MQ=37)
                                                                                      tGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa               >  2:170106/1‑100 (MQ=255)
                                                                                       gTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCACCAGCGa               >  1:1197632/1‑99 (MQ=37)
                                                                                        tAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt             >  2:1474302/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:3043967/1‑100 (MQ=255)
                                                                                         aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGATTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:2397734/1‑100 (MQ=255)
                                                                                         aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCAAAGGTTCAATATCGCAACTTTTTGATCGCCATCGAtt             >  1:1362652/1‑99 (MQ=255)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:2977919/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGCATGTCTACAGCTGAAAAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTACAATATCGGCCCTTTTTGATCgaaagcaatac            >  2:2960378/1‑90 (MQ=255)
                                                                                            gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCATTTTTGATCGCCAGCGATTCGTc         >  1:64006/1‑100 (MQ=255)
                                                                                            gCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         >  1:2568817/1‑100 (MQ=39)
                                                                                              cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc       >  2:1778533/1‑100 (MQ=255)
                                                                                              cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGCGATTCGTCAc       >  2:1165096/1‑100 (MQ=255)
                                                                                                  gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTTATTGTTATTGTTTCGTCACCGTg   >  2:2896292/1‑100 (MQ=255)
                                                                                                   aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  1:803074/100‑1 (MQ=255)
                                                                                                   aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  1:2978190/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520401‑4520599

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