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,2720AC93.5% 166.0 / NA 62T393P (ACC→CCC) fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base C (31/27);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAAT  >  NC_000913/4546025‑4546514
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
ggTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAcccc                                                                                                                                                                                                                                                      <  1:402556/251‑1 (MQ=255)
        ttAAACAAAATGGGAATGACATTTATAATAGTACGGCGCCACCGGAGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGAGCATACTCGTTATTCCATTACAGCAGGAGAATACCGAAAAGGAAATGAGCAGCCGGACAAAACCCGCTTTTTc                                                                                                                                                                                                                                              >  2:579563/1‑251 (MQ=255)
              aaaaTGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGt                                                                                                                                                                                                                                        >  1:185064/1‑251 (MQ=255)
                   gggTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGGAATGCGCAGCAGGAAAAACCCCGCTTTTTCCa                                                                                                                                                                                                                                            >  2:288555/1‑242 (MQ=255)
                   gggTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCa                                                                                                                                                                                                                                            <  1:288555/242‑1 (MQ=255)
                                ataataGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCtt                                                                                                                                                                                                                      >  1:444674/1‑251 (MQ=255)
                                 tactaGTACGGTGCCACCGGGGCCTTTTACCATCAACGATCTCTATGCCGCAGGTAATAGTGGAGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGAGCATACTCGTTATTCCATTACGGCAGGAGAATACCGGCGGGGACATGAGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTc                                                                                                                                                                                                                     >  2:408371/4‑251 (MQ=255)
                                  aataGTACGGTGCCACAGGGGCCTTTTACCATCAACGATATCTATGCAGCAGGTAATAGTGGTGACTTGCAGGTAACGAGCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGccctc                                                                                                                                                                                                                     <  2:330641/250‑4 (MQ=255)
                                        aCGGTGCCACCGGGGCCTTTTTCCAACAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCGGATTTTTACCGTACCCTATTCGTCAGTCCCTCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTTGAAATGCGCAGCATGAAAAACCCCTCTTTTTCCAGAGTACATTACTCCACGTCCTTCCggctg                                                                                                                                                                                                               <  2:435472/250‑1 (MQ=255)
                                        aCGGTGCCACCGGGGCCTTTTACCATCCACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCggctgg                                                                                                                                                                                                              <  2:368621/251‑1 (MQ=255)
                                        aCGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCggctg                                                                                                                                                                                                               <  1:461878/250‑1 (MQ=255)
                                         cGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGACATGCGCAGCAGGGAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCggctggc                                                                                                                                                                                                             >  1:329219/1‑251 (MQ=255)
                                                aCCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGGGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCg                                                                                                                                                                                                                   >  2:119308/1‑238 (MQ=255)
                                                aCCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCg                                                                                                                                                                                                                   <  1:119308/238‑1 (MQ=255)
                                                     tgCCTTTTACAATCAACGATCTCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACAATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACAGTACCCTATTCGTCAGTCCCGCTTTTGCAACGGGAAGGGCATACTCGTTATTCCATTACGGCAGGGAAATACCGTTGTGGGAAAGCGCAGCTGGAAAAAACACGCATTTTGCAGAGAACAGTACTCCAACGCCCTCCGGATGGCTGGACAattatt                                                                                                                                                                                                >  2:462790/2‑247 (MQ=255)
                                                         ttttACAATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCCGAAAGAAAAACCCCCCTTTTTCCAGAGAACATTACTCCACGGCCATCCGGCTGGGTGGACCATATAtggtg                                                                                                                                                                                            >  2:104256/1‑251 (MQ=255)
                                                                  cAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTACCGATCAAAGCGGCTGACGGCAGCACGCAGATTTTTACCGTACACTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAACACCCCGCTTTTTCCAGAGAACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAAc                                                                                                                                                                                   >  1:287831/1‑251 (MQ=255)
                                                                    aCGATATCTATGCCGAAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCATCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATCCTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTg                                                                                                                                                                                 <  2:329219/251‑1 (MQ=255)
                                                                        tatCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCgg                                                                                                                                                                             >  1:73809/1‑251 (MQ=255)
                                                                               gCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTa                                                                                                                                                                      <  1:423425/251‑1 (MQ=255)
                                                                                      gTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCt                                                                                                                                                               >  1:563856/1‑251 (MQ=255)
                                                                                        aaTAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCttt                                                                                                                                                             <  1:140206/251‑1 (MQ=255)
                                                                                             gggtgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACGCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTTCTTTTAAtt                                                                                                                                                        <  2:389136/250‑1 (MQ=255)
                                                                                                  aCTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTGTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGt                                                                                                                                                   <  1:239952/251‑1 (MQ=255)
                                                                                                       cAGGTAACGATCAAAGCGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCgg                                                                                                                                              >  1:586775/1‑251 (MQ=255)
                                                                                                            aaCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGaaa                                                                                                                                         >  1:583535/1‑250 (MQ=255)
                                                                                                             aCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTGACCGGTCCCTATTCGTCAGTCCCGCTTTTGCAACGGTAAGGGCATACTCGGTAATCCATTACGTCAGGAGGATGCCGTAGTTGAAATTCGCAGCAGGAAAAACCCCGCTGTTTCCAGAGTACATTACTCCACGGCCTTCCGGCGGGCTGGACAATATATTGTGTATCGCAACTGGCGGATCGGTATCGTGCTTTTAATTTCGGTATCGGGtaaaa                                                                                                                                       <  2:338951/251‑1 (MQ=255)
                                                                                                               taTCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTCTTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGACTTCCGTAGTGGAAATACGCAGCCAGAACAACCCCGCTTTTTCCAGAGTACATTCCTCCCCGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCCACTGGCGGCTCATTTTCTTGCTTTTAGTTTCCGTATCgatcggc                                                                                                                                         >  2:246692/2‑242 (MQ=255)
                                                                                                               gATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATGCGTCAGGCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTAt                                                                                                                                                 <  1:246692/240‑1 (MQ=255)
                                                                                                                aTCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATAGTGCTTTTAATTTCGGTATCGGGAAAAACAt                                                                                                                                   >  1:601115/1‑251 (MQ=255)
                                                                                                                  cAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATgg                                                                                                                                 >  1:302155/1‑251 (MQ=255)
                                                                                                                       agGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCGTTTTCCAGAGTACATTACTCCCCGGCCTTCCGGCTGGCTGGACAATATATGGGGGAACGCCACTGGCGGATCGTTATCGTGCTTTTAATTTCGTGATCGGGACAAACAGTggggc                                                                                                                              >  1:294201/1‑249 (MQ=255)
                                                                                                                        ggCTGACGGCAGCACGCAGATTTTTACAGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGAATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGGCCATTCCTCCACGGCCTTCCGGCTGGCTGGGCCATATATGGTGGAACGACACTGGCGGATCGTTATTGTGCTTTTAATTTCGGACTCGGGCAAAACATGGGGGCAc                                                                                                                           >  1:625996/1‑251 (MQ=255)
                                                                                                                         gCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCATGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACt                                                                                                                          <  1:140207/251‑1 (MQ=255)
                                                                                                                          gtGACGGCAGCACGCGGCTTTTTACCTTACCCTATTCGTCTGTCCCGCTTTTGCAACGTGAAGGGCATACTCTTTATTCCATTGCGGCAGGAGAATACCGTTGTGGAAATGCGCAGCAGGAAAAGCCCCGCTTTTTCCAGATTACATTACTCCACGGCCTTCCGGCTGGCTTGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGCAAATCATGGGGTCACTg                                                                                                                         <  2:625996/250‑1 (MQ=255)
                                                                                                                               ggCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGCAAAACATGGGGGCACTGGgcgc                                                                                                                    >  1:347500/1‑251 (MQ=255)
                                                                                                                                   gcaCGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGCAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGACCTGCGCAGAAGGCAAAACACCGCTTTTTCCAGAGACCATTACTCAACGGCCTTCCGGCTGGCTGGACAATATCTGGAGGAACGCAACTGGCGGATCGTTCTCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGctctt                                                                                                                >  1:95471/1‑250 (MQ=255)
                                                                                                                                    caCGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCCCGGCCTTCCGGCTGGCTGGCCAATATATGGTGGTACGCCAATGGCCGGTCGTTATCGGGTTTTTAATTTCGGTATCGGGAGAAACATGGGGGACCTGGGCGgtcttt                                                                                                               >  2:180064/1‑249 (MQ=255)
                                                                                                                                        cAGATTTTTACCTTACCCTATTCGTCAGTCCCTCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATGTCGGTATCGGGAAAATCATGGGGGCACTGGGCGCtctgtct                                                                                                             <  2:9102/249‑1 (MQ=255)
                                                                                                                                         ccttcactccgctcttccgatcttCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACAGGGGGGCACTGGGCGCTCTGTCtgtg                                                                                                          <  2:552348/228‑1 (MQ=255)
                                                                                                                                          gATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTgg                                                                                                         <  1:111008/251‑1 (MQ=255)
                                                                                                                                           aTTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATc                                                                                                                                                                          <  1:517018/187‑1 (MQ=255)
                                                                                                                                             ttttACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGat                                                                                                       >  1:193397/1‑250 (MQ=255)
                                                                                                                                                      aCCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAgg                                                                                             >  1:94618/1‑251 (MQ=255)
                                                                                                                                                               gTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACACCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACCATATATGGTGGAACGCAACAGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTcc                                                                                    >  1:308298/1‑250 (MQ=255)
                                                                                                                                                                tCAGTCCCGCTTTTGCCACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCtgtg                                                                                                          >  1:552348/1‑228 (MQ=255)
                                                                                                                                                                    tCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACtt                                                                              >  1:138248/1‑251 (MQ=255)
                                                                                                                                                                          ttttGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGctggctgg                                                                                                                                                                                                          <  1:392088/125‑1 (MQ=255)
                                                                                                                                                                           tttGCAACTTGAAGGACATACGCGTTATTCAATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTGTAATTTCGTTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGGGGATATGCCGCAGGCTAATTCcatacgtgccgat                                                                        <  2:446507/250‑12 (MQ=255)
                                                                                                                                                                             tGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGa                                                                      <  1:555855/250‑1 (MQ=255)
                                                                                                                                                                              gCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGAcc                                                                    >  1:482626/1‑250 (MQ=255)
                                                                                                                                                                               cAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCAGCTCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCCCTTGGCGCTCTGTCTGTGGATATGGAGCAGGCTAATTCCACACTTCCCGATGAc                                                                     >  1:419397/1‑249 (MQ=255)
                                                                                                                                                                                     gAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGGGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTCATTCCACACTTCCCGCTGACAGTCAga                                                              >  1:286925/1‑249 (MQ=255)
                                                                                                                                                                                       aGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCAACGGCCTTCCGGCTGGCTGGAAAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGCAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGACAGCATg                                                           >  1:386349/1‑251 (MQ=255)
                                                                                                                                                                                           cATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACgg                                                       <  1:104256/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTATTCCATTACGGCAGGAGTAGACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGTGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACATTTCCCGATGACAGTCAGCATGACGGACAATCg                                               <  1:554294/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTATTCCATTACGGCAGGAGAATACCGTAGTGGAACTGCGCATCAGGAAATAACCCGCTTTTTACAGAGTACATTACTCCACGGCCCTCCGGCTGGCTGGACCAATTATGGTGGAACGCCACTGGCGGCTCGTTTTCGGGATTTTAATTTCGGGCTAGGGACAAAACTGGGTGCACGGGGCGTTCTCTCTCGGGGTATCACGAAGGCATAATTCACACTTTCCGCTGACAGTCAGCATGACGCACAATc                                                >  2:335968/1‑250 (MQ=255)
                                                                                                                                                                                                        ccATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAACACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCCGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATAGGTGCgg                                         >  1:316111/1‑250 (MQ=255)
                                                                                                                                                                                                                               aGTGGAACTACGCAGCAGGAAAAACCCAGAATTTTCCAGAGTACATTAATCCACGGCATTACGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGAATAGGGAAGAAAATGGGGGCAATGGACACTGTGTGTGTGGATATGACGCATGCTCATTCAACACTTCCCGATGACAGTCAGCATGACGGGACATCGGTGAGCTTTATCTCTAACAAATCGCtac                   >  2:674367/1‑248 (MQ=255)
                                                                                                                                                                                                                                 tGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATg                >  1:620148/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aTGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCa           <  1:4675/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAgg         <  2:338110/251‑1 (MQ=255)
                                                                                                                                                                                                                                               agGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTTGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAat  <  2:563856/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGACAATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAAT  >  NC_000913/4546025‑4546514

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: