Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,251 C→T 15.3% intergenic (+26/+35) rpmE → / ← yiiX 50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,2510CT15.3% 153.7 / 19.4 67intergenic (+26/+35)rpmE/yiiX50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (5/5);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-04

GAGCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAA  >  NC_000913/4126982‑4127511
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
ggggcaccgggCCTTAAACCAGGGTTTTCCCAGGAAAAAAGTGATTCACCCGAAAAACGAAGAAATTACTGCTAGCTGCTATTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACgg                                                                                                                                                                                                                                                                  <  4:451672/265‑1 (MQ=255)
  gCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACg                                                                                                                                                                                                                                                                   >  2:336142/1‑271 (MQ=255)
   cGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCt                                                                                                                                                                                                                                                              <  2:165237/275‑1 (MQ=255)
                                 taaaaaaGATAGTCACCCAGAATACGAAGAAATTACTGCGAGCTGCTCTGTCGGTAACGGAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAc                                                                                                                                                                                                                                <  1:15663/274‑1 (MQ=255)
                                                   gAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAACCAGCGGCCAt                                                                                                                                                                                                              >  2:221799/1‑275 (MQ=255)
                                                          gaagaaAATACTGCTAACTTATCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGGCCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                                                                                                                                       <  2:479590/275‑1 (MQ=255)
                                                           aagaaATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGc                                                                                                                                                                                                      <  2:41243/275‑1 (MQ=255)
                                                                                        aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGa                                                                                                                                                                                                                           >  2:41785/1‑225 (MQ=255)
                                                                                         aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                                                                                                                                         <  2:212438/274‑1 (MQ=255)
                                                                                                                       aTGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAg                                                                                                                                          <  2:3527/275‑1 (MQ=255)
                                                                                                                         gACCTGCACCTCGACGTGGCCAGCGAGTGCCCCCCTTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTTACCGCTCCAACAAGCGTTTCAACATCACGTGCGGCACGTGAGACCGGTTTTCCGAAACACACTCAGCCCTGCGTTCTTGTTTTGATCGAATCTTTGTCAGCTGACCTCGACACCAGCGGCAATTCTTTGACGATCGCGGCTCGTGGTGGCGCATCGAAACTGGCTACGCGCGAGACCACCGTTACCTCagtcg                                                                                                                                        >  4:468660/1‑271 (MQ=255)
                                                                                                                         gACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGcgccgc                                                                                                                                                                                                                                                                      <  6:68378/149‑1 (MQ=255)
                                                                                                                              gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGGATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAata                                                                                                                                   >  2:491115/1‑275 (MQ=255)
                                                                                                                               aaCCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCCAGACCACCGTTTCCTCAAGCGGAatat                                                                                                                                  >  2:41389/1‑275 (MQ=255)
                                                                                                                                      aCGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc                                                                                                                           >  1:38135/1‑275 (MQ=255)
                                                                                                                                           tgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGGGGTGAGTTGTGGCGCATCCAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATTTTCTTGCCGTAAc                                                                                                                      >  2:309660/1‑275 (MQ=255)
                                                                                                                                            gcagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGACAACAGCGGCCATTCTTTGGCGACCGTGGTGATGTGTTGAGCATCGAAAACCGCTTGTGGCGCGACCACCGTTTCCTCACGCGGGATATTCTTGCCGtgaag                                                                                                                     >  4:325001/1‑271 (MQ=255)
                                                                                                                                               gcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCa                                                                                                                                                                                                                                                                     <  1:84973/128‑1 (MQ=255)
                                                                                                                                                 aaGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                        >  2:547072/1‑155 (MQ=255)
                                                                                                                                                  aGTGCCCCCCGTGTTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGAATTTCCGAAAAGAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCttt                                                                                                               <  2:415776/275‑1 (MQ=255)
                                                                                                                                                   gTGCCACCCGTTCTTCACGGNCNAAAGGCGGGTTGTTCCTACCGGTGGCCGTGTTGACCGCTTCAACANGCGTTTCAACATCCCGGGTAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGAACGGCGCTTTTTTTGNGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc                                                                                                                      <  3:38135/267‑1 (MQ=255)
                                                                                                                                                     gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTa                                                                                                                                                                                                                                    <  4:547067/155‑1 (MQ=255)
                                                                                                                                                     gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAg                                                                                                            >  2:136161/1‑275 (MQ=255)
                                                                                                                                                     gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTCGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGTCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCttt                                                                                                               >  2:292164/1‑272 (MQ=255)
                                                                                                                                                      ccACCCTTTCTGAACTTGCAAAAAGGGTGATGTTGGTACCGGTGGCAGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTGTTCGTGCCTGAATTTTGGGCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGTCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGc                                                                                                           <  4:292160/275‑1 (MQ=255)
                                                                                                                                                        aCCCGTTCTTCACTGGCAAACAGCGTGATGGTGCTATCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCt                                                                                                          <  1:50043/274‑1 (MQ=255)
                                                                                                                                                                ttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAcc                                                                                                                                                                                                                               >  5:469885/1‑149 (MQ=255)
                                                                                                                                                                    ccggCGAACTGCGTGAGGTTGCTACCGTTGGCCTTGTTGACCGCGTCAACAAGCGTTTCAGCACCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCTCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACtt                                                                                              <  4:247232/272‑1 (MQ=255)
                                                                                                                                                                       gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGATTATTCTTGCCGTAACGGTCTTTGAGCTTTGCCTGAACTTGTgg                                                                                          >  2:137001/1‑275 (MQ=255)
                                                                                                                                                                       gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGACGCAAGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGTCCATTCTTTGGCGACAGTGGGGAGTTGTGGCGCATCCAAAACGGCTTGTTGCGAGAACACCGTTTCCGCGAGCGGAATGTTATTTCCGAACGGTTTTTTGTCGTTTGCCTGaaatggttg                                                                                          >  4:430540/1‑268 (MQ=255)
                                                                                                                                                                         aaaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttt                                                                                                                                                                                                                                                           >  6:160691/1‑112 (MQ=255)
                                                                                                                                                                         aaaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttt                                                                                                                                                                                                                                                           <  5:160691/112‑1 (MQ=255)
                                                                                                                                                                               cGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAATCAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCtta                                                                                                  >  2:247236/1‑257 (MQ=255)
                                                                                                                                                                               cGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTg                                                                                                                                                                                                        >  4:528558/1‑157 (MQ=255)
                                                                                                                                                                                 tGATTTTGCTACCGGTGTCCGTGTTGACCGCTTCAACAAGCGTTTCGACATCCCGGGCGGCAAATAAGACCTGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTGGTGGCTGAGTTTTGGTCGTATTACCACGAAAACAGCGGCCAt                                                                                                                                                                                                              >  8:115520/1‑149 (MQ=255)
                                                                                                                                                                                 tGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCGAt                                                                                                                                                                                                              >  5:271059/1‑149 (MQ=255)
                                                                                                                                                                                 tGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGTCGAATCGAAAACGTCTTGTGGCGAGAGCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTACTTTGATCTTTGCCTGAACTTGTTGGCTGGTCAg                                                                                 >  1:136562/1‑274 (MQ=255)
                                                                                                                                                                                   aTGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGa                                                                                                                                                                                                    <  2:528563/157‑1 (MQ=255)
                                                                                                                                                                                       tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATCCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGaa                                                                           <  2:337906/274‑1 (MQ=255)
                                                                                                                                                                                         cTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCg                                                                                                                          >  2:327558/1‑225 (MQ=255)
                                                                                                                                                                                          tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATGCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                                                                                                                                     <  7:115520/149‑1 (MQ=255)
                                                                                                                                                                                             cGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgtg                                                                                                                                                                                                                                                     <  6:12099/98‑1 (MQ=255)
                                                                                                                                                                                             cGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgtg                                                                                                                                                                                                                                                     >  5:12099/1‑98 (MQ=255)
                                                                                                                                                                                              ggTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATATTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGTCCACCGTTTCCTCTAGCGGACTTTTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCAGCTCCGGTCGAACACTTTc                                                                    >  2:124617/1‑274 (MQ=255)
                                                                                                                                                                                               gTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCCAGCGGAATATTCTTGCCGTAAAGTTCTTTGAGCTTTGCCTGAACTTTTGGGCTGCGCAGGTCGAACTCTTTCa                                                                   >  2:481716/1‑274 (MQ=255)
                                                                                                                                                                                                                 taatcgatcGTGTCAACAACCGGGGCGGGCAATAGAACCGGATTTCCTAAAAAAAGGCGGTACGGCGCTTGTTTTGTGCCTGAATTTTGGTCTGAGTACCACGAAAACAGCTGCCATTCTTTTGCGATCGTGTTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCGCGCCCa                                                 <  4:124616/266‑1 (MQ=255)
                                                                                                                                                                                                                           tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                                                                                                                                    <  6:469885/149‑1 (MQ=255)
                                                                                                                                                                                                                            ttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCtt                                                                                                                                                                   >  6:248458/1‑149 (MQ=255)
                                                                                                                                                                                                                               aaCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtgtg                                                                                                                                                                >  5:477586/1‑149 (MQ=255)
                                                                                                                                                                                                                                aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGAGGCGCATCGAAAACGGCTTGTgg                                                                                                                                                               >  7:265648/1‑149 (MQ=255)
                                                                                                                                                                                                                                 cATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGcgcg                                <  2:430545/275‑1 (MQ=255)
                                                                                                                                                                                                                                  ataccGGGTATCAAATAAGACCGGATTTCCTAAAAAAGCGCCGGACGGCGCTGTTTTTGTGCCTGAATTTTGGTAGTATTACCACGAAAACAGCAGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGATCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGt                               <  1:5227/272‑1 (MQ=255)
                                                                                                                                                                                                                                  agcccGGGCAGCAAATAAGACCGGTTTTCCGAAAAAAGCGCCGCACGGCTCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATGCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGt                               <  1:107557/273‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTTGCACATCTAAAACGGCTTGTGGCGAGACAACCGTTTCCTCCAGCGGAATATTCTTGCCGTACCCTTCTTTGAGCTTTGCCTGAACTTGTGCGCTGCCCTGGTCTAAGTCTTTCCACTTCTGCTTGTCTCCCTCGCCCATTTCCAGCTCGTTCTGGTACACCTTTCAcacc           >  4:146156/1‑274 (MQ=255)
                                                                                                                                                                                                                                                         gATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAcc                                                                                                                                                                                                  <  6:219398/89‑1 (MQ=255)
                                                                                                                                                                                                                                                         gATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAcc                                                                                                                                                                                                  >  5:219398/1‑89 (MQ=255)
                                                                                                                                                                                                                                                             gcctaaaaagGGGACGCACGGCGCTTTTTTTGTGCTGGAATTTAGGGCGTAGTACACCGAAAACAGAGGCCAGGCTTGGGCGACCGGGGTGAGTTGTGGCGCATCGAAAACGGCTTGGGGCGAGACCACGTTTTCCTCAAGCGGAATATTGTGGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACGCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTGTGGTAAACCTTCCACACCACTTCCg    <  4:309656/271‑1 (MQ=255)
                                                                                                                                                                                                                                                               cGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGaa  <  2:325005/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc                                                                                                                            <  6:271059/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc                                                                                                                            <  5:248458/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                                                                                                                                               <  6:345643/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                                                                                                                                               >  5:345643/1‑128 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                              <  7:299486/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                              >  8:299486/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                     aaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCg                                                                                                                          >  8:72120/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                              <  7:181541/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                              >  8:181541/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                       gcgcCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTa                                                                                                                        >  5:105767/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc                                                                                                                      <  7:72120/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGtt                                                                                                                   >  8:176856/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                    
GAGCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAA  >  NC_000913/4126982‑4127511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: