Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 16.7% intergenic (‑180/+167) 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,5041.T16.7% 320.8 / 69.8 102intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (47/38);  new base T (7/10);  total (54/48)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 106 total aligned reads displayed.

TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520406‑4520603
                                                                                                   |                                                                                                   
tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGc                                                                                                    >  1:2164836/1‑100 (MQ=255)
     ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCgg                                                                                               <  2:1272504/100‑1 (MQ=38)
     ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCgg                                                                                               >  2:1865405/1‑100 (MQ=255)
     ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCgg                                                                                               <  1:796386/100‑1 (MQ=39)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                              <  2:412826/100‑1 (MQ=255)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGc                                                                                              <  2:2338252/100‑1 (MQ=39)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGc                                                                                              <  2:2686090/100‑1 (MQ=39)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                             <  2:208636/100‑1 (MQ=255)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  2:1716603/100‑1 (MQ=38)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTa                                                                                          <  2:2502630/100‑1 (MQ=38)
           aGATATTTTATTCTTTTACATGAAATAACCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  1:1122089/100‑1 (MQ=255)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCTAGCTGTTCCCCCATACTGTAAAACCGGAGACATGCTCCGGCTTTa                                                                                          >  1:1352346/1‑99 (MQ=255)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  2:1401147/100‑1 (MQ=255)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                         <  2:233231/100‑1 (MQ=38)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTTCCCCATACTACAAAGCCGGAGACAGGCTCCgg                                                                                               >  1:2523651/1‑93 (MQ=255)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                        <  2:125419/100‑1 (MQ=255)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                        <  1:691585/100‑1 (MQ=255)
            gATATTTTATTCTCTTACATGAAATAAACCAACATACCATTTAAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                        >  2:1967624/1‑100 (MQ=255)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  2:2451909/100‑1 (MQ=38)
                ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAt                                                                                    <  2:253695/100‑1 (MQ=255)
                   tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTAATCCCACCTACTGAAAAGCCGGAGAAATACTACGGCAAAACAGTATTAc                                                                                 >  2:1696358/1‑100 (MQ=255)
                    aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACt                                                                                <  1:2549540/100‑1 (MQ=255)
                    aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTAtgggg                                                                                <  2:457403/100‑5 (MQ=25)
                     ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                               <  2:2904495/100‑1 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAACGTCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  1:1433283/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             <  2:2123667/100‑1 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  1:1631650/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGAGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACt                                                                                >  2:2703473/1‑97 (MQ=255)
                       cTTTTACATCATATACACCAGCATTCCATTTGATCGCTCGCCCTCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTGACAGTATTa                                                                                  <  1:2165717/95‑1 (MQ=37)
                         tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                           >  1:154420/1‑100 (MQ=255)
                         tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                           <  1:90219/100‑1 (MQ=255)
                          ttACATGAAATAAACCAGCATACCATTTGAGCGATCGACCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAt                                                                          >  1:215349/1‑100 (MQ=255)
                           tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATa                                                                         >  2:2681518/1‑100 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                       >  1:384953/1‑100 (MQ=255)
                              aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                      >  2:1277405/1‑100 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     <  2:1107846/100‑1 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     >  2:508579/1‑100 (MQ=255)
                                  aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                  >  1:615331/1‑100 (MQ=255)
                                   aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGACTTTACAGTATTACTGAGAATACAccc                                                                    <  1:1061428/97‑1 (MQ=255)
                                    tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                <  2:2241549/100‑1 (MQ=255)
                                     aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                               >  2:1016486/1‑100 (MQ=255)
                                      aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                              <  1:2803082/100‑1 (MQ=255)
                                       aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                             <  1:1328258/100‑1 (MQ=255)
                                            cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTCCACTATTACTGAGAATACACCCATCAAGCATAAt                                                        >  1:2751691/1‑100 (MQ=255)
                                              tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                      <  2:2310356/100‑1 (MQ=255)
                                              tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                      >  2:644728/1‑100 (MQ=255)
                                               aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                     >  2:2579184/1‑100 (MQ=255)
                                                  aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                  >  2:1534204/1‑100 (MQ=255)
                                                   tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATCCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                  <  1:2597241/99‑1 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               >  2:1821532/1‑100 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAGACACCCATCAAGCATAATCGAAGGTTc                                               <  1:1387014/100‑1 (MQ=255)
                                                        gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            >  1:2260212/1‑100 (MQ=255)
                                                         cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                           <  2:575124/100‑1 (MQ=255)
                                                            tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                        <  2:1847787/100‑1 (MQ=255)
                                                            tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                        <  1:828589/100‑1 (MQ=255)
                                                             cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                       >  1:2701225/1‑100 (MQ=255)
                                                             cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAGCAAGCATAATCGAAGGTTCAATATCg                                        <  2:1371499/99‑1 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      >  2:2331138/1‑100 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      >  1:416775/1‑100 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      >  1:923483/1‑100 (MQ=255)
                                                               cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  2:2845061/1‑100 (MQ=255)
                                                                ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                       >  1:1628899/1‑97 (MQ=255)
                                                                    tGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                                >  2:295345/1‑100 (MQ=255)
                                                                     gTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                               >  1:844156/1‑100 (MQ=255)
                                                                      ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              <  2:192894/100‑1 (MQ=255)
                                                                       tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  2:416777/100‑1 (MQ=255)
                                                                       tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  2:1158828/100‑1 (MQ=255)
                                                                        cccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            >  2:1839513/1‑100 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            <  1:1897475/99‑1 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACATGCTCCGCCTTTACAGTATTACTGAGAATACACCCATCAAGCATACTCGAAGGTTCAATATCGCCCCTTTTTGagc                           >  1:1076139/1‑98 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          <  2:1400889/100‑1 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          >  2:439736/1‑100 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          >  2:760535/1‑100 (MQ=255)
                                                                            cATACTGTAAAGCCGGAGACATGCTCCGACTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                        >  1:1887567/1‑100 (MQ=255)
                                                                             aTACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                        <  1:930286/99‑1 (MQ=255)
                                                                              tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGATCAATATCGCCCCTTTTTGATCGCCAg                      >  1:468265/1‑100 (MQ=255)
                                                                               aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                     <  1:1487489/100‑1 (MQ=255)
                                                                                cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  1:2330951/1‑100 (MQ=255)
                                                                                cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATTGCCTGc                     >  2:1998221/1‑99 (MQ=37)
                                                                                cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                     >  1:2735277/1‑99 (MQ=38)
                                                                                  gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGGGa                   >  2:2652351/1‑99 (MQ=255)
                                                                                  gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTCTTGATCGCAAACGAt                  >  1:527417/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:620269/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:281693/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:971474/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:934567/1‑100 (MQ=38)
                                                                                     aaGCCGGAGACATGTTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               <  1:997888/100‑1 (MQ=255)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  1:1348299/1‑100 (MQ=255)
                                                                                     aaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  2:2687627/1‑100 (MQ=38)
                                                                                        ccGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa            >  2:1558549/1‑100 (MQ=255)
                                                                                        ccGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa            <  2:1190495/100‑1 (MQ=255)
                                                                                             gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg       >  1:506142/1‑100 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  1:1664166/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  2:743438/100‑1 (MQ=255)
                                                                                               cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGTTCGCCAGCGATTCGTCACCGTGc      <  1:1844920/99‑1 (MQ=255)
                                                                                               cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  2:1360381/1‑100 (MQ=255)
                                                                                                 tGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGa   >  1:2187493/1‑100 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGGATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  2:2848797/1‑100 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  1:1457847/1‑100 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  1:2482878/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520406‑4520603

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: