Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,546,272 A→C T393P (ACC→CCC)  fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,546,2720AC94.3% 144.5 / ‑5.3 53T393P (ACC→CCC) fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (27/23);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTT  >  NC_000913/4546034‑4546522
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
tAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTCACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAATGGAAATGCCCAGCAGGAAAAACCCCCCTTTTTcc                                                                                                                                                                                                                                                   >  2:60421/1‑251 (MQ=255)
 aaaCAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTcc                                                                                                                                                                                                                                                   >  1:79589/1‑250 (MQ=255)
     aaaaTGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAAAGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTCACGATCAAAGAGGCTGACGACAGCACGAAGATTTTTACCGTACCCTATTCGTCCGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGAAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGt                                                                                                                                                                                                                                             >  1:306477/1‑251 (MQ=255)
               tGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGTAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAATGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTACAG‑GGTACATTTCTcc                                                                                                                                                                                                                                   >  2:212471/1‑251 (MQ=255)
                      tataATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCt                                                                                                                                                                                                                            >  1:328152/1‑251 (MQ=255)
                            aGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGTGAATACCGTAGGGGAAATGCGCAGCAGGAAAAACCCCCCTTTTTCCAG‑AGGAAATTACTCCACCGCCTTCCggc                                                                                                                                                                                                                      >  2:55259/1‑251 (MQ=255)
                                    gCCACCGGGGCCTTTTCCCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATAACGTAGTGGCAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGctggctggg                                                                                                                                                                                                              >  1:354686/1‑250 (MQ=255)
                                    gCCACCGGGGCCTTTTACCATCAACGATAACTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑a                                                                                                                                                                                                             <  1:494391/251‑1 (MQ=255)
                                       gccGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGAATACCGAAGTGGAAATGCGCAGCAGGGAAAACCCCGCATTTTCCAG‑AGTACATTACTCCACGTCCTTCCggcgtgccgtgaa                                                                                                                                                                                                            >  2:293082/2‑240 (MQ=255)
                                         cGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGCGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAat                                                                                                                                                                                                         >  1:454719/1‑250 (MQ=255)
                                         cGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTCGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTAACGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGCAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGACCTCCGGCTGCATGG‑ACact                                                                                                                                                                                                         >  1:23552/1‑248 (MQ=255)
                                              ccTTTTTCCATCAACGATATCTATGCCGCAGGTAAGATTGGTGACTGGCAGGTAACGAGCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTGTCCAG‑AGTACATTACGCCACGGCCTTCCGGCTGGCTGG‑ACAATATAtgg                                                                                                                                                                                                   <  2:475145/251‑1 (MQ=255)
                                                gttttccagcAACGATATCTAGGCCGCAGGTAATAGTGTGGATTTTCAGGTATCGATCAAAGAGGCTGACGGCAGCACGCATATTTTTACCGTTCCCTAATCGTCAGTCCCGCTTTTGCAACGTTAAGGGCATACTCGTTATTCCATGACGGCAGGAGAATACCGTATTGGATATGCGCAGCCGGAAAATCCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCGGGCTGG‑CCAAGATAtgctc                                                                                                                                                                                                 <  2:280508/242‑4 (MQ=255)
                                                  ttACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATAtggtgg                                                                                                                                                                                                <  2:454719/250‑1 (MQ=255)
                                                    aCCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAAc                                                                                                                                                                                             >  1:291445/1‑251 (MQ=255)
                                                      cATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAACGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGATTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑GCAATATATGGTGGCACGc                                                                                                                                                                                           >  1:313073/1‑251 (MQ=255)
                                                        tCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGAAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGAAGCAGGAAAAAACCCGCTTTTTCCAG‑AGTACATTACTCAACGGCCTTCCGGCTGGCTGGGACA‑TATATGGTGGAACGCa                                                                                                                                                                                          >  1:337118/1‑250 (MQ=255)
                                                            cGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTgg                                                                                                                                                                                     <  1:276331/251‑1 (MQ=255)
                                                                    aTGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGAt                                                                                                                                                                                <  1:55259/248‑1 (MQ=255)
                                                                      gCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATCTATTGTGGAACGCAACTAGCGGATCGtt                                                                                                                                                                            >  1:456234/1‑250 (MQ=255)
                                                                                  aGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCCGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGGGGAACGCACCTGGCGGATCGTTAACGGGCTTTTa                                                                                                                                                                >  1:472900/1‑250 (MQ=255)
                                                                                     ggtgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAAtg                                                                                                                                                             <  2:40447/250‑2 (MQ=255)
                                                                                       tgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCCACTGGCGGCTCGTTATCGTGCTTTTAATTT‑Cg                                                                                                                                                         >  1:78747/1‑251 (MQ=255)
                                                                                       tgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATAACGTAGTGGGAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAA‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑GCAATATTTGGGGGAACGCAACTGGCGGGTTGGTTTCGGGCTTTTTATTT‑Cg                                                                                                                                                         >  2:157217/1‑251 (MQ=255)
                                                                                         aCTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCAAGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGt                                                                                                                                                       >  1:434317/1‑251 (MQ=255)
                                                                                         aCTTGCAGGTAACGATCAAAGAGGCGGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATAGGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGt                                                                                                                                                       <  2:152291/251‑1 (MQ=255)
                                                                                            tGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTAACCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTa                                                                                                                                                      <  2:154859/249‑1 (MQ=255)
                                                                                                   aaCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGaaa                                                                                                                                              >  1:123595/1‑250 (MQ=255)
                                                                                                             gagGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGAATACCCTATTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑GGTACATTACTCAACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGGACGCCACTGGCGGGTCGTTATTGTGCTTTTAATTT‑TCGTTTCGGGGAAAACATtggggg                                                                                                                                   >  2:384450/1‑250 (MQ=255)
                                                                                                              gggCTGACGGCTGCACGCAGATTTTTTCCGTGACCTATTCGTCAGTCCCTCTGTTGCAACGTGAAGGGCATAATCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCGGCAGGAAAAACCCCGCTTGTTCCAG‑AGTACATTACTACTCGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGTTATCGCGAAAACCATGGGGGCa                                                                                                                                  <  2:404617/250‑1 (MQ=255)
                                                                                                                  tGACGGCAGCACGCAGCTTTTTACCGTACCCTATTCGTCCGTCCCGCTTTTGCAACGTGAGGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑CATACATTACTCCACGGCCTTCCGGCTGGCTTT‑ACCATATATGGTGGAACGCCACTGGCGTGTCGTTATCGTGCTTTTAATTT‑CATTATCGGGAAAAACATggggcacctgg                                                                                                                              >  1:447645/1‑244 (MQ=255)
                                                                                                                    ccGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGgc                                                                                                                            <  2:291445/250‑1 (MQ=255)
                                                                                                                            ccGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGAGCGTTATCGTGCTTTTAATTG‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCtctgtc                                                                                                                    <  2:457380/250‑1 (MQ=255)
                                                                                                                                             aCCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAg                                                                                                    <  1:383332/250‑1 (MQ=255)
                                                                                                                                                 tATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCCGGAGAATACCGTAGTGGAAAAGCGCAGCCGGAAAAACCCCCCTTTTTCCAG‑AGTACATTAATCCACCGCCTTACTGCTGGCTGG‑AACATATATGGTGGCACGCAACTTGCGGCTCGTTATCGTTCTTTTAATTTTCGG‑ATCGGGAAAAACTTGGGGGCAATGGGCGctgagtcgagcgaacgtcggtgggggg                                                                                                >  2:279424/1‑225 (MQ=255)
                                                                                                                                                 tATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTAGCGGGAAAAACATGGGGGCACTGGGCGct                                                                                                                         <  1:279424/225‑1 (MQ=255)
                                                                                                                                                    tCGTCAGGCCCGCTTTTGCAACGAGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAACACCCTCTTTTTCCAG‑AGTAAATTACTCCACGGCCTTCCGGCTGGCTGG‑CCCATTTATTGGGGAGCGCAACAGGAGGGGTGGTATCGTGCTTTTACTTT‑CTGGGTCGGGGACAACATGGAGGCACATGGCGCTCTGTCTGTGGGTCTGACGCAAGCGAAttt                                                                                            >  2:373533/1‑250 (MQ=255)
                                                                                                                                                      gTCAGTCCCGCGTGTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGAGCGTTATCGTGCTTTTAATTT‑CGGGATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          <  2:456234/251‑1 (MQ=255)
                                                                                                                                                         aGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGGACGTTATCGTGCTTTTAATTT‑CGGGATCGGGTAAAACCTGGGGGCACTGGGCGCTCTGTCTGTGGGTATGGAGCAGGCTAATTCCacac                                                                                       >  2:445028/1‑251 (MQ=255)
                                                                                                                                                          gTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACt                                                                                      <  1:445028/251‑1 (MQ=255)
                                                                                                                                                           tCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACtt                                                                                     >  1:309651/1‑251 (MQ=255)
                                                                                                                                                                                               ccATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCg                                                  <  1:157217/250‑1 (MQ=255)
                                                                                                                                                                                                 aTTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCCCGGCCTTCCGGCTGGCTGG‑GCAATATATGGTGGGACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTAACGGGAAAAACATGGGGGCACCGGGCGCTCTTTCTGTGGGTATGACGCAGGACAATTTCACACTTTCCCATGACAGTCCGCATGAAGGGCAATCCGGGCGttt                                               >  2:205324/1‑251 (MQ=255)
                                                                                                                                                                                                     cGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCAAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGCTGCCAGTCAGCATGACGGACAATAGGTGCGTTTtctc                                           >  2:321837/1‑251 (MQ=255)
                                                                                                                                                                                                         aGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCtata                                       <  1:412833/251‑1 (MQ=255)
                                                                                                                                                                                                          ggAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCAGCACGTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCtata                                       >  1:146021/1‑250 (MQ=255)
                                                                                                                                                                                                            agaATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          >  2:34832/1‑197 (MQ=255)
                                                                                                                                                                                                            agaATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          <  1:34832/197‑1 (MQ=255)
                                                                                                                                                                                                              aaTACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCAACGGCCTTCCGGCTGGCTGG‑CCAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAACACAAGGGGGCACTGGGCGCTCTGTCTGTGGTTATTACGCAGGCTAATTCCACACTTCCCGGTGACAGTAAGCATGACGGAAAAACGGAGCGttttatttatacaaa                                   >  2:274210/1‑239 (MQ=255)
                                                                                                                                                                                                               aTACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAAt                                 <  1:306664/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATc                                <  2:79589/250‑1 (MQ=255)
                                                                                                                                                                                                                        gtgaaaagcgacgcagGAAAAACCCCGCTTTTTACAGGAG‑ACATTACGCCACGGCCTGCCGGCTGGCTGG‑TAAATTGATGGTGTAACGCAACTGGGGGAGCGGTATCGTGCTTTGAATTT‑CGGTAACGGGTAAAACATGGGGGCACTGGGAGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCTTTTTTTCTATAACAAATCGCTCAATg                        <  2:447645/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaaTGCGCAGCAGGAAAAACCCCTCTTATTCCAG‑AGTACATTACCCCGCGGCCTTCCGGCTGGTTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCGAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAAt                     <  2:339935/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGtt  >  2:213827/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAG‑AGTACATTACTCCACGGCCTTCCGGCTGGCTGG‑ACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTT‑CGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTT  >  NC_000913/4546034‑4546522

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑