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

GGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAAT  >  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)
                                                     tgCCTTTTACAATCAACGATCTCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACAATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACAGTACCCTATTCGTCAGTCCCGCTTTTGCAACGGGAAGGGCATACTCGTTATTCCATTACGGCAGGGAAATACCGTTGTGGGAAAGCGCAGCTGGAAAAAACACGCATTTTGCAGAGAACAGTACTCCAACGCCCTCCGGATGGCTGGAC‑Aattatt                                                                                                                                                                                                >  2:462790/2‑247 (MQ=255)
                                                         ttttACAATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCCGAAAGAAAAACCCCCCTTTTTCCAGAGAACATTACTCCACGGCCATCCGGCTGGGTGGAC‑CATATAtggtg                                                                                                                                                                                            >  2:104256/1‑251 (MQ=255)
                                                                  cAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTACCGATCAAAGCGGCTGACGGCAGCACGCAGATTTTTACCGTACACTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAACACCCCGCTTTTTCCAGAGAACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAAc                                                                                                                                                                                   >  1:287831/1‑251 (MQ=255)
                                                                    aCGATATCTATGCCGAAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCATCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATCCTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTg                                                                                                                                                                                 <  2:329219/251‑1 (MQ=255)
                                                                        tatCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCgg                                                                                                                                                                             >  1:73809/1‑251 (MQ=255)
                                                                               gCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTa                                                                                                                                                                      <  1:423425/251‑1 (MQ=255)
                                                                                      gTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCt                                                                                                                                                               >  1:563856/1‑251 (MQ=255)
                                                                                        aaTAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCttt                                                                                                                                                             <  1:140206/251‑1 (MQ=255)
                                                                                             gggtgACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACGCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTTCTTTTAAtt                                                                                                                                                        <  2:389136/250‑1 (MQ=255)
                                                                                                  aCTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTGTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGt                                                                                                                                                   <  1:239952/251‑1 (MQ=255)
                                                                                                       cAGGTAACGATCAAAGCGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCgg                                                                                                                                              >  1:586775/1‑251 (MQ=255)
                                                                                                            aaCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGa‑aa                                                                                                                                         >  1:583535/1‑250 (MQ=255)
                                                                                                             aCGATCAAAGAGGCTGACGGCAGCACGCAGATTTTGACCGGTCCCTATTCGTCAGTCCCGCTTTTGCAACGGTAAGGGCATACTCGGTAATCCATTACGTCAGGAGGATGCCGTAGTTGAAATTCGCAGCAGGAAAAACCCCGCTGTTTCCAGAGTACATTACTCCACGGCCTTCCGGCGGGCTGGAC‑AATATATTGTGTATCGCAACTGGCGGATCGGTATCGTGCTTTTAATTTCGGTATCGGGt‑aaaa                                                                                                                                       <  2:338951/251‑1 (MQ=255)
                                                                                                               taTCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTCTTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGACTTCCGTAGTGGAAATACGCAGCCAGAACAACCCCGCTTTTTCCAGAGTACATTCCTCCCCGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCCACTGGCGGCTCATTTTCTTGCTTTTAGTTTCCGTATCgatcggc                                                                                                                                         >  2:246692/2‑242 (MQ=255)
                                                                                                               gATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATGCGTCAGGCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTAt                                                                                                                                                 <  1:246692/240‑1 (MQ=255)
                                                                                                                aTCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATAGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑t                                                                                                                                   >  1:601115/1‑251 (MQ=255)
                                                                                                                  cAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑Tgg                                                                                                                                 >  1:302155/1‑251 (MQ=255)
                                                                                                                       agGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCGTTTTCCAGAGTACATTACTCCCCGGCCTTCCGGCTGGCTGGAC‑AATATATGGGGGAACGCCACTGGCGGATCGTTATCGTGCTTTTAATTTCGTGATCGGGA‑CAAACAGTggggc                                                                                                                              >  1:294201/1‑249 (MQ=255)
                                                                                                                        ggCTGACGGCAGCACGCAGATTTTTACAGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGAATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGGCCATTCCTCCACGGCCTTCCGGCTGGCTGGGC‑CATATATGGTGGAACGACACTGGCGGATCGTTATTGTGCTTTTAATTTCGGACTCGGGC‑AAAACA‑TGGGGGCAc                                                                                                                           >  1:625996/1‑251 (MQ=255)
                                                                                                                         gCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCATGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACt                                                                                                                          <  1:140207/251‑1 (MQ=255)
                                                                                                                          gtGACGGCAGCACGCGGCTTTTTACCTTACCCTATTCGTCTGTCCCGCTTTTGCAACGTGAAGGGCATACTCTTTATTCCATTGCGGCAGGAGAATACCGTTGTGGAAATGCGCAGCAGGAAAAGCCCCGCTTTTTCCAGATTACATTACTCCACGGCCTTCCGGCTGGCTTGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGC‑AAATCA‑TGGGGTCACTg                                                                                                                         <  2:625996/250‑1 (MQ=255)
                                                                                                                               ggCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGC‑AAAACA‑TGGGGGCACTGGgcgc                                                                                                                    >  1:347500/1‑251 (MQ=255)
                                                                                                                                   gcaCGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGCAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGACCTGCGCAGAAGGCAAAACACCGCTTTTTCCAGAGACCATTACTCAACGGCCTTCCGGCTGGCTGGAC‑AATATCTGGAGGAACGCAACTGGCGGATCGTTCTCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGctctt                                                                                                                >  1:95471/1‑250 (MQ=255)
                                                                                                                                    caCGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCCCGGCCTTCCGGCTGGCTGGCC‑AATATATGGTGGTACGCCAATGGCCGGTCGTTATCGGGTTTTTAATTTCGGTATCGGGA‑GAAACA‑TGGGGGACCTGGGCGgtcttt                                                                                                               >  2:180064/1‑249 (MQ=255)
                                                                                                                                        cAGATTTTTACCTTACCCTATTCGTCAGTCCCTCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATGTCGGTATCGGGA‑AAATCA‑TGGGGGCACTGGGCGCtctgtct                                                                                                             <  2:9102/249‑1 (MQ=255)
                                                                                                                                         ccttcactccgctcttccgatcttCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑GGGGGGCACTGGGCGCTCTGTCtgtg                                                                                                          <  2:552348/228‑1 (MQ=255)
                                                                                                                                          gATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTgg                                                                                                         <  1:111008/251‑1 (MQ=255)
                                                                                                                                           aTTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATc                                                                                                                                                                          <  1:517018/187‑1 (MQ=255)
                                                                                                                                             ttttACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGat                                                                                                       >  1:193397/1‑250 (MQ=255)
                                                                                                                                                      aCCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAgg                                                                                             >  1:94618/1‑251 (MQ=255)
                                                                                                                                                               gTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACACCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑CATATATGGTGGAACGCAACAGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTcc                                                                                    >  1:308298/1‑250 (MQ=255)
                                                                                                                                                                tCAGTCCCGCTTTTGCCACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCtgtg                                                                                                          >  1:552348/1‑228 (MQ=255)
                                                                                                                                                                    tCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACtt                                                                              >  1:138248/1‑251 (MQ=255)
                                                                                                                                                                          ttttGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGctggctgg                                                                                                                                                                                                          <  1:392088/125‑1 (MQ=255)
                                                                                                                                                                           tttGCAACTTGAAGGACATACGCGTTATTCAATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTGTAATTTCGTTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGGGGATATGCCGCAGGC‑TAATTCcatacgtgccgat                                                                        <  2:446507/250‑12 (MQ=255)
                                                                                                                                                                             tGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGa                                                                      <  1:555855/250‑1 (MQ=255)
                                                                                                                                                                              gCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGAcc                                                                    >  1:482626/1‑250 (MQ=255)
                                                                                                                                                                               cAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCAGCTCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCCCTTGGCGCTCTGTCTGTGGATATGGAGCAGGC‑TAATTCCACACTTCCCGATGAc                                                                     >  1:419397/1‑249 (MQ=255)
                                                                                                                                                                                     gAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGGGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TCATTCCACACTTCCCGCTGACAGTCAga                                                              >  1:286925/1‑249 (MQ=255)
                                                                                                                                                                                       aGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCAACGGCCTTCCGGCTGGCTGGAA‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGC‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGACAGCATg                                                           >  1:386349/1‑251 (MQ=255)
                                                                                                                                                                                           cATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACgg                                                       <  1:104256/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTATTCCATTACGGCAGGAGTAGACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGTGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACATTTCCCGATGACAGTCAGCATGACGG‑ACAATCg                                               <  1:554294/251‑1 (MQ=255)
                                                                                                                                                                                                  gTTATTCCATTACGGCAGGAGAATACCGTAGTGGAACTGCGCATCAGGAAATAACCCGCTTTTTACAGAGTACATTACTCCACGGCCCTCCGGCTGGCTGGACCAAT‑TATGGTGGAACGCCACTGGCGGCTCGTTTTCGGGATTTTAATTTCGGGCTAGGGACAAAAC‑‑TGGGTGCACGGGGCGTTCTCTCTCGGGGTATCACGAAGGCATAATTC‑ACACTTTCCGCTGACAGTCAGCATGACGC‑ACAATc                                                >  2:335968/1‑250 (MQ=255)
                                                                                                                                                                                                        ccATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AACACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCCGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATAGGTGCgg                                         >  1:316111/1‑250 (MQ=255)
                                                                                                                                                                                                                               aGTGGAACTACGCAGCAGGAAAAACCCAGAATTTTCCAGAGTACATTAATCCACGGCATTACGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGAATAGGGA‑AGAAAA‑TGGGGGCAATGGACACTGTGTGTGTGGATATGACGCATGC‑TCATTCAACACTTCCCGATGACAGTCAGCATGACGGGACA‑TCGGTGAGCTTTATCTCTAACAAATCGCtac                   >  2:674367/1‑248 (MQ=255)
                                                                                                                                                                                                                                 tGGAAATGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATg                >  1:620148/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aTGCGCAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCa           <  1:4675/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcAGCAGGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAgg         <  2:338110/251‑1 (MQ=255)
                                                                                                                                                                                                                                               agGAAAAACCCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTTGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAat  <  2:563856/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GGTCACTATTAAACAAAATGGGTATGACATTTATAATAGTACGGTGCCACCGGGGCCTTTTACCATCAACGATATCTATGCCGCAGGTAATAGTGGTGACTTGCAGGTAACGATCAAAGAGGCTGACGGCAGCACGCAGATTTTTACCGTACCCTATTCGTCAGTCCCGCTTTTGCAACGTGAAGGGCATACTCGTTATTCCATTACGGCAGGAGAATACCGTAGTGGAAATGCGCAGCAGGAAAAAACCCGCTTTTTCCAGAGTACATTACTCCACGGCCTTCCGGCTGGCTGGAC‑AATATATGGTGGAACGCAACTGGCGGATCGTTATCGTGCTTTTAATTTCGGTATCGGGA‑AAAACA‑TGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGC‑TAATTCCACACTTCCCGATGACAGTCAGCATGACGG‑ACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAAT  >  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: ‑