Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 16.3% 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.3% 293.5 / 55.0 104intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (47/40);  new base . (7/10);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

Only 100 of 107 total aligned reads displayed.

GGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520402‑4520599
                                                                                                  |                                                                                                    
ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  1:1226403/100‑1 (MQ=255)
ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     <  2:1943529/100‑1 (MQ=255)
 gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    <  2:2634312/100‑1 (MQ=255)
 gCATAATGTTATACAAATATTTTTGTCTTTTACACGAAATAATCCATCAGAAAATTTGAGCGATCGCCCGTTGTTCCCCCCTCCTGTAAAGCCGGAGACa                                                                                                    <  1:1390145/100‑1 (MQ=255)
  cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   <  2:1335522/100‑1 (MQ=255)
  cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  1:1119860/1‑100 (MQ=255)
  cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  2:1037381/1‑100 (MQ=255)
    tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑c                                                                                                >  1:2164836/1‑100 (MQ=255)
         ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCgg                                                                                           <  2:1272504/100‑1 (MQ=38)
         ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCgg                                                                                           >  2:1865405/1‑100 (MQ=255)
         ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCgg                                                                                           <  1:796386/100‑1 (MQ=39)
          tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                          <  2:412826/100‑1 (MQ=255)
          tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                          <  2:2686090/100‑1 (MQ=39)
          tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                          <  2:2338252/100‑1 (MQ=39)
           ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCt                                                                                         <  2:208636/100‑1 (MQ=255)
             aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                       <  2:1716603/100‑1 (MQ=38)
              cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTa                                                                                      <  2:2502630/100‑1 (MQ=38)
               aGATATTTTATTCTTTTACATGAAATAACCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  1:1122089/100‑1 (MQ=255)
               aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCTAGCTGTTCCCCCATACTGTAAAACCGGAGACATG‑CTCCGGCTTTa                                                                                      >  1:1352346/1‑99 (MQ=255)
               aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                     <  2:1401147/100‑1 (MQ=255)
               aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTAc                                                                                     <  2:233231/100‑1 (MQ=38)
                gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTTCCCCATACTACAAAGCCGGAGACAGG‑CTCCgg                                                                                           >  1:2523651/1‑93 (MQ=255)
                gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                    <  2:125419/100‑1 (MQ=255)
                gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                    <  1:691585/100‑1 (MQ=255)
                gATATTTTATTCTCTTACATGAAATAAACCAACATACCATTTAAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                    >  2:1967624/1‑100 (MQ=255)
                 atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                   <  2:2451909/100‑1 (MQ=38)
                    ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                                <  2:253695/100‑1 (MQ=255)
                       tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTAATCCCACCTACTGAAAAGCCGGAGAAATA‑CTACGGCAAAACAGTATTAc                                                                             >  2:1696358/1‑100 (MQ=255)
                       tATTCGTTTACATGAGATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCTCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                             <  1:202819/100‑1 (MQ=255)
                        aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACt                                                                            <  1:2549540/100‑1 (MQ=255)
                        aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAtgggg                                                                            <  2:457403/100‑5 (MQ=25)
                         ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                           <  2:2904495/100‑1 (MQ=255)
                           cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAACGTCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         >  1:1433283/1‑100 (MQ=255)
                           cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         >  1:1631650/1‑100 (MQ=255)
                           cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                         <  2:2123667/100‑1 (MQ=255)
                           cTTTTACATGAAATAAACCAGCATACCATTTGAGAGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACt                                                                            >  2:2703473/1‑97 (MQ=255)
                           cTTTTACATCATATACACCAGCATTCCATTTGATCGCTCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTGACAGTATTa                                                                              <  1:2165717/95‑1 (MQ=37)
                            ttttACATGAAATAAACCAGCATACCATTTGTGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                        <  2:2811701/100‑1 (MQ=255)
                             tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                       >  1:154420/1‑100 (MQ=255)
                             tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                       <  1:90219/100‑1 (MQ=255)
                              ttACATGAAATAAACCAGCATACCATTTGAGCGATCGACCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAAt                                                                      >  1:215349/1‑100 (MQ=255)
                               tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATa                                                                     >  2:2681518/1‑100 (MQ=255)
                                 cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                   >  1:384953/1‑100 (MQ=255)
                                  aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATacac                                                                  >  2:1277405/1‑100 (MQ=255)
                                   tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 >  2:508579/1‑100 (MQ=255)
                                   tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1107846/100‑1 (MQ=255)
                                      aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                              >  1:615331/1‑100 (MQ=255)
                                       aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGACTTTACAGTATTACTGAGAATACAccc                                                                <  1:1061428/97‑1 (MQ=255)
                                        tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                            <  2:2241549/100‑1 (MQ=255)
                                         aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                           >  2:1016486/1‑100 (MQ=255)
                                           aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                         <  1:1328258/100‑1 (MQ=255)
                                                cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTCCACTATTACTGAGAATACACCCATCAAGCATAAt                                                    >  1:2751691/1‑100 (MQ=255)
                                                  tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                  >  2:644728/1‑100 (MQ=255)
                                                   aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                 >  2:2579184/1‑100 (MQ=255)
                                                      aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                              >  2:1534204/1‑100 (MQ=255)
                                                       tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATC‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                              <  1:2597241/99‑1 (MQ=255)
                                                         tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           >  2:1821532/1‑100 (MQ=255)
                                                         tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAAGACACCCATCAAGCATAATCGAAGGTTc                                           <  1:1387014/100‑1 (MQ=255)
                                                          gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                          >  2:1210766/1‑100 (MQ=255)
                                                            gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                        >  1:2260212/1‑100 (MQ=255)
                                                             cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                       <  2:575124/100‑1 (MQ=255)
                                                               aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                      >  1:499611/1‑99 (MQ=255)
                                                                tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                    <  2:1847787/100‑1 (MQ=255)
                                                                tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                    <  1:828589/100‑1 (MQ=255)
                                                                 cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   >  1:2701225/1‑100 (MQ=255)
                                                                 cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAGCAAGCATAATCGAAGGTTCAATATCg                                    <  2:1371499/99‑1 (MQ=255)
                                                                  gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  1:923483/1‑100 (MQ=255)
                                                                  gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  1:416775/1‑100 (MQ=255)
                                                                  gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  2:2331138/1‑100 (MQ=255)
                                                                   cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 >  2:2845061/1‑100 (MQ=255)
                                                                    ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   >  1:1628899/1‑97 (MQ=255)
                                                                        tGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                            >  2:295345/1‑100 (MQ=255)
                                                                         gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                           >  1:844156/1‑100 (MQ=255)
                                                                          ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          <  2:192894/100‑1 (MQ=255)
                                                                           tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                         <  2:416777/100‑1 (MQ=255)
                                                                           tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                         <  2:1158828/100‑1 (MQ=255)
                                                                            cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                        >  2:1839513/1‑100 (MQ=255)
                                                                             ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                        <  1:1897475/99‑1 (MQ=255)
                                                                             ccccATACTGTAAAGCCGGAGACATG‑CTCCGCCTTTACAGTATTACTGAGAATACACCCATCAAGCATACTCGAAGGTTCAATATCGCCCCTTTTTGagc                       >  1:1076139/1‑98 (MQ=255)
                                                                              cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  2:760535/1‑100 (MQ=255)
                                                                              cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  2:1400889/100‑1 (MQ=255)
                                                                              cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      >  2:439736/1‑100 (MQ=255)
                                                                                cATACTGTAAAGCCGGAGACATG‑CTCCGACTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    >  1:1887567/1‑100 (MQ=255)
                                                                                 aTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                    <  1:930286/99‑1 (MQ=255)
                                                                                  tACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGATCAATATCGCCCCTTTTTGATCGCCAg                  >  1:468265/1‑100 (MQ=255)
                                                                                   aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 <  1:1487489/100‑1 (MQ=255)
                                                                                    cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                >  1:2330951/1‑100 (MQ=255)
                                                                                    cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATTGCCTGc                 >  2:1998221/1‑99 (MQ=37)
                                                                                    cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 >  1:2735277/1‑99 (MQ=38)
                                                                                      gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGGGa               >  2:2652351/1‑99 (MQ=255)
                                                                                      gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTCTTGATCGCAAACGAt              >  1:527417/1‑100 (MQ=255)
                                                                                       tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAAAATCGCAACTTTTTCATCCACACCGAt              >  2:1994946/1‑99 (MQ=21)
                                                                                        aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:620269/1‑100 (MQ=255)
                                                                                        aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  2:934567/1‑100 (MQ=38)
                                                                                         aaGCCGGAGACATG‑TTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           <  1:997888/100‑1 (MQ=255)
                                                                                         aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           >  1:1348299/1‑100 (MQ=255)
                                                                                         aaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg           >  2:2687627/1‑100 (MQ=38)
                                                                                            ccGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa        <  2:1190495/100‑1 (MQ=255)
                                                                                                 gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg   >  1:506142/1‑100 (MQ=255)
                                                                                                  aCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc  <  2:743438/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: ‑