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

TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTT  >  NC_000913/4546034‑4546522
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
tAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTCACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAATGGAAATGCCCAGCAGGAAAAACCCCCCTTTTTcc                                                                                                                                                                                                                                                   >  2:60421/1‑251 (MQ=255)
 aaaCAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTcc                                                                                                                                                                                                                                                   >  1:79589/1‑250 (MQ=255)
     aaaaTGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAAAGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTCACGATCAAAGAGGCTGACGACAGCACGAAGATTTTTACCGTACCCTATTCGTCCGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGAAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGt                                                                                                                                                                                                                                             >  1:306477/1‑251 (MQ=255)
               tGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGTAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAATGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTACAGGGTACATTTCTcc                                                                                                                                                                                                                                   >  2:212471/1‑251 (MQ=255)
                      tataATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCt                                                                                                                                                                                                                            >  1:328152/1‑251 (MQ=255)
                            aGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGTGAATACCGTAGGGGAAATGCGCAGCAGGAAAAACCCCCCTTTTTCCAGAGGAAATTACTCCACCGCCTTCCggc                                                                                                                                                                                                                      >  2:55259/1‑251 (MQ=255)
                                    gCCACCGGGGCCTTTTCCCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATAACGTAGTGGCAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGctggctggg                                                                                                                                                                                                              >  1:354686/1‑250 (MQ=255)
                                    gCCACCGGGGCCTTTTACCATCAACGATAACTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGa                                                                                                                                                                                                             <  1:494391/251‑1 (MQ=255)
                                       gccGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGAATACCGAAGTGGAAATGCGCAGCAGGGAAAACCCCGCATTTTCCAGAGTACATTACTCCACGTCCTTCCggcgtgccgtgaa                                                                                                                                                                                                            >  2:293082/2‑240 (MQ=255)
                                         cGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGCGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAat                                                                                                                                                                                                         >  1:454719/1‑250 (MQ=255)
                                         cGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTCGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTAACGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGCAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGACCTCCGGCTGCATGGACact                                                                                                                                                                                                         >  1:23552/1‑248 (MQ=255)
                                              ccTTTTTCCATCAACGATATCTATGCCGCAGGTAAGATTGGTGACTGGCAGGTAACGAGCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTGTCCAGAGTACATTACGCCACGGCCTTCCGGCTGGCTGGACAATATAtgg                                                                                                                                                                                                   <  2:475145/251‑1 (MQ=255)
                                                gttttccagcAACGATATCTAGGCCGCAGGTAATAGTGTGGATTTTCAGGTATCGATCAAAGAGGCTGACGGCAGCACGCATATTTTTACCGTTCCCTAATCGTCAGTCCCGCTTTTGCAACGTTAAGGGCATACTCGTTATTCCATGACGGCAGGAGAATACCGTATTGGATATGCGCAGCCGGAAAATCCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCGGGCTGGCCAAGATAtgctc                                                                                                                                                                                                 <  2:280508/242‑4 (MQ=255)
                                                  ttACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATAtggtgg                                                                                                                                                                                                <  2:454719/250‑1 (MQ=255)
                                                    aCCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAAc                                                                                                                                                                                             >  1:291445/1‑251 (MQ=255)
                                                      cATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAACGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGATTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGGCAATATATGGTGGCACGc                                                                                                                                                                                           >  1:313073/1‑251 (MQ=255)
                                                        tCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGAAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGAAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCAACGGCCTTCCGGCTGGCTGGGACATATATGGTGGAACGCa                                                                                                                                                                                          >  1:337118/1‑250 (MQ=255)
                                                            cGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTgg                                                                                                                                                                                     <  1:276331/251‑1 (MQ=255)
                                                                    aTGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGAt                                                                                                                                                                                <  1:55259/248‑1 (MQ=255)
                                                                      gCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATCTATTGTGGAACGCAACTAGCGGATCGtt                                                                                                                                                                            >  1:456234/1‑250 (MQ=255)
                                                                                  aGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCCGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGGGGAACGCACCTGGCGGATCGTTAACGGGCTTTTa                                                                                                                                                                >  1:472900/1‑250 (MQ=255)
                                                                                     ggtgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAAtg                                                                                                                                                             <  2:40447/250‑2 (MQ=255)
                                                                                       tgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCCACTGGCGGCTCGTTATCGTGCTTTTAATTTCg                                                                                                                                                         >  1:78747/1‑251 (MQ=255)
                                                                                       tgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATAACGTAGTGGGAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAAAGTACATTACTCCACGGCCTTCCGGCTGGCTGGGCAATATTTGGGGGAACGCAACTGGCGGGTTGGTTTCGGGCTTTTTATTTCg                                                                                                                                                         >  2:157217/1‑251 (MQ=255)
                                                                                         aCTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCAAGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGt                                                                                                                                                       >  1:434317/1‑251 (MQ=255)
                                                                                         aCTTGCAGGTAACGATCAAAGAGGCGGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATAGGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGt                                                                                                                                                       <  2:152291/251‑1 (MQ=255)
                                                                                            tGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTAACCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTa                                                                                                                                                      <  2:154859/249‑1 (MQ=255)
                                                                                                   aaCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGaaa                                                                                                                                              >  1:123595/1‑250 (MQ=255)
                                                                                                             gagGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGAATACCCTATTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGGGTACATTACTCAACGGCCTTCCGGCTGGCTGGACAATATATGGTGGGACGCCACTGGCGGGTCGTTATTGTGCTTTTAATTTTCGTTTCGGGGAAAACATtggggg                                                                                                                                   >  2:384450/1‑250 (MQ=255)
                                                                                                              gggCTGACGGCTGCACGCAGATTTTTTCCGTGACCTATTCGTCAGTCCCTCTGTTGCAACGTGAAGGGCATAATCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCGGCAGGAAAAACCCCGCTTGTTCCAGAGTACATTACTACTCGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGTTATCGCGAAAACCATGGGGGCa                                                                                                                                  <  2:404617/250‑1 (MQ=255)
                                                                                                                  tGACGGCAGCACGCAGCTTTTTACCGTACCCTATTCGTCCGTCCCGCTTTTGCAACGTGAGGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGCATACATTACTCCACGGCCTTCCGGCTGGCTTTACCATATATGGTGGAACGCCACTGGCGTGTCGTTATCGTGCTTTTAATTTCATTATCGGGAAAAACATggggcacctgg                                                                                                                              >  1:447645/1‑244 (MQ=255)
                                                                                                                    ccGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGgc                                                                                                                            <  2:291445/250‑1 (MQ=255)
                                                                                                                            ccGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGAGCGTTATCGTGCTTTTAATTGCGGTATCGGGAAAAACATGGGGGCACTGGGCGCtctgtc                                                                                                                    <  2:457380/250‑1 (MQ=255)
                                                                                                                                             aCCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAg                                                                                                    <  1:383332/250‑1 (MQ=255)
                                                                                                                                                 tATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCCGGAGAATACCGTAGTGGAAAAGCGCAGCCGGAAAAACCCCCCTTTTTCCAGAGTACATTAATCCACCGCCTTACTGCTGGCTGGAACATATATGGTGGCACGCAACTTGCGGCTCGTTATCGTTCTTTTAATTTTCGGATCGGGAAAAACTTGGGGGCAATGGGCGctgagtcgagcgaacgtcggtgggggg                                                                                                >  2:279424/1‑225 (MQ=255)
                                                                                                                                                 tATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTAGCGGGAAAAACATGGGGGCACTGGGCGct                                                                                                                         <  1:279424/225‑1 (MQ=255)
                                                                                                                                                    tCGTCAGGCCCGCTTTTGCAACGAGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAACACCCTCTTTTTCCAGAGTAAATTACTCCACGGCCTTCCGGCTGGCTGGCCCATTTATTGGGGAGCGCAACAGGAGGGGTGGTATCGTGCTTTTACTTTCTGGGTCGGGGACAACATGGAGGCACATGGCGCTCTGTCTGTGGGTCTGACGCAAGCGAAttt                                                                                            >  2:373533/1‑250 (MQ=255)
                                                                                                                                                      gTCAGTCCCGCGTGTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGAGCGTTATCGTGCTTTTAATTTCGGGATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          <  2:456234/251‑1 (MQ=255)
                                                                                                                                                         aGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGGACGTTATCGTGCTTTTAATTTCGGGATCGGGTAAAACCTGGGGGCACTGGGCGCTCTGTCTGTGGGTATGGAGCAGGCTAATTCCacac                                                                                       >  2:445028/1‑251 (MQ=255)
                                                                                                                                                          gTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACt                                                                                      <  1:445028/251‑1 (MQ=255)
                                                                                                                                                           tCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACtt                                                                                     >  1:309651/1‑251 (MQ=255)
                                                                                                                                                                                               ccATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCg                                                  <  1:157217/250‑1 (MQ=255)
                                                                                                                                                                                                 aTTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCCCGGCCTTCCGGCTGGCTGGGCAATATATGGTGGGACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTAACGGGAAAAACATGGGGGCACCGGGCGCTCTTTCTGTGGGTATGACGCAGGACAATTTCACACTTTCCCATGACAGTCCGCATGAAGGGCAATCCGGGCGttt                                               >  2:205324/1‑251 (MQ=255)
                                                                                                                                                                                                     cGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCAAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGCTGCCAGTCAGCATGACGGACAATAGGTGCGTTTtctc                                           >  2:321837/1‑251 (MQ=255)
                                                                                                                                                                                                         aGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCtata                                       <  1:412833/251‑1 (MQ=255)
                                                                                                                                                                                                          ggAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCAGCACGTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCtata                                       >  1:146021/1‑250 (MQ=255)
                                                                                                                                                                                                            agaATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          >  2:34832/1‑197 (MQ=255)
                                                                                                                                                                                                            agaATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCca                                                                                          <  1:34832/197‑1 (MQ=255)
                                                                                                                                                                                                              aaTACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCAACGGCCTTCCGGCTGGCTGGCCAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAACACAAGGGGGCACTGGGCGCTCTGTCTGTGGTTATTACGCAGGCTAATTCCACACTTCCCGGTGACAGTAAGCATGACGGAAAAACGGAGCGttttatttatacaaa                                   >  2:274210/1‑239 (MQ=255)
                                                                                                                                                                                                               aTACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAAt                                 <  1:306664/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATc                                <  2:79589/250‑1 (MQ=255)
                                                                                                                                                                                                                        gtgaaaagcgacgcagGAAAAACCCCGCTTTTTACAGGAGACATTACGCCACGGCCTGCCGGCTGGCTGGTAAATTGATGGTGTAACGCAACTGGGGGAGCGGTATCGTGCTTTGAATTTCGGTAACGGGTAAAACATGGGGGCACTGGGAGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCTTTTTTTCTATAACAAATCGCTCAATg                        <  2:447645/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaaTGCGCAGCAGGAAAAACCCCTCTTATTCCAGAGTACATTACCCCGCGGCCTTCCGGCTGGTTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCGAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAAt                     <  2:339935/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGtt  >  2:213827/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
TAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTT  >  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: