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.

TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  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)
          ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCgg                                                                                           <  2:2911167/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                          <  2:1856768/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                          <  2:794034/100‑1 (MQ=39)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                        <  2:795995/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                        <  2:104223/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  2:796265/100‑1 (MQ=38)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  1:2215314/100‑1 (MQ=38)
              aCATATATTTTATTCTTTTACATGAAATAAACCATCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  2:929328/100‑1 (MQ=37)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                       <  1:1704532/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                       <  1:1520644/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  2:325829/100‑1 (MQ=38)
                 tcttccgatctCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                    <  1:1950119/90‑1 (MQ=25)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                    <  1:567151/100‑1 (MQ=38)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                    <  1:1737504/100‑1 (MQ=38)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATAACATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  1:34933/99‑1 (MQ=255)
                  atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                   <  1:910657/100‑1 (MQ=38)
                  atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                   <  1:1457707/100‑1 (MQ=38)
                    atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                 >  1:908986/1‑100 (MQ=255)
                    actttATTCTTTTACATGAAATAAACCAGCATACCATTTGGTCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCTGATACATG‑CTCCGGCTTTACAGt                                                                                  <  2:1481318/97‑1 (MQ=255)
                     ttttATTCTTTTACATGAAATAAACCAGCATGCCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGTCGTAGACATG‑CTCCGGCTTTACAGTa                                                                                 <  2:2682777/99‑1 (MQ=255)
                          ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACt                                                                            >  2:2184237/1‑99 (MQ=255)
                          ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                           >  2:1529746/1‑100 (MQ=255)
                           tCTTTTACATGAAATCAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                          <  1:3000453/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAGATATG‑CTCCGGCTTTACAGTATTACTgag                                                                         <  2:702473/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         >  1:122035/1‑100 (MQ=255)
                             ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        <  1:2824730/100‑1 (MQ=255)
                             ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        <  1:2840008/100‑1 (MQ=255)
                             ttttACATGAAATAAACCAGAATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        >  2:161558/1‑100 (MQ=255)
                              tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                       <  2:352656/100‑1 (MQ=255)
                              tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAAACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                       >  1:502247/1‑100 (MQ=255)
                                taactGAACTAAACCAGCATACCATTTGTGCGATCGCCCGCTGTTGCCCCATACTGTAAAGCGTTAGACATG‑CTCCGGTTTTACAGTATGACTgag                                                                         <  2:2904867/92‑1 (MQ=255)
                                 aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                    >  1:1812591/1‑100 (MQ=255)
                                  cATGAAATAAACCAGTATACCATTTGAGCGATCGCTCACTGTTCCCCCATACTGTAAATCCTGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  2:2660536/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACGGTATTACTGAGAATaca                                                                   >  1:2170524/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  2:959972/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCACCATACTGTAAAGCCGGAGACCTG‑CTACGGCTTTACAGTATTATTGAGAATcca                                                                   >  2:245667/1‑100 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAATCCGGAGACATG‑CTCCAGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1136347/100‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATacac                                                                  <  2:407344/99‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CGCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1996617/100‑1 (MQ=255)
                                     gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                >  2:2073817/1‑100 (MQ=255)
                                      aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               >  1:2225646/1‑100 (MQ=255)
                                      aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               >  2:2254686/1‑100 (MQ=255)
                                       aaTACACCAGCATACCATTTGAGCGATCGCCCGCTGTTCACCCATACTGTAAAGCAGGAGACATA‑CTCCGCCTTTACAGTATTACAGAAAAAACACCCAt                                                              >  2:2154147/1‑100 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGGTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          <  1:2219203/100‑1 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  1:3033402/1‑100 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  1:2331508/1‑100 (MQ=255)
                                            aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCAGACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                         <  2:2331481/100‑1 (MQ=255)
                                             ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                        <  1:245670/100‑1 (MQ=255)
                                                  aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                   >  1:1496197/1‑100 (MQ=255)
                                                   tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  >  1:1082136/1‑100 (MQ=255)
                                                   tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  <  1:2195176/100‑1 (MQ=255)
                                                    aCCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  <  2:2225623/99‑1 (MQ=255)
                                                    aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                 >  2:673321/1‑100 (MQ=255)
                                                     ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                >  1:2789301/1‑100 (MQ=255)
                                                       aTTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                              <  2:533105/100‑1 (MQ=255)
                                                        tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                >  2:1306386/1‑97 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:2706918/100‑1 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:2618685/100‑1 (MQ=255)
                                                            aGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                         <  2:230406/100‑1 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   >  2:1217750/1‑100 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAGTACACTCATCAAGCATAATCGAAGGTTCAATATCGc                                   <  2:2114975/99‑1 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  <  1:625068/100‑1 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  2:943851/1‑100 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  2:2905147/1‑100 (MQ=255)
                                                                    cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATT‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 <  1:902725/100‑1 (MQ=255)
                                                                      cGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                               <  1:1325952/100‑1 (MQ=255)
                                                                         ttctcagtaaTACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                            <  1:1573612/91‑1 (MQ=21)
                                                                              ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       >  1:2245500/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  1:2689461/100‑1 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  1:3019787/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  2:1705731/1‑100 (MQ=255)
                                                                                ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGACCCTTTTTGATCGc                     >  2:2067240/1‑100 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    >  1:2907344/1‑100 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    <  2:2265585/100‑1 (MQ=255)
                                                                                 cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    >  1:81511/1‑100 (MQ=255)
                                                                                   tACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  1:2169233/1‑100 (MQ=255)
                                                                                    aCTGTAAAGCCGGAG‑CGTGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGAATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  1:586548/1‑99 (MQ=37)
                                                                                     cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  1:23834/1‑100 (MQ=37)
                                                                                      tGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa               >  2:170106/1‑100 (MQ=255)
                                                                                       gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCACCAGCGa               >  1:1197632/1‑99 (MQ=37)
                                                                                        tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt             >  2:1474302/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:3043967/1‑100 (MQ=255)
                                                                                         aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGATTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:2397734/1‑100 (MQ=255)
                                                                                         aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCAAAGGTTCAATATCGCAACTTTTTGATCGCCATCGAtt             >  1:1362652/1‑99 (MQ=255)
                                                                                         aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:2977919/1‑100 (MQ=38)
                                                                                         aaaGCCGGAG‑CATGTCTACAGCTGAAAAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTACAATATCGGCCCTTTTTGATCgaaagcaatac            >  2:2960378/1‑90 (MQ=255)
                                                                                            gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCATTTTTGATCGCCAGCGATTCGTc         >  1:64006/1‑100 (MQ=255)
                                                                                            gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         >  1:2568817/1‑100 (MQ=39)
                                                                                              cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc       >  2:1778533/1‑100 (MQ=255)
                                                                                              cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGCGATTCGTCAc       >  2:1165096/1‑100 (MQ=255)
                                                                                                  gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTTATTGTTATTGTTTCGTCACCGTg   >  2:2896292/1‑100 (MQ=255)
                                                                                                   aCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  1:803074/100‑1 (MQ=255)
                                                                                                   aCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  1:2978190/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  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: ‑