Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 16.0% 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.0% 573.2 / 107.9 200intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (82/86);  new base . (18/14);  total (100/100)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 208 total aligned reads displayed.

TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520401‑4520599
                                                                                                   |                                                                                                    
tGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTATAAAGCCGgaga                                                                                                      >  1:2062054/1‑100 (MQ=255)
 ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  2:1887206/100‑1 (MQ=255)
   cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  2:85292/1‑100 (MQ=255)
   cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  1:617544/1‑100 (MQ=255)
    aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATg                                                                                                  <  1:1879751/100‑1 (MQ=255)
      gATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑Ct                                                                                               >  2:3498331/1‑100 (MQ=255)
        ttttATACAGATATTTTATTCTTTTACATGAAATAACCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTcc                                                                                             <  2:1191976/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                          >  1:1365322/1‑100 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                          <  2:3907422/100‑1 (MQ=255)
            ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCt                                                                                         <  1:3320703/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  2:567090/100‑1 (MQ=38)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  2:1935271/100‑1 (MQ=38)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                       >  2:3984931/1‑100 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                       <  2:831811/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  1:3474015/100‑1 (MQ=38)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCAATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  1:3191638/100‑1 (MQ=38)
               cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                      <  1:683985/100‑1 (MQ=255)
               cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                      >  1:2078443/1‑100 (MQ=255)
               cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                      >  1:737139/1‑100 (MQ=255)
                aGATATTTTATTCTTTTACATGACATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTGCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                      <  1:2712336/99‑1 (MQ=255)
                aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  2:1752692/100‑1 (MQ=255)
                aGATATTTTATTCTATTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  1:2908898/100‑1 (MQ=255)
                 gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCTATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATA‑CTCCGGCTTTACa                                                                                    >  1:1000322/1‑100 (MQ=255)
                  atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCTATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                   <  1:788404/100‑1 (MQ=37)
                    atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCTGGCTTTACAGTa                                                                                 <  2:2813794/100‑1 (MQ=255)
                     ttttATTCTTTTACATGAAATAAACCAGCAAACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCTGGCTTTACAGTAt                                                                                >  1:2809342/1‑100 (MQ=255)
                        tATTCGTTTACATGACATAAACCATCATACCATGTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                             <  1:3924553/100‑1 (MQ=39)
                          ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                           <  2:3334990/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         >  1:4160870/1‑100 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         >  1:3249561/1‑100 (MQ=255)
                                  naTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  1:2089937/2‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  2:1861375/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  1:4052353/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  1:2142451/1‑100 (MQ=255)
                                  cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  1:2579690/1‑100 (MQ=255)
                                     gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                >  2:4269183/1‑100 (MQ=255)
                                       aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                              >  2:1618851/1‑100 (MQ=255)
                                         tAAACCAGCATACCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATc                                                             <  1:1628520/99‑1 (MQ=255)
                                          aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                           >  2:1892727/1‑100 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          <  1:1274455/100‑1 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          <  2:2389444/100‑1 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  2:3283160/1‑100 (MQ=255)
                                             ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                        <  2:1634311/100‑1 (MQ=255)
                                               aGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATa                                                      <  2:2081471/100‑1 (MQ=255)
                                                gCATACCATTTGAGCGATCGCTCGCTGTTTCCGCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                       >  1:1731611/1‑98 (MQ=255)
                                                  aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                   <  1:1760567/100‑1 (MQ=255)
                                                   tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  >  2:3846010/1‑100 (MQ=255)
                                                    aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                   <  1:715520/98‑1 (MQ=255)
                                                     ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                <  2:1106740/100‑1 (MQ=255)
                                                     ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                >  2:2507580/1‑100 (MQ=255)
                                                       aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                              <  2:2307736/100‑1 (MQ=255)
                                                         ttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGtt                                            >  1:2689163/1‑100 (MQ=255)
                                                          tctgCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:2490048/97‑1 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATTCCCCCATCAAGCATAATAGGAGGtt                                            >  2:3355051/1‑99 (MQ=38)
                                                           gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATATACCCATCAAGCATAATCGAAGGTTCa                                          <  2:2603425/100‑1 (MQ=255)
                                                             gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                        >  1:748173/1‑100 (MQ=255)
                                                             gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                        >  2:3360445/1‑100 (MQ=255)
                                                               ggtCGCCCGCTGTGCACACATACGGTAAAGCCAGAGACATG‑CGACGGCTTTACAGACTTACTCAGAATAGACCCATCCAGCATAATCGAAGGTCCAat                                        >  1:2174326/3‑98 (MQ=255)
                                                               gagcGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                          <  1:1240412/93‑1 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   <  2:4160892/100‑1 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAGG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   <  1:342037/100‑1 (MQ=255)
                                                                          gTTCCCCCATACTGTAAGGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAAGACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           >  1:194257/1‑100 (MQ=255)
                                                                          gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTCTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           <  1:3442669/100‑1 (MQ=255)
                                                                           ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGTTTCAATATCGCCCCTTTTTg                          <  1:2603506/100‑1 (MQ=255)
                                                                           ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          <  1:3586676/100‑1 (MQ=255)
                                                                           ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          <  1:3726704/100‑1 (MQ=255)
                                                                           ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          <  2:2166542/100‑1 (MQ=255)
                                                                            tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                         >  2:3185901/1‑100 (MQ=255)
                                                                             cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                        >  2:2228182/1‑100 (MQ=255)
                                                                              ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       >  1:1787042/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCACGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  1:3022573/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  1:578943/100‑1 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  2:3296825/100‑1 (MQ=255)
                                                                                    aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTGTTTGATCGCCAGc                 <  2:2901342/100‑1 (MQ=255)
                                                                                     cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCTTCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  <  2:3643977/98‑1 (MQ=255)
                                                                                     cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATACTCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 >  1:3817463/1‑99 (MQ=255)
                                                                                     cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 >  1:846546/1‑99 (MQ=38)
                                                                                     cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  2:708420/1‑100 (MQ=38)
                                                                                     cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  2:1226003/1‑100 (MQ=38)
                                                                                      tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGATCAATATCGCCCCTTTTTGATCGCCAGCGa               >  2:2356903/1‑100 (MQ=37)
                                                                                      tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAAAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATAGCCAGCGa               >  1:2744510/1‑100 (MQ=37)
                                                                                       gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCCATATCGCCCCTTTTTGATCCCCCGc                 >  1:2403821/1‑97 (MQ=255)
                                                                                       gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  1:906494/1‑100 (MQ=255)
                                                                                       gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  1:1207936/1‑100 (MQ=38)
                                                                                        tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt             <  2:617546/100‑1 (MQ=255)
                                                                                        tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCGATATCGCCCCTTTTTGATCGCCAGCGAtt             >  1:1528068/1‑100 (MQ=38)
                                                                                        tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  1:2946064/1‑99 (MQ=38)
                                                                                         aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:324339/1‑100 (MQ=38)
                                                                                         aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAACGAAGGTTCAATATCGACCCTTTTTGATCGCCAGCGACTc            >  2:1627107/1‑100 (MQ=37)
                                                                                          aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           <  2:1048786/100‑1 (MQ=255)
                                                                                          aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           >  2:1560435/1‑100 (MQ=255)
                                                                                           aGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt          >  1:1227775/1‑100 (MQ=38)
                                                                                            gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         >  2:1379197/1‑100 (MQ=255)
                                                                                            gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         <  1:1084353/100‑1 (MQ=255)
                                                                                              cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc       <  2:1507283/100‑1 (MQ=255)
                                                                                              cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc       <  2:3255175/100‑1 (MQ=255)
                                                                                               ggAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAcc      >  2:3434533/1‑100 (MQ=255)
                                                                                                gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCg     <  1:486399/100‑1 (MQ=255)
                                                                                                 ataCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCg     <  1:512593/97‑1 (MQ=255)
                                                                                                  gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg   >  2:4109241/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑