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

GCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGT  >  NC_000913/4520355‑4520648
                                                                                                                                                 |                                                                                                                                                     
gCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATg                                                                                                                                                   <  1:1733392/150‑1 (MQ=255)
   cATCTTAACCCTGTTTACCCGGAGTGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGATCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTc                                                                                                                                               >  2:1170126/1‑150 (MQ=255)
   cATCTTAACCCTGTTTACCCGGAGAGTTGATTCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTATTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTc                                                                                                                                               <  2:524943/150‑1 (MQ=255)
         aaCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGGTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                                                                         <  1:682712/150‑1 (MQ=255)
          aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                                                                        <  1:899776/150‑1 (MQ=255)
          aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATTCTCCGGCttt                                                                                                                                        >  1:1722934/1‑150 (MQ=255)
          aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                                                                        <  2:1451029/150‑1 (MQ=255)
          aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                                                                        <  1:463847/150‑1 (MQ=255)
            ccTGTTTACCCGGATAGTTGATGCAGATTGACAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                                                                      <  1:1278792/150‑1 (MQ=255)
             cTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                                                                     <  2:590030/150‑1 (MQ=255)
              tGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAg                                                                                                                                    <  1:415376/150‑1 (MQ=255)
               gTTTTCCCGGAGAGGTGATGCAGATTGCCAGTGGCATGATTTTATTCAGATATTTTATTCTTTTACATGAAATATACCAGCATACCAATTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGTGACATGCTCCGGCTTTACAGt                                                                                                                                   <  2:912815/150‑1 (MQ=255)
               gTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGt                                                                                                                                   <  1:36867/150‑1 (MQ=255)
                  tACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAtt                                                                                                                                >  1:1798399/1‑150 (MQ=255)
                      cGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                                                                            <  1:1850416/150‑1 (MQ=255)
                        gagaGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTATTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                                                                          <  2:45988/150‑1 (MQ=255)
                          gagTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGTGaa                                                                                                                        <  1:1293703/150‑1 (MQ=255)
                            gTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATa                                                                                                                      <  1:368316/150‑1 (MQ=255)
                             ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATac                                                                                                                     >  1:1780291/1‑150 (MQ=255)
                              tGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCTGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                                                                    >  2:1804781/1‑150 (MQ=255)
                               gATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                                                                   >  2:531740/1‑150 (MQ=255)
                               gATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                                                                   <  2:377968/150‑1 (MQ=255)
                               gATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCACCATACTTTAAAGCCGGAGACATTCTCCGGCTTTACAGTATTAATGAGAAAacac                                                                                                                   >  2:830553/1‑150 (MQ=255)
                                 tGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                 >  2:780148/1‑150 (MQ=255)
                                 tGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                 >  1:208702/1‑150 (MQ=255)
                                   cAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                                                               >  1:890089/1‑150 (MQ=255)
                                     gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                                                             >  2:271782/1‑150 (MQ=255)
                                       ttGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                                                                           <  2:1904067/150‑1 (MQ=255)
                                       ttGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                                                                           >  2:1393533/1‑150 (MQ=255)
                                        tGCCAGTGTCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                                                                          <  2:1600591/150‑1 (MQ=255)
                                         gCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                                                                         >  2:813153/1‑150 (MQ=255)
                                         gCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                                                                         >  2:1816650/1‑150 (MQ=255)
                                           cAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATa                                                                                                       <  2:1521377/150‑1 (MQ=255)
                                             gTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAt                                                                                                     >  1:453767/1‑150 (MQ=255)
                                               ggCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTACCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                                                                   >  1:18495/1‑150 (MQ=255)
                                                 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                                                                 >  1:406486/1‑150 (MQ=255)
                                                  aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                >  1:752041/1‑150 (MQ=255)
                                                   tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                                                               >  2:971501/1‑150 (MQ=255)
                                                    gATTTTATACAGATATTTTATTCTTTTACATGATATAAACCAGCATACCATTTGAGCGCTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGt                                                                                              >  2:1280589/1‑150 (MQ=255)
                                                       tttATACAGATAATTTATTCTTATACATGAAATAAACCAGCATACCATTTGAGCGATCACCCGCTGATCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                                                                           <  1:789406/150‑1 (MQ=255)
                                                        ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                                                                          <  2:263187/150‑1 (MQ=255)
                                                         tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                                                                         <  2:13831/150‑1 (MQ=255)
                                                           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCAGCTGTTCCCCCATACAGTAAAGCCGGAGAAATGCTCAGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAACGAAGGATCAatat                                                                                       >  1:1241504/1‑150 (MQ=255)
                                                            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                                                                      <  1:1674312/150‑1 (MQ=255)
                                                              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                                                                    <  1:1425187/150‑1 (MQ=255)
                                                                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                                                                  <  2:1505364/150‑1 (MQ=255)
                                                                       tttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGATGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                                                                           <  2:824700/149‑1 (MQ=255)
                                                                        ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                                                                          >  2:830836/1‑150 (MQ=255)
                                                                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGATCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGATGGTTCAATATCGCCCCTTTTTGaa                                                                         >  1:168873/1‑149 (MQ=255)
                                                                                cATGAAATAAACCATCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                                                                  <  2:904052/150‑1 (MQ=255)
                                                                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                                                                  <  2:1169665/150‑1 (MQ=255)
                                                                                 aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGCg                                                                 >  1:654365/1‑150 (MQ=255)
                                                                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                                                                >  2:608783/1‑150 (MQ=255)
                                                                                   gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                                                               <  1:602658/150‑1 (MQ=255)
                                                                                    aaaTAAACCAGCATACCATTTGAGCGATCGCCGGCTGTTCCCCCATACTGTAATGCCGGAGACATGCTCCGGCTTTACAGTATTTCTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCGCCTTTTTGATCGCCAGCGttt                                                              >  2:817649/1‑150 (MQ=255)
                                                                                     aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCATTTTTGATCGCCAGCGATTc                                                             >  2:1266233/1‑150 (MQ=255)
                                                                                     aaTAAACCAGCATACAATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                                                             <  1:930159/150‑1 (MQ=255)
                                                                                      aTAAACCAGCATACCATTCGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg                                                            >  2:5603/1‑150 (MQ=255)
                                                                                        aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATAAAACATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc                                                          >  2:1886783/1‑150 (MQ=255)
                                                                                                   ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAAACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTAACCGTGCAGAc                                               >  2:1578349/1‑150 (MQ=255)
                                                                                                      tttGATCGATCGCCCGCTGTTCCCCCAAACTGTAAAGCCGGATACATTCTCCGTCTTTACAGTATTACTGAGAATACACCCTTCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGt                                            <  1:1280805/150‑1 (MQ=255)
                                                                                                      tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGt                                            >  2:1184807/1‑150 (MQ=255)
                                                                                                        tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTg                                          >  2:1088836/1‑150 (MQ=255)
                                                                                                         gAGCGATCGCCCGCTGTTCCCCCATCCTTTACAGCCTGCGACATGCTCCGGCTTTACAGTATTACTGAGAATACAACCATCAAGCATTATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGCTTCGTCACCGTGCCGCCGGTTGc                                         <  1:271782/150‑1 (MQ=255)
                                                                                                          aGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGcc                                        >  2:661040/1‑150 (MQ=255)
                                                                                                           gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCg                                       >  2:82424/1‑150 (MQ=255)
                                                                                                              aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACCCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCCTCACCGTGCAGACGGTTGCCGtcc                                    >  1:1786512/1‑149 (MQ=255)
                                                                                                              aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACCCCCATCAAGCATAATCGAAGGTTCAATATCGCACCAATTTGATCGCCAGCGATTCGACACCCTCCAGACGGTTGCCGgca                                    >  2:1281658/1‑150 (MQ=255)
                                                                                                               tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGACGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcac                                   <  2:462035/150‑1 (MQ=255)
                                                                                                                cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcact                                  >  2:800327/1‑150 (MQ=255)
                                                                                                                cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcact                                  <  1:1300176/150‑1 (MQ=255)
                                                                                                                  cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                >  2:706715/1‑150 (MQ=255)
                                                                                                                  cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                <  1:471483/150‑1 (MQ=255)
                                                                                                                  cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                >  1:652425/1‑150 (MQ=255)
                                                                                                                  cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                <  1:709027/150‑1 (MQ=255)
                                                                                                                  cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAACGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                <  1:1363821/150‑1 (MQ=255)
                                                                                                                   ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactcac                               <  1:129897/150‑1 (MQ=255)
                                                                                                                       tGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaaa                           <  2:1798399/150‑1 (MQ=255)
                                                                                                                          tcgcccATACTTTAAAGCCGGAGACATGCTCAGGCTTTACAGTATTACTGATAATACACCCATCAAGCATAATCGAATGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTa                        <  2:975507/147‑1 (MQ=255)
                                                                                                                           ccaccAAACTGTAAAGCCCGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCTTCAAGCAAAATCGAAGGTTCAATATCGACACTTTTTGATCGCCATCTATTCGTAAACGTGCAGACCCTAGCCGTAACTCACTAAAatat                       >  2:23225/4‑150 (MQ=255)
                                                                                                                              ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATggg                    <  2:92545/150‑1 (MQ=255)
                                                                                                                              ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATggg                    <  1:1707568/150‑1 (MQ=255)
                                                                                                                                 tACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGa                 >  2:1443775/1‑150 (MQ=255)
                                                                                                                                   cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGTTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACa               >  2:1031089/1‑150 (MQ=255)
                                                                                                                                    tGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAg              <  1:495828/150‑1 (MQ=255)
                                                                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGATTACTCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGCCAgg             >  2:1191015/1‑149 (MQ=255)
                                                                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCACTATCGCCCCTTTTTGATCGCCAGCGCTTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGTCAGt             >  1:233147/1‑150 (MQ=255)
                                                                                                                                      tAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTg            >  1:1272639/1‑150 (MQ=255)
                                                                                                                                       aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAATCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa           <  2:1209405/150‑1 (MQ=255)
                                                                                                                                       aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa           <  1:1861902/150‑1 (MQ=255)
                                                                                                                                       aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa           <  2:783025/150‑1 (MQ=255)
                                                                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa           >  2:639965/1‑150 (MQ=255)
                                                                                                                                        aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTgta          >  1:715649/1‑148 (MQ=255)
                                                                                                                                         aGCCGGAGCATGTCTCCGGCTTTACAGTATTTCTGAGAATACACCCATCAAGCATAATCGAAGGATCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTGACGAAATTATGGGGGACAGTGaac         >  2:1270221/1‑150 (MQ=255)
                                                                                                                                          gCCGGAGCCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCGTCACGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGGTGCCGTCGCGCACTAAATTATGGGGGTCGGTGaaga        >  2:1603315/1‑148 (MQ=255)
                                                                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaaca        >  1:20138/1‑150 (MQ=255)
                                                                                                                                           ccGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaacaa       <  1:1422817/150‑1 (MQ=255)
                                                                                                                                              gagaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAg    <  1:1072011/150‑1 (MQ=255)
                                                                                                                                                gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGt  >  1:1576623/1‑150 (MQ=255)
                                                                                                                                                gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGt  >  1:139454/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGT  >  NC_000913/4520355‑4520648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: