Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 19.8% 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.8% 324.8 / 83.8 121intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (50/46);  new base . (11/13);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 127 total aligned reads displayed.

TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520401‑4520599
                                                                                                   |                                                                                                    
 ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                     >  2:2395207/1‑100 (MQ=255)
  gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    >  2:1353533/1‑100 (MQ=255)
   cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  1:1491185/1‑100 (MQ=255)
     tgcttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGc                                                                                                >  1:1149188/4‑100 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTACAGCCGGAGACATGCCCCGGc                                                                                          <  2:2438073/100‑1 (MQ=255)
           tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                          >  1:2161/1‑100 (MQ=255)
            ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                         <  2:1391109/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCTGCTGTTCCCCGATACTGTAAAGCCGGAGGCATGCTCCGGCtt                                                                                        <  1:1757070/100‑1 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                        >  2:2530991/1‑100 (MQ=255)
             taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                        <  2:783461/100‑1 (MQ=38)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCNGGCttt                                                                                       <  2:2289966/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                       <  1:1917898/100‑1 (MQ=255)
              aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                       <  2:2188065/100‑1 (MQ=255)
               cAGATATTTTCTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATTCTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                      <  2:1713097/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)
                aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                     <  2:1803789/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)
                   tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGt                                                                                  <  1:1729476/100‑1 (MQ=38)
                    atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                 <  1:1835183/100‑1 (MQ=255)
                            cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                         <  1:1140053/100‑1 (MQ=255)
                             ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                        <  2:1690591/100‑1 (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)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  1:1971807/100‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 <  2:1388985/100‑1 (MQ=255)
                                    tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                 >  1:2273044/1‑100 (MQ=255)
                                      aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGGTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                               >  2:1395718/1‑100 (MQ=255)
                                         tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                            >  2:308824/1‑100 (MQ=255)
                                         tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                            >  1:704607/1‑100 (MQ=255)
                                         tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                            <  1:1390157/100‑1 (MQ=255)
                                           aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                          >  1:370295/1‑100 (MQ=255)
                                             ccAGCATACCATTTGAGCGACCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                        >  2:1585666/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                                               >  1:366836/1‑100 (MQ=255)
                                                      cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                               <  1:2106255/100‑1 (MQ=255)
                                                      cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                               <  2:254392/100‑1 (MQ=255)
                                                        tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGt                                             <  1:1690006/100‑1 (MQ=255)
                                                         ttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGtt                                            <  1:927449/100‑1 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           >  1:2289632/1‑100 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           <  1:160377/100‑1 (MQ=255)
                                                          tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                           >  1:1182391/1‑100 (MQ=255)
                                                           gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                          >  2:1751239/1‑100 (MQ=255)
                                                           gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                          >  1:1624549/1‑100 (MQ=255)
                                                              cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                       <  1:1675638/100‑1 (MQ=255)
                                                               gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                      >  1:113450/1‑100 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   <  2:2384009/100‑1 (MQ=255)
                                                                  cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                   >  1:1401878/1‑100 (MQ=255)
                                                                   gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                  >  1:1287917/1‑100 (MQ=255)
                                                                    nccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 >  1:555302/2‑100 (MQ=255)
                                                                    cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 <  1:10697/100‑1 (MQ=255)
                                                                    cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 >  2:1919171/1‑100 (MQ=255)
                                                                     ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                 >  2:2305259/1‑99 (MQ=255)
                                                                        attctcagtaaTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGGATACAGCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                             >  2:1063777/11‑100 (MQ=18)
                                                                           ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          >  1:2316793/1‑100 (MQ=255)
                                                                           ctcagtaaTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                          >  1:1092801/8‑100 (MQ=21)
                                                                            tCCCGCATACTGTAAAGCCGGATACACTATCCGGCTTTACATTATTACTGAGAATACACCCATCAAGCATAGTCGAAGGTTCAATATCGCCCCttttt                           <  2:1987693/98‑1 (MQ=255)
                                                                             cccccATACTGTAAAACCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                        >  1:2104530/1‑100 (MQ=255)
                                                                              ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                       >  1:733683/1‑100 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  2:370291/100‑1 (MQ=255)
                                                                               cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                      <  2:2048751/100‑1 (MQ=255)
                                                                                ccATACTGTAAAGACGGAGACATGATCCGGCTTTACAGTATTACTGAGAATAAAACCATCAAGCATAATAGAAGGCTCAATATTGTCACATTTTg                          >  2:2089208/1‑95 (MQ=255)
                                                                                   tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  2:2000785/1‑100 (MQ=255)
                                                                                   tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATAAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                  >  1:898399/1‑100 (MQ=255)
                                                                                    aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 <  2:347240/100‑1 (MQ=255)
                                                                                     cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                 >  2:2296972/1‑99 (MQ=38)
                                                                                      tGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa               >  2:57325/1‑100 (MQ=255)
                                                                                      tGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa               >  1:906667/1‑100 (MQ=38)
                                                                                       gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  1:420335/1‑100 (MQ=255)
                                                                                       gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              <  2:848686/100‑1 (MQ=255)
                                                                                         aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:759279/1‑100 (MQ=255)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt              >  2:1629187/1‑98 (MQ=38)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:984328/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:127143/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:1072778/1‑100 (MQ=38)
                                                                                         aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc            >  1:1072552/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)
                                                                                                   |                                                                                                    
TGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGC  >  NC_000913/4520401‑4520599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: