Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,240 G→A L21L (TTG→TTA creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,2400GA97.3% 272.4 / ‑0.8 76L21L (TTG→TTAcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (30/43);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

AGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAA  >  NC_000913/4637995‑4638483
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTcc                                                                                                                                                                                                                                                <  1:255589/251‑1 (MQ=255)
agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTcc                                                                                                                                                                                                                                                <  1:97603/251‑1 (MQ=255)
agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCGGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTcc                                                                                                                                                                                                                                                <  2:216009/251‑1 (MQ=255)
   agcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGAt                                                                                                                                                                                                                                             <  2:198981/251‑1 (MQ=255)
       gTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGACGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGATGTTGTTAATTCCGATAAtt                                                                                                                                                                                                                                         >  2:224742/1‑250 (MQ=255)
       gTCTGGGGGTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCTTGTTTGGTGCAGTATTGCTGTTGTTAATACCGATAATg                                                                                                                                                                                                                                         <  2:565225/251‑1 (MQ=255)
              gTTGGCGTTCGTCAGGGAGGTCGACCGGTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGGTGATTc                                                                                                                                                                                                                                  <  2:451159/251‑1 (MQ=255)
                tGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCAGATTCg                                                                                                                                                                                                                                 <  1:154362/250‑1 (MQ=255)
                tGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGGGGCGTGTTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCGTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCGGTTGTTAATTCCGATAATGCTGATTCg                                                                                                                                                                                                                                 <  2:522434/250‑1 (MQ=255)
                 ggCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGAGGCGTGCTGGCCTCGCTTCGACTTAACCGTCCCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGGGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCgg                                                                                                                                                                                                                                >  1:548360/1‑250 (MQ=255)
                   cGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGc                                                                                                                                                                                                                               >  2:522539/1‑249 (MQ=255)
                           aGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGt                                                                                                                                                                                                                     <  2:128779/251‑1 (MQ=255)
                            gTGAGGTCGTCCGTTTGTTTAACGGCGAAGTCACGCGGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTCCTCCTTCACTTCACATAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTGTAAATCCCCCCTGGTCGGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGGGGTTTAAGCCGATAATGCTGATTCGGCATGTGATTGTc                                                                                                                                                                                                                    <  2:538660/251‑1 (MQ=255)
                                        gtttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTGCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAt                                                                                                                                                                                                        <  2:49195/251‑1 (MQ=255)
                                           tgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTa                                                                                                                                                                                                      <  1:255773/250‑1 (MQ=255)
                                               tAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGAGTGCTGGCCTCGCTTAGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGCCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGACTGTTTGGTGCAGTATTGCTGTTGTTAATTCAGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGt                                                                                                                                                                                                  >  2:84992/1‑250 (MQ=255)
                                                aaCGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGt                                                                                                                                                                                                  >  1:550014/1‑249 (MQ=255)
                                                aaCGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTa                                                                                                                                                                                                 <  1:137077/250‑1 (MQ=255)
                                                      gAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCTCTTCGACTTCACCGTCACTTCACATAGCTGCAAATTCTTCCCACATAGTCTTCGTATCCGGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATGCCGATAATGCTTATTCGGGAGGTGAGTGTCGAACGGGCTTATTAc                                                                                                                                                                                                     <  2:244431/240‑1 (MQ=255)
                                                      gAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGGTTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGAGAtgt                                                                                                                                                                                          >  2:348354/1‑251 (MQ=255)
                                                      gAAGTCACGCTGCGCAACGTGCAGGAAGGTGACGTGCTCCCCTCGCTTCGACTTCCCCGTCACTTCACATCGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGCGATTGTCGAACGTGCTGATTAc                                                                                                                                                                                                     >  1:244431/1‑240 (MQ=255)
                                                        aGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTCTGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTgg                                                                                                                                                                                        >  2:463926/1‑251 (MQ=255)
                                                             cGCTGCGCAACGTGCAGGAAGGGTGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGAt                                                                                                                                                                                   <  2:565410/251‑1 (MQ=255)
                                                                    cAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTc                                                                                                                                                                            <  1:201435/251‑1 (MQ=255)
                                                                        gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGcc                                                                                                                                                                         >  1:272094/1‑250 (MQ=255)
                                                                           cAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCaa                                                                                                                                                                       >  1:317192/1‑249 (MQ=255)
                                                                             ggAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGt                                                                                                                                                                    >  1:381942/1‑250 (MQ=255)
                                                                                  gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGc                                                                                                                                                              <  1:522539/251‑1 (MQ=255)
                                                                                             gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaaaa                                                                                                                                                   <  1:504251/251‑1 (MQ=255)
                                                                                              ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaaaa                                                                                                                                                   >  2:81026/1‑250 (MQ=255)
                                                                                                  gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGAtt                                                                                                                                                                                                                       >  1:384079/1‑178 (MQ=255)
                                                                                                  gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGAGTCGGCAGGTGAtt                                                                                                                                                                                                                       <  2:384079/178‑1 (MQ=255)
                                                                                                     tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACAAGCGGGCCGCCAAAACTCGTTg                                                                                                                                           >  2:234608/1‑251 (MQ=255)
                                                                                                        aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGCTGTGGAAAATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGc                                                                                                                                        >  2:121507/1‑251 (MQ=255)
                                                                                                        aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGc                                                                                                                                        >  1:545335/1‑251 (MQ=255)
                                                                                                               cggcACTTCACATAGCTTCAAATTCTTCCCACATATGCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTGGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACGGTAGCGATGTGGGATATGCGATGCGCCAAAGTACCAGCGGGCCTCAAAAACTCGTTTGGGCGCTCAt                                                                                                                                 <  2:351734/248‑1 (MQ=255)
                                                                                                               cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCAt                                                                                                                                 >  2:42495/1‑251 (MQ=255)
                                                                                                                        acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTcc                                                                                                                        <  1:84992/251‑1 (MQ=255)
                                                                                                                         caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTGCCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                       <  1:41803/251‑1 (MQ=255)
                                                                                                                            aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCCAAGGAGAAGACTATGTTGTAATCCCCCATGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCtgt                                                                                                                     <  1:314997/250‑1 (MQ=255)
                                                                                                                                 cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGAGTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGAGGCGATTCGCCAAAGTCCCAGCGGGCCGCAAAAACTCGTTTGGCCGCTCATCGCTATTCCTGTGACCg                                                                                                                <  2:585229/250‑1 (MQ=255)
                                                                                                                                   aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGCAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGc                                                                                                             >  1:361495/1‑251 (MQ=255)
                                                                                                                                    aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCt                                                                                                            <  1:43789/251‑1 (MQ=255)
                                                                                                                                           ccACATAGTCCTCGTAGCCTGCGGCCAGTGCAAAGGAGAAGTCTATGTGGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTGGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACg                                                                                                     <  2:308481/251‑1 (MQ=255)
                                                                                                                                              caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGTTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTg                                                                                                  >  1:326034/1‑251 (MQ=255)
                                                                                                                                               aTAGTCTTCGTAGCCTGCTGCCATTGCAAAGTAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTATCCTGTTTGTTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCTATTCGCCATATTACCAGCGGGCCGCAAAAACTCGTGGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                 <  1:350316/251‑1 (MQ=255)
                                                                                                                                                  gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACTGTGCAgg                                                                                              <  1:463926/251‑1 (MQ=255)
                                                                                                                                                            ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaa                                                                                    <  2:315582/251‑1 (MQ=255)
                                                                                                                                                                   ccATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTg                                                                             <  1:234608/251‑1 (MQ=255)
                                                                                                                                                                      ttGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                            >  1:397902/1‑249 (MQ=255)
                                                                                                                                                                         cAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCTAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCgg                                                                       <  1:42495/251‑1 (MQ=255)
                                                                                                                                                                          aaaGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAACTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTc                                                                                                                                               >  2:351820/1‑178 (MQ=255)
                                                                                                                                                                          aaaGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTc                                                                                                                                               <  1:351820/178‑1 (MQ=255)
                                                                                                                                                                              gagaAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACg                                                                  <  1:314805/251‑1 (MQ=255)
                                                                                                                                                                              gagaAGACTATGGTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGGCCGCTCATCGCTAGTCCTGTGACCGAGCTTTATGCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACg                                                                  <  2:272094/251‑1 (MQ=255)
                                                                                                                                                                              gagaAGACTATGGTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATGCCGATAATGCTGATGCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATGCCTGTGACCGAGCGTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACg                                                                  <  2:99773/251‑1 (MQ=255)
                                                                                                                                                                               agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGCACGTGCTGATTACCGTAGCGAtgt                                                                                                                                                                                          >  1:560536/1‑130 (MQ=255)
                                                                                                                                                                               agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTTCCGTAGCGAtgt                                                                                                                                                                                          <  2:560536/130‑1 (MQ=255)
                                                                                                                                                                                               aTCCTCCCTTTTCTGGAAAATGACTAGCCTGTTTGGCGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCTCACAAACTCGGTGGGCCGCTCATCGATATTCCTGTGCACGCGCTTTATTCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTGTTATCCAttttt                                                 >  2:328519/1‑251 (MQ=255)
                                                                                                                                                                                                 ccccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTg                                                >  1:103357/1‑250 (MQ=255)
                                                                                                                                                                                                 ccccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTg                                                <  2:103357/250‑1 (MQ=255)
                                                                                                                                                                                                   ccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGtt                                             <  2:301634/251‑1 (MQ=255)
                                                                                                                                                                                                     ccttggcTGGAAAATGACTAGCCTGGGTGGTGCAGTATTGCGGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGAGGGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGAGGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTAGTCCTGTGACCGAGCTTTATACGGTGCAGGTATAGTATAAATCCGTGGAGCGGAATCGAAGTTTTATCCATTTTTTGTTAc                                           <  2:333174/245‑1 (MQ=255)
                                                                                                                                                                                                        gTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTTCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTTCCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTGGGTTACCt                                         <  2:217233/250‑1 (MQ=255)
                                                                                                                                                                                                        gTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATCACGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACAGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGATATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                        >  1:478763/1‑251 (MQ=255)
                                                                                                                                                                                                         ttCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAGACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                        >  1:472855/1‑250 (MQ=255)
                                                                                                                                                                                                               aaaaTGACTAGCCTGTTTGGTGCCGTATTGCTGCTGTTAATTCCGATAATGATGATTCGGCAGGTGATTGTCGAACGTGCTGATTACAGTAGCGATGTGGAAGATGCGATTCGACAAAGTATCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACATGAGTCAtt                                 >  1:494730/1‑251 (MQ=255)
                                                                                                                                                                                                                aaaTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGATACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttg                                <  2:548360/251‑1 (MQ=255)
                                                                                                                                                                                                                              tttGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCa                   >  2:392926/1‑250 (MQ=255)
                                                                                                                                                                                                                                tGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGTCCGCt                                                                                                                                    <  1:397642/135‑1 (MQ=255)
                                                                                                                                                                                                                                tGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCt                                                                                                                                    >  2:397642/1‑135 (MQ=255)
                                                                                                                                                                                                                                     cAGTATTGCTGTTGTTACTTCCGCTCATGCTGATTCGGCAAGTGATTGTCGAACGTGCTGATTACCGTAACGATGTGGAAGATGCGATTCGCCTAAGTACCAGAGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaaaa                                                                                  >  2:330841/1‑180 (MQ=255)
                                                                                                                                                                                                                                     cAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAAc           <  1:81026/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGTATTGCTGTGGGTAATTCCGAGAATGCTGATTCGGCAGGTGATTGTCGAACTTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAATTACCAGCGGGCCGCAAAAACTCGTTGAGCGGCTCAGCGATATTCCGGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaaaa                                                                                  <  1:330841/180‑1 (MQ=255)
                                                                                                                                                                                                                                       gTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGt         >  1:548818/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gttgttAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGgaagaa  <  2:317192/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAA  >  NC_000913/4637995‑4638483

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