Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,541,261 C→T A102V (GCT→GTT)  fimB → Type 1 fimbriae regulatory protein FimB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,541,2610CT98.3% 200.3 / NA 60A102V (GCT→GTT) fimBType 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/31);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGT  >  NC_000913/4541027‑4541505
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
cagcaAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGt                                                                                                                                                                                                                                       >  2:638139/1‑250 (MQ=255)
      aTACCGGGCCTCATGCAGCACGTAATTATTATCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTAttttt                                                                                                                                                                                                                                 >  1:681837/1‑250 (MQ=255)
        aCCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTAt                                                                                                                                                                                                                               <  1:164680/250‑1 (MQ=255)
                     cagcaCGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTACGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAAt                                                                                                                                                                                                                 >  1:34244/1‑251 (MQ=255)
                                tattGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCg                                                                                                                                                                                                      >  1:434090/1‑251 (MQ=255)
                                tattGTCTGACTTTGCTTGGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGAGTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAGCTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGGATTTTTATCACGTAAGGGGAATCCGCTTtctc                                                                                                                                                                                                       <  2:535517/250‑1 (MQ=255)
                                  ttGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCAATCGCTTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGCATCCGCTTTCTCGGc                                                                                                                                                                                                    >  1:443260/1‑251 (MQ=255)
                                            ttgctttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGtttt                                                                                                                                                                                           <  1:47681/250‑1 (MQ=255)
                                                tttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCat                                                                                                                                                                                      >  2:591855/1‑251 (MQ=255)
                                                        aTTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAtctc                                                                                                                                                                              <  1:638139/251‑1 (MQ=255)
                                                              ggTTGCCGGGCGAGTGAAATTTGTCTATTGAAGATTTCGGATATTGATCTTAAGGCAAAGGGTATATATATCCATCGATTAAAAAAAGGCTTTGCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTc                                                                                                                                                                        <  2:409088/251‑1 (MQ=255)
                                                                tttCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTCAGGCAAAGTGTATATATATCCATCGATTAAACAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATCTTATCTCGACTTCCg                                                                                                                                                                      >  1:657858/1‑251 (MQ=255)
                                                                    cGGGCGAGGGAAATTTGTCGATTGAGGATTTCGGCTATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCAACCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGCATCCGCTTTCTCGGCACCAGTTTTACCATATTATCTCGACTTCCggtgg                                                                                                                                                                  >  2:603022/1‑251 (MQ=255)
                                                                          aGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGc                                                                                                                                                            <  1:399888/251‑1 (MQ=255)
                                                                           gTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCTCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAGCTCGACTTCCGGTGGTAATGcc                                                                                                                                                           <  2:651363/251‑1 (MQ=255)
                                                                                         ttGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                                                              >  2:200833/1‑250 (MQ=255)
                                                                                          tGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCAATCGATTCAAAAAAGGCTTTTAAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAACAAACTGGTTGAGTATCCGTACTTCGTACCAGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                                                              >  1:77506/1‑249 (MQ=255)
                                                                                            aGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGAGTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGTGTATTTTTATCTCGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAt                                                                                                                                          <  2:153378/251‑1 (MQ=255)
                                                                                                     gATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                                                    <  2:72727/248‑1 (MQ=255)
                                                                                                     gATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                                                    >  1:72727/1‑248 (MQ=255)
                                                                                                             tCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTAc                                                                                                                         >  2:20488/1‑251 (MQ=255)
                                                                                                              cTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACg                                                                                                                        >  1:374295/1‑251 (MQ=255)
                                                                                                                 aaGGCAAAGTGTATATATAACCAGCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGGATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAGTCATCCGCACATGTTACGcc                                                                                                                      <  2:508390/250‑1 (MQ=255)
                                                                                                                 aaGGCAAAGTGGATATAAAACCATCGATTAAAAAAAAGCTTTTCAACAACGCATCCGCTATTGAATAAAGAAGTTCAGGCTATAAAAAACTGGTTGAGTATCCGTACTTCGTTCACGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGGAATGCCGGGGATTATCTGGAGATTCATCCGCACATGTTACGcg                                                                                                                      <  2:217514/250‑2 (MQ=255)
                                                                                                                    gCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGTGGAAGCCGTTTTCTCGGCAACAGTTTTACCATGTTATCTCGACTTCCGGTTGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTc                                                                                                                  <  1:61118/251‑1 (MQ=255)
                                                                                                                      aaaGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgt                                                                                                                <  1:106528/251‑1 (MQ=255)
                                                                                                                        aGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTGTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtgt                                                                                                              <  1:275064/251‑1 (MQ=255)
                                                                                                                               atatatCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGACTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTg                                                                                                       >  2:342931/1‑251 (MQ=255)
                                                                                                                                 atatCCATCGAGTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCTGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGGCACTGGATAGTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGc                                                                                                      <  2:459588/250‑1 (MQ=255)
                                                                                                                                  tatCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTACGGGGAACCCGCTTTCTCGGCAACAGTTTTACCTTATTATCTCGACTTCCGGTGGTATTGCCGGGTTGTCATTGGAGATTCTTCCGCACTTGTTCCGCCATTCGTGGGGtttggttt                                                                                                    >  2:564149/1‑246 (MQ=255)
                                                                                                                                   atCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTCATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt                                                                                                   >  2:218047/1‑251 (MQ=255)
                                                                                                                                      cATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                                                                >  2:177697/1‑251 (MQ=255)
                                                                                                                                             tAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                                                         >  1:678756/1‑251 (MQ=255)
                                                                                                                                               aaaaaaGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATggg                                                                                       >  1:560093/1‑251 (MQ=255)
                                                                                                                                                 aaaaGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTGGGCGAATATGGGaa                                                                                     <  2:85660/251‑1 (MQ=255)
                                                                                                                                                  aaaGGCTTTTCAACAACGCACACGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCAGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCACCAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGTCGAATATGGGAAt                                                                                    >  1:626818/1‑251 (MQ=255)
                                                                                                                                                         tttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                                                        >  2:282044/1‑240 (MQ=255)
                                                                                                                                                         tttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                                                        <  1:282044/240‑1 (MQ=255)
                                                                                                                                                         tttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAGAGATa                                                                               <  2:351720/249‑1 (MQ=255)
                                                                                                                                                              acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAAc                                                                                                                                                                                                 >  1:307560/1‑130 (MQ=255)
                                                                                                                                                              acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt                                                                        <  1:218047/251‑1 (MQ=255)
                                                                                                                                                              acaacGCACCCGCTAGTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAAc                                                                                                                                                                                                 <  2:307560/130‑1 (MQ=255)
                                                                                                                                                                    cacacGCTATTGAATAAAGAAGTTCAGGCTGTAAAAAACTGGTTGAGTATTCGTCTTTCGTTCCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCTGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGTTTTGGCGAATAGGGGAATAGATACGCGAcgtagcca                                                                  <  2:681837/247‑8 (MQ=255)
                                                                                                                                                                           tATTGAATAAAGAAGTCCAGGATTCAAAAAACTGGTTGAGCATCAGTAATTCGTACACGCATGCTGAGAGCGAGTGGGTATTTTTCTCCCGTAAGGGGACTCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                                                           >  2:47083/1‑251 (MQ=255)
                                                                                                                                                                            aTTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                                                           <  1:200833/250‑1 (MQ=255)
                                                                                                                                                                              tgaataccgaagTTCAGGCTTCAAAAAACTGGTTGAGAATCTGTACTTCGTACCCGCATGTTGTGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTATTGCCGGGTTGACACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATTTGGGAATAGATACGCGACTTATCCAGGATTATCtt                                                        >  2:308704/9‑251 (MQ=255)
                                                                                                                                                                                  tAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGAGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                                                    >  1:504555/1‑251 (MQ=255)
                                                                                                                                                                                        agTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGc                                               >  1:349063/1‑250 (MQ=255)
                                                                                                                                                                                           tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAat                                            <  1:20488/250‑1 (MQ=255)
                                                                                                                                                                                              ggCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTc                                        <  1:591855/251‑1 (MQ=255)
                                                                                                                                                                                               gCTTTAAAAAACTGGTTGAGTATTCGTACTTCGTACCCGCATGTTGAGAGCGAGTCTGTATTTTTATCATGTAAGAGGAATCCTCTTTCTCGGCAACAGTTTTTCCATATTATCTCGACTTCCGGTGGTAATGCCGGGGTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTTGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCAtggttatattcg                                       <  1:723202/251‑12 (MQ=255)
                                                                                                                                                                                                  ttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATAACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc                                     >  2:700719/1‑250 (MQ=255)
                                                                                                                                                                                                  ttAAAAAACTGGTTGAGTATCCGTAATTCGTACCTGCATGTTGAGAGCGCGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACCGTTTTACCATATTCTCTCGCCTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCAATTTGCGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGAATTATCCAGGATTATCTTGGGCATCGCAATATTCGTc                                     >  2:703975/1‑250 (MQ=255)
                                                                                                                                                                                                   tAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGTAATAGATACGCTACTTATCCAGGATGATCTGTGTCATCGCAATATTCGTCa                                    <  2:374295/250‑1 (MQ=255)
                                                                                                                                                                                                   tAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTGTCTCGGCAACAGTTTTACCATATTAACTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa                                    <  2:626818/250‑1 (MQ=255)
                                                                                                                                                                                                            ggTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTGCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGGTTTGCTTTGGCGAATAGTGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTgg                          <  2:34244/251‑1 (MQ=255)
                                                                                                                                                                                                                          aCTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACTTAAGGGTAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGTTAATGCGGTGTTGTCACTGGAGATTCATCCGCACATGTTGCGCCTTTCGTGTGGTTTGGCTTTGGCGAATATGGGAATAGAGACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCaa            <  1:703975/251‑1 (MQ=255)
                                                                                                                                                                                                                             tCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAAt           <  1:177697/249‑1 (MQ=255)
                                                                                                                                                                                                                               gTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGc         >  2:224033/1‑249 (MQ=255)
                                                                                                                                                                                                                                    cGCATGTTGAGAGCGAGTGGGTATTGTTATCACGTTAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCTGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATAGGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTTATACTGTCTGGTATACCGCCAGCAATGCAGGGCGt  <  1:64158/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
CAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGT  >  NC_000913/4541027‑4541505

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