Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 19.2% 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.T19.2% 368.6 / 104.5 125intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (48/53);  new base T (11/13);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 128 total aligned reads displayed.

TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520406‑4520603
                                                                                                   |                                                                                                   
tgcttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGc                                                                                                    >  1:1149188/4‑100 (MQ=255)
    tttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCg                                                                                                <  1:1353540/100‑1 (MQ=255)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTACAGCCGGAGACATGCCCCGGc                                                                                              <  2:2438073/100‑1 (MQ=255)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                              >  1:2161/1‑100 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAGCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAGGCCGGAGCATGTCTCCGGCt                                                                                             <  2:1595532/100‑1 (MQ=38)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                             <  2:1391109/100‑1 (MQ=255)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCTGCTGTTCCCCGATACTGTAAAGCCGGAGGCATGCTCCGGCtt                                                                                            <  1:1757070/100‑1 (MQ=255)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                            <  1:64523/100‑1 (MQ=38)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                            <  2:783461/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCNGGCttt                                                                                           <  2:2289966/100‑1 (MQ=255)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  1:1917898/100‑1 (MQ=255)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  1:258874/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  2:2365579/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  2:2188065/100‑1 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  2:521018/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  1:981007/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTa                                                                                          <  2:714741/100‑1 (MQ=38)
           aGATATTTTATTCTTTTACATGAAATAACCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                         <  2:508987/100‑1 (MQ=38)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  2:2044582/100‑1 (MQ=255)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                         <  2:2463317/100‑1 (MQ=38)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGGTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  2:1912851/100‑1 (MQ=255)
            gctatTTTATCCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                        <  2:1578542/98‑1 (MQ=37)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                        <  2:1553824/100‑1 (MQ=255)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  2:676796/100‑1 (MQ=38)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  1:695243/100‑1 (MQ=38)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATATGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAg                                                                                       <  2:296398/100‑1 (MQ=255)
              tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGt                                                                                      >  2:1140048/1‑100 (MQ=255)
              tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGt                                                                                      <  1:1729476/100‑1 (MQ=38)
               atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                     <  1:1835183/100‑1 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             <  2:1872213/100‑1 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             <  1:1140053/100‑1 (MQ=255)
                        ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                            <  2:1690591/100‑1 (MQ=255)
                         tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGACaa                                                                           >  2:1549315/1‑100 (MQ=255)
                          ttACATGAAATAAACCATCATACCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                           <  1:2433923/99‑1 (MQ=255)
                          ttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAt                                                                          <  1:1572571/100‑1 (MQ=255)
                           tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCGGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAt                                                                          <  1:2156946/99‑1 (MQ=255)
                            aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAAATGCTCCGGCTTTACAGTATTACTGAGAATa                                                                         >  1:2384023/1‑99 (MQ=255)
                             cacgAAATAACCCAGCATACCATTAGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGGCGGAGACAGGCTCCGGCTTTACAGAATTACTGCGAATa                                                                         <  2:2369464/95‑1 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTGCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                       <  2:1834708/100‑1 (MQ=255)
                             cATGAAATAAACCAGAATACCATTTGAGCGATCGCCCGATGTACCCCCATACTGTAAAGCAGGAGACATGCTCGGGCTTTACAGTATTACTGAGTATaca                                                                       >  2:1772393/1‑100 (MQ=255)
                              aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                      >  2:2461003/1‑100 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     <  2:1388985/100‑1 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     <  1:1971807/100‑1 (MQ=255)
                                   aTAAACCAGCATACCATTTGAGCGATCGCCCGATGTTCCCCCATACTCTAAAGCCGGCGAAATGATCCAACGTTACACTACTACTGAGAATACACACATc                                                                 >  2:2203943/1‑100 (MQ=255)
                                    tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                >  1:704607/1‑100 (MQ=255)
                                    tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                >  2:308824/1‑100 (MQ=255)
                                    tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                <  1:1390157/100‑1 (MQ=255)
                                      aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                              >  1:370295/1‑100 (MQ=255)
                                         cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                           <  1:2326987/100‑1 (MQ=255)
                                         cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                           >  1:1197398/1‑100 (MQ=255)
                                              tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                      <  1:1512959/100‑1 (MQ=255)
                                              tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                      >  1:1385889/1‑100 (MQ=255)
                                               aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                     >  1:2048762/1‑100 (MQ=255)
                                                 cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                   <  2:254392/100‑1 (MQ=255)
                                                 cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                   <  1:2106255/100‑1 (MQ=255)
                                                   tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGt                                                 <  1:1690006/100‑1 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               >  1:2289632/1‑100 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               <  2:2555194/100‑1 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               <  1:160377/100‑1 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               >  1:1182391/1‑100 (MQ=255)
                                                      gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                              >  1:1624549/1‑100 (MQ=255)
                                                      gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                              >  2:1751239/1‑100 (MQ=255)
                                                         cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                           <  1:1675638/100‑1 (MQ=255)
                                                          gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                          >  1:113450/1‑100 (MQ=255)
                                                             cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                       >  1:1401878/1‑100 (MQ=255)
                                                               cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  2:1919171/1‑100 (MQ=255)
                                                               cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     <  1:10697/100‑1 (MQ=255)
                                                                ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  2:2305259/1‑99 (MQ=255)
                                                                      ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              >  1:2316793/1‑100 (MQ=255)
                                                                      ctcagtaaTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              >  1:1092801/8‑100 (MQ=21)
                                                                       tCCCGCATACTGTAAAGCCGGATACACTATCCGGCTTTACATTATTACTGAGAATACACCCATCAAGCATAGTCGAAGGTTCAATATCGCCCCttttt                               <  2:1987693/98‑1 (MQ=255)
                                                                       tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  1:927589/100‑1 (MQ=255)
                                                                        cccccATACTGTAAAACCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            >  1:2104530/1‑100 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           >  1:733683/1‑100 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          <  2:370291/100‑1 (MQ=255)
                                                                           ccATACTGTAAAGACGGAGACATGATCCGGCTTTACAGTATTACTGAGAATAAAACCATCAAGCATAATAGAAGGCTCAATATTGTCACATTTTg                              >  2:2089208/1‑95 (MQ=255)
                                                                              tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  2:2000785/1‑100 (MQ=255)
                                                                              tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  2:1718110/1‑100 (MQ=255)
                                                                              tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATAAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  1:898399/1‑100 (MQ=255)
                                                                                cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  1:798345/1‑100 (MQ=255)
                                                                                 tGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                   >  1:906667/1‑100 (MQ=38)
                                                                                  gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  <  2:848686/100‑1 (MQ=255)
                                                                                  gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  1:420335/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  2:1629187/1‑98 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:1072552/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:1072778/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:984328/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:127143/1‑100 (MQ=38)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               <  1:2517545/100‑1 (MQ=255)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  1:1360650/1‑100 (MQ=255)
                                                                                     aaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  1:1046448/1‑100 (MQ=38)
                                                                                      aGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt              >  2:1621369/1‑100 (MQ=255)
                                                                                         cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc           <  1:144866/100‑1 (MQ=255)
                                                                                             gaCATGCTCCGCCTTTACAGTATTACTGAGAATACGCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg       <  2:734973/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  2:2302180/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAGCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  1:960384/100‑1 (MQ=255)
                                                                                               cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  2:1072257/1‑100 (MQ=255)
                                                                                               cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     <  1:478703/100‑1 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  2:658770/1‑100 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  2:683114/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: