Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 16.3% 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.3% 543.2 / 119.0 172intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (71/73);  new base T (14/14);  total (85/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-02

Only 100 of 175 total aligned reads displayed.

TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520406‑4520603
                                                                                                   |                                                                                                   
tGATTTTATACAGACATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATCCTGTAAAGCCGGAGACATGc                                                                                                    <  1:1223156/100‑1 (MQ=255)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                              <  2:1720084/100‑1 (MQ=255)
      tataCAGATATTGTATTCTTTTACATCAAATAAACCAGCATACCATATGGGCGATCGCACGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCgg                                                                                               <  2:1133464/99‑1 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                              >  2:1661392/1‑99 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                             <  1:627193/100‑1 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                             <  1:964957/100‑1 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                             <  2:1156136/100‑1 (MQ=255)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                            <  1:329207/100‑1 (MQ=38)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                            <  1:2493946/100‑1 (MQ=38)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtn                                                                                            <  1:1061355/100‑2 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  1:2953561/100‑1 (MQ=255)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  2:386318/100‑1 (MQ=255)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  1:1846515/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  2:892580/100‑1 (MQ=38)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  1:1550404/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  2:715298/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  2:1934819/1‑100 (MQ=255)
          cAGATATTTTATTCTTGTACATGAAATAAACCAGCATACCATTTGAGCGGTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  2:736501/99‑1 (MQ=37)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                         <  1:1674778/100‑1 (MQ=38)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          <  2:1867260/98‑1 (MQ=255)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                        <  1:2607887/100‑1 (MQ=38)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  1:3789361/100‑1 (MQ=37)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAg                                                                                       <  1:3671015/100‑1 (MQ=255)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAg                                                                                       <  2:1768013/100‑1 (MQ=255)
              tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTTATCTATCGCCCGCTGTTCCACCATACTTTAAAGCCGGAGACATGCTCCGTCTTTATAGt                                                                                      >  2:1898521/1‑100 (MQ=255)
                   tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATACTGTAAAGCCGGACACATGCTCCGGCTTCACAGTATTAc                                                                                 >  2:3134269/1‑100 (MQ=255)
                   tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                                 <  2:1767304/100‑1 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  2:364600/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  2:2942974/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  1:2325947/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             <  2:2737115/100‑1 (MQ=255)
                            aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATac                                                                        >  1:1904236/1‑100 (MQ=255)
                            aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATac                                                                        <  2:195472/100‑1 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                       >  2:3449753/1‑100 (MQ=255)
                              aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                      >  2:618447/1‑100 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     >  1:2699889/1‑100 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     >  1:68018/1‑100 (MQ=255)
                               tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     >  2:2015036/1‑100 (MQ=255)
                                gAAATAAACCAGCATACCATTTGAGCGCTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                    <  2:3033737/100‑1 (MQ=255)
                                gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                    >  1:983706/1‑100 (MQ=255)
                                  aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                  >  1:2753806/1‑100 (MQ=255)
                                     aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                               >  1:2818151/1‑100 (MQ=255)
                                      aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                              <  1:1840749/100‑1 (MQ=255)
                                        ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                            >  1:2659097/1‑100 (MQ=255)
                                             aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                       >  2:2557368/1‑100 (MQ=255)
                                             aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                       >  1:1851681/1‑100 (MQ=255)
                                             aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATCCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                       <  2:3495659/100‑1 (MQ=255)
                                              tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                      <  1:3672876/100‑1 (MQ=255)
                                                   tatGAGCGAACGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                   <  2:2252583/96‑1 (MQ=255)
                                                    ttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGt                                                 >  2:431599/1‑99 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               <  2:3275811/100‑1 (MQ=255)
                                                     tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               >  1:493629/1‑100 (MQ=255)
                                                      gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                              <  2:2785288/100‑1 (MQ=255)
                                                      gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAGGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAGACTCGAAGGTTCa                                              >  2:1220759/1‑100 (MQ=255)
                                                        gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            <  1:1866315/100‑1 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatagatcg                                      >  1:1964610/1‑95 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      <  2:1424713/100‑1 (MQ=255)
                                                              gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                       <  2:3664376/99‑1 (MQ=255)
                                                               cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  1:2314986/1‑100 (MQ=255)
                                                                   cTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGCATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtt                                  >  2:1122933/1‑99 (MQ=255)
                                                                   cTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                 <  1:107437/100‑1 (MQ=255)
                                                                     gTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                               <  1:3694000/100‑1 (MQ=255)
                                                                     gTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATCACTGAGAATACACTCATCAAGCATAATCGAAGGTTCAATATTGCCCCtttt                                >  1:1446255/1‑99 (MQ=255)
                                                                      ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              <  1:1898494/100‑1 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           <  2:2749667/99‑1 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          >  2:538935/1‑100 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          >  1:350918/1‑100 (MQ=255)
                                                                            aaTACTGTAAAGCCGGAGACATGCACCGGCTAGATAGTATTACTGAGTATACCCCCGTCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  2:2991853/94‑1 (MQ=255)
                                                                             aTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCa                       >  1:3360693/1‑100 (MQ=37)
                                                                              tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAATACTGAGAATACATCCGTCAAGCATAATTTAAGGTTCAAAATCGCCCCTATCTGATCGCCAg                      >  2:2388404/1‑100 (MQ=37)
                                                                              tACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  1:2233189/1‑100 (MQ=37)
                                                                               aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAgn                     <  2:1025704/100‑2 (MQ=255)
                                                                               aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                     <  1:112292/100‑1 (MQ=255)
                                                                               aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                     <  2:493624/100‑1 (MQ=255)
                                                                                 tttAAAGCCGTAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCTTCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    <  2:1251683/97‑1 (MQ=255)
                                                                                 tGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                   >  1:1479051/1‑100 (MQ=255)
                                                                                  gTAAAGCCGGATACATGCTCCGGCTTTACAGTATTACTGAGAATACACACATCAAGCATAAACGAAGGTTCAATATCTCCCCTTTTAGACCCCCAGCGAt                  >  2:75571/1‑100 (MQ=255)
                                                                                  gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  1:3729046/1‑100 (MQ=255)
                                                                                  gTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  1:78367/1‑100 (MQ=38)
                                                                                   tAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATCATCGAAGGTTCAATATCGCCCCtttttttt                            >  2:1441165/1‑86 (MQ=25)
                                                                                    aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:3198659/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:2781379/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:3725545/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:626497/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:1633695/1‑100 (MQ=38)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  1:399896/1‑100 (MQ=255)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               <  2:3762643/100‑1 (MQ=255)
                                                                                     aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGCTCAATATCGCCCCTTTTTGATCGCCAGAGATTCg               >  1:1589276/1‑100 (MQ=255)
                                                                                        ccGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa            <  1:1606562/100‑1 (MQ=255)
                                                                                         cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc           <  1:2122370/100‑1 (MQ=255)
                                                                                          ggAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAcc          >  2:2970152/1‑100 (MQ=255)
                                                                                          ggAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAcc          >  2:2657463/1‑100 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg       >  1:2406114/1‑99 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  2:3674360/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  2:3197793/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  1:2803773/100‑1 (MQ=255)
                                                                                              aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  1:32143/100‑1 (MQ=255)
                                                                                               cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  2:2200841/1‑100 (MQ=255)
                                                                                                 tGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGATCAATATAGCCCCTTTTTGATCGCAAGCGATTCGTCACCGTGCAGa   >  2:3371656/1‑100 (MQ=255)
                                                                                                  gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  1:571189/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: