Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 19.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.19.0% 408.8 / 100.0 153intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (69/51);  new base . (18/11);  total (91/62)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

Only 100 of 159 total aligned reads displayed.

GGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520402‑4520599
                                                                                                  |                                                                                                    
ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCACCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  2:285678/100‑1 (MQ=255)
 gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    >  2:1189240/1‑100 (MQ=255)
 gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    <  2:2318194/100‑1 (MQ=255)
  cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  1:910669/1‑100 (MQ=255)
   aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATg                                                                                                  >  1:2285774/1‑100 (MQ=255)
    tGATTTTATACAGATATTTTATTCTTTTACATGTAATAAACCAGCATACCATTTGAGCGATCGTCCGCTGTTCCCCCATACTCTAAAGCCGGAGACATg                                                                                                  >  2:1532697/1‑99 (MQ=255)
    tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑c                                                                                                >  1:1794905/1‑100 (MQ=255)
       ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTcc                                                                                             <  1:717114/100‑1 (MQ=255)
          tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                          <  2:3401765/100‑1 (MQ=255)
           ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                         <  1:216044/100‑1 (MQ=38)
            taCAGATATTTTATTCTTTTACATGATATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  2:2946106/100‑1 (MQ=38)
            taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  2:639265/100‑1 (MQ=38)
            taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                        <  1:2088751/100‑1 (MQ=38)
             tcAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  1:2805745/99‑1 (MQ=38)
             aCAGGTATTTTATTCTTTTACGTGAAATATACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCGTACTGTAAAGCCGGAGATATG‑CTCCGGCttt                                                                                       <  2:3863146/100‑1 (MQ=255)
             aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  2:3475708/100‑1 (MQ=38)
             aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  2:3927765/100‑1 (MQ=38)
              cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                      >  2:162040/1‑100 (MQ=255)
              cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTa                                                                                      >  2:2805728/1‑100 (MQ=38)
               aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  1:3427644/100‑1 (MQ=255)
               aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTAc                                                                                     <  1:855608/100‑1 (MQ=38)
                gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                    >  2:2410885/1‑100 (MQ=255)
                gATATTTTATTCTATTACATGAAATAAACCAGCATACCCTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                    >  2:2431074/1‑100 (MQ=255)
                 atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAg                                                                                   <  2:194301/100‑1 (MQ=255)
                  tatTTTATTCTTTTACATGAAATAAACCCGCATACCATTTGAGCGATCGCACGCTGTTCCCCCATACTGTCAAGCCGGAGACATG‑CTCCGGCTTTACAg                                                                                   >  1:1696375/1‑99 (MQ=255)
                  tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGt                                                                                  <  1:1232435/100‑1 (MQ=255)
                    ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAt                                                                                <  2:3605879/100‑1 (MQ=37)
                      ttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTa                                                                              >  2:3852622/1‑100 (MQ=255)
                       tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                             <  1:2424601/100‑1 (MQ=255)
                         ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                           <  2:1345611/100‑1 (MQ=255)
                            ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        <  2:1481613/100‑1 (MQ=255)
                            ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        >  1:485804/1‑100 (MQ=255)
                               tACATGAAATAAACCAGCATACCATTTGAGCGACCGCCCGCTGTTCCCCCATAATGTAAAGCCGAACACATG‑CTCCGGCTTTACAGTAGTACTGAGAAt                                                                      <  1:194562/99‑1 (MQ=255)
                                aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                    <  1:199881/100‑1 (MQ=255)
                                aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                    <  1:3349172/100‑1 (MQ=255)
                                aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                    >  1:1533565/1‑100 (MQ=255)
                                 cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACTTG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   <  1:3744673/100‑1 (MQ=255)
                                   tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1047146/100‑1 (MQ=255)
                                   tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 >  1:2049616/1‑100 (MQ=255)
                                   tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1972837/100‑1 (MQ=255)
                                    gAAATAAAGCAGCATACCATTTTAGCGATCGCCCGCTGTACCCCCATACTGTAACGCCGGAGCCCTG‑CTCCCGCTTTACCCACTTACTGGGAATACAccc                                                                >  2:3685915/1‑100 (MQ=255)
                                    gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCGGGAGACATG‑CTCCGGCTTTACAATATTACGGAGAATACAccc                                                                >  2:329134/1‑100 (MQ=255)
                                    gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                >  2:1033094/1‑100 (MQ=255)
                                     aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               <  2:2538671/100‑1 (MQ=255)
                                      aagataCCATCATACCATTTTACCGGTCTCCCGCTTGTCCCCCATACTGTAAATCCGGAGACATG‑CTCCGTCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:2758075/92‑1 (MQ=255)
                                      aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               <  1:3790287/99‑1 (MQ=255)
                                      aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                              <  2:1915798/100‑1 (MQ=255)
                                      aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                              <  2:3672166/100‑1 (MQ=255)
                                        tAAAACAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGACGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                            <  2:3759105/100‑1 (MQ=255)
                                             cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                       >  1:1259970/1‑100 (MQ=255)
                                                 aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                   >  1:2015341/1‑100 (MQ=255)
                                                  tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  >  2:1466939/1‑100 (MQ=255)
                                                   aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                 <  1:2727219/100‑1 (MQ=255)
                                                    ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                <  1:2114054/100‑1 (MQ=255)
                                                         tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  2:3791397/100‑1 (MQ=255)
                                                          gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                          >  1:2102202/1‑100 (MQ=255)
                                                              gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATACTCGAAGGTTCAaaat                                      >  1:1848874/1‑100 (MQ=255)
                                                              gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                      >  1:3224249/1‑100 (MQ=255)
                                                               aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                     <  1:1135240/100‑1 (MQ=255)
                                                                tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                    <  1:2341001/100‑1 (MQ=255)
                                                                  gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  1:1886679/1‑100 (MQ=255)
                                                                   cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 >  2:648446/1‑100 (MQ=255)
                                                                      gCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtt                              <  1:3636382/100‑1 (MQ=255)
                                                                         gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           <  1:681309/100‑1 (MQ=255)
                                                                         gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           >  1:587628/1‑100 (MQ=255)
                                                                         gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           >  1:2472389/1‑100 (MQ=255)
                                                                          ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          <  1:3567892/100‑1 (MQ=255)
                                                                            cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                         <  2:3419578/99‑1 (MQ=255)
                                                                            cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                        >  2:2331/1‑100 (MQ=255)
                                                                             ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       <  2:52253/100‑1 (MQ=255)
                                                                             ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       >  1:2535892/1‑100 (MQ=255)
                                                                              cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  1:122853/1‑100 (MQ=255)
                                                                               ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  2:2208075/99‑1 (MQ=255)
                                                                                 aTACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCa                   >  1:339918/1‑100 (MQ=37)
                                                                                   aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAATCATACTCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 <  2:2151546/100‑1 (MQ=255)
                                                                                   aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  <  2:2434021/99‑1 (MQ=255)
                                                                                   aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 <  1:1093096/100‑1 (MQ=255)
                                                                                   aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 <  1:3189455/100‑1 (MQ=255)
                                                                                    cTGTACAGCCGGAG‑CATGTCTCGGGCTTTACAGTATTAATGAGAGTACAATCATCCCGTATCGTCGAAGGTCCAAGATCGccgc                                >  2:3102084/1‑82 (MQ=2)
                                                                                    cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  2:762574/1‑100 (MQ=38)
                                                                                      gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGCATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTTATcccc                    >  2:3217121/1‑94 (MQ=255)
                                                                                      gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  2:1200551/1‑100 (MQ=38)
                                                                                      gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTCGCTCGCCAGCGAt              >  2:3636055/1‑100 (MQ=37)
                                                                                      gTAAAGCAGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  1:2063539/1‑100 (MQ=255)
                                                                                        aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGACGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:1095009/1‑100 (MQ=255)
                                                                                        aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:1129870/1‑100 (MQ=255)
                                                                                        aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:3878914/1‑100 (MQ=38)
                                                                                        aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAACGTCAGCGAt              >  2:304522/1‑98 (MQ=37)
                                                                                        aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCATTTTGATCGCCAGCGATTc            >  2:3199683/1‑100 (MQ=38)
                                                                                        aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAAAAAAAATATAGCCCCTTTTTGATCGCCAGCGATTc            >  1:987573/1‑100 (MQ=25)
                                                                                         aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           <  1:2611314/100‑1 (MQ=255)
                                                                                          aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt          <  2:3827187/100‑1 (MQ=255)
                                                                                           gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         <  1:493022/100‑1 (MQ=255)
                                                                                           gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc         >  1:1426014/1‑100 (MQ=255)
                                                                                             cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc       >  2:3306479/1‑100 (MQ=255)
                                                                                              ggAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAcc      <  1:3794542/100‑1 (MQ=255)
                                                                                                agaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGt    <  2:507460/100‑1 (MQ=255)
                                                                                                 gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg   <  2:2237688/100‑1 (MQ=255)
                                                                                                 gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg   >  1:1710994/1‑100 (MQ=255)
                                                                                                  aCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  1:3841734/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520402‑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: ‑