Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,364,405 C→A 11.6% S534* (TCA→TAA)  yhcD → putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,364,4050CA11.6% 237.8 / 2.0 83S534* (TCA→TAA) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base C (61/11);  new base A (2/9);  total (63/20)
Fisher's exact test for biased strand distribution p-value = 2.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTTGGTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTG  >  NC_000913/3364176‑3364641
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
cTTGGTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCg                                                                                                                                                                                                                                      >  2:1016715/1‑238 (MQ=255)
  tGGGTTAAACCGTGTCTAACCCAATTAGGGGGGGTTTCAACATCAGTGGCTAAAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTAAATATAAACGGGAACAAAGCTTTCAAACCCGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTAAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGACCAACCACCGcc                                                                                                                                                                                                                                    <  2:1090284/239‑1 (MQ=255)
   gggttaaaactgtacaaccccaaTTCGGGGTGGTTTCAACACCAGGGACTAAAAAACATAGAGTTCCGCGGGTGAAAATGCCTAAAGCAACTTTAATATTAACTGGAACAAGGCTTCAAACCACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGCCGATGGTAATCTTTTCTACGTCAACATCAGTATTCCATTTGGAAGATCAAACACCGCca                                                                                                                                                                                                                                   <  2:783361/217‑1 (MQ=255)
   ggTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGCAAATAATGAAGACGATGGTGATCTGTTCTACGTCACCATCAGTATTCCATTTGGCAGATCAAACACCGCca                                                                                                                                                                                                                                   >  2:990134/1‑238 (MQ=255)
   ggTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGGGCTTCAGGACATAATGAAGACGATGCTGATCTGTTCAACGTAAACATCAGTATTCCATATGGAAAATCAAACACCGCca                                                                                                                                                                                                                                   >  2:990133/1‑238 (MQ=255)
    gTTTAAACTTGTCAAACCTCATTGTGGGTGTTTCCAACATCAGTGGCTAAAAAAAATAATGTTACGAAGGTGACAATGACCAAACAAACTTAAATATTAACTGAAAAAAAGCGTTAAAACACGCCACGGTTCCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCcaca                                                                                                                                                                                                                                 <  2:471542/239‑1 (MQ=255)
               gTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGATCAAACACTCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGCATTCCAATTGTCAAATCAAAAACACCCACACTgaaaaca                                                                                                                                                                                                                        >  1:752924/1‑232 (MQ=255)
               gTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAATCCAGCCACCTTAATATTACCTGGAACAACGCGTTCAAACACGCAACGGTTTCCGCCAACTGGCTGCATCAACTTAGCGCTTCAGAAAATAATGACGACGCTGGGGACCTGTTCTACGTCAAATTCAGTATTCACTTTGCCAAGCCAAACACACCcaccagattaaaa                                                                                                                                                                                                                        >  1:1195189/1‑228 (MQ=255)
                          ttCTGGGTGGTTTCAACATCAGTGGCTATAACACATATAGTTACGACGTTGATCATGACTTAAGCATCCTTAATATTAATTGGAACAAAGCGTTCAAACACGCCCCGGTTTCCGCTAACAGGATGCATCTACTTAGTGCTTCAGAAAATAATAAACACGCTGGTTCTCTGTCTTACACACACCTCAGGATTCCATTTCGAAGATCAAAAACCTCCACACTgataaccccccataacgaa                                                                                                                                                                                                            >  1:139280/1‑221 (MQ=255)
                          ttCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAAGGACTCAAGCAACCTTAATATTAACTGGAACAAACCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACATAAAAATCAGAATTCCATTTGGCAAATCAAAAACCGCCACACTgaaaaacaccaatacaaaa                                                                                                                                                                                                            >  1:882427/1‑221 (MQ=255)
                               ggtggtTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGCTTCCGTCAACAGGCAGCATCAACTTATTGCTTCAGAAAATAATGAAGACGATGGTGAACTGTTCTACGTAAACATCAGTATTCCATTTGGCAGATCGAACACCGCCACACTGAATACTCGCCATGAcaaacaaaa                                                                                                                                                                                                       >  1:318440/1‑231 (MQ=255)
                                 gttttcaaaCTTCAGTGGCTAAAAAAAGTATAGTTACGACGGTAAAAATGCCTCAATCAACTTAAATATTAACTGAAACAAGCCGTTAAAACACGCAGGGGTTTCCTTAACCGGGCAGCATCAACTGTGGCCTATAAAAAATAATAAATACGATGGGAATTTTTTCTCCGTCAACAGAAGTATCCTTTTTGGCAGAAAAAACACCGCCACACTGTATACTCGCCATGACGATCATaaa                                                                                                                                                                                                      <  2:318440/231‑1 (MQ=255)
                                    tttCAACATCAGGGGTTAAAAAACATATAGTTACGAGGGTGACAATGACTCAAGCAACTTTAATATGAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGAAGCATAAATTTAGTGTTAAAAAAAATAATGAAGCCGATGGTGATCTGTTCTACGAAAACATCAGTATTCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAcc                                                                                                                                                                                                   <  2:159566/238‑1 (MQ=255)
                                      taaaataaagggccTAAAAAAAATATATTAACGAGGGGAAAAAAGACTAAAAAACCTTAAAATTAACTGGAAAAAAGGTTTGCACACCCCCCGAGTTTCCCGTAACCTGGCACCAACAACTTGATGCTTAAAAAAAAAATGAAGACGAGGGTAATCTGTTTTCCGTAAACACCGTTTTTCCTTTTGGCAGATAAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCAc                                                                                                                                                                                                <  1:942768/226‑1 (MQ=255)
                                        aaCATCGGAGGATAAAAAACATAATGTAAGGGCGGTAAAAATCACTAAAGCAACGTTATTATTAACGGAAACAAGGTGTAAAAACCCGCCACGTTTTCCGTAAACTGGCACCATCATTTTAGTGCTTCAGAAAAAAATGAAGACGATTGTGATTTGTTCTACATCATAAACATTATTCCATTTGGAAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTa                                                                                                                                                                                              <  2:139280/239‑1 (MQ=255)
                                              aGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTCTGGAAc                                                                                                                                                                                        >  2:1024377/1‑239 (MQ=255)
                                              aGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGAATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                                                                                                                                                                                        >  2:326181/1‑239 (MQ=255)
                                                tGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATAATAAAACCCACTATGGAACTg                                                                                                                                                                                      >  2:221919/1‑239 (MQ=255)
                                                  gCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgt                                                                                                                                                                                    >  2:919932/1‑239 (MQ=255)
                                                              ataGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACTCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCATCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACACCAGCATTCAATTTGGCAGATAAACCACCGCCAAACTGAATACTCCCCATAACCATCCTAAAACCAACAATggcaatgcacaacaagaaaag                                                                                                                                                                        >  1:1256543/1‑220 (MQ=255)
                                                              ataGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc                                                                                                                                                                        >  2:405453/1‑239 (MQ=255)
                                                              ataGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc                                                                                                                                                                        >  2:709747/1‑239 (MQ=255)
                                                              ataGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGATTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATTACGACCATAAAACCCACTATGGAATCGGGGTAATGGGAg                                                                                                                                                                          >  2:1218728/1‑237 (MQ=255)
                                                              ataGTTACGACGGGGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGCCGATGGTGATCTGTTCTACTTCAACATCAGTATTCCATTTGGCAGATCAAACACCGACCCACTGGATACTCGCCATGACGATCATAAAACCCACCATGGAACAGGGTTCAAggggaag                                                                                                                                                                        >  1:1116688/1‑235 (MQ=255)
                                                                    acgacgGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAGATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGCCATGGGACTCGTCTCa                                                                                                                                                                  >  2:447407/1‑239 (MQ=255)
                                                                     cgacgGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCAAACTGTATACTCGCCATCACGACAATAAAACCCCCAAAGGAACACGAGCCATGGGAgacggtccaa                                                                                                                                                                 >  1:637010/1‑233 (MQ=255)
                                                                        cgGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtata                                                                                                                                                                                                                          >  2:1091597/1‑179 (MQ=255)
                                                                          gTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                            >  2:542467/1‑239 (MQ=255)
                                                                          gTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATAATAAAACCCACTATGGCACCGGTGTCATGGGAatcggcccgacagca                                                                                                                                                            >  1:367794/1‑228 (MQ=255)
                                                                           tgaaaatacccaaaacAACCTGAATATGACTGGAAACAAAGCGTAGAAACCAGCCACGCGTTCCCACAAGGGGCACCTCAAACTTAGTGATTCAAAAAAAAAGAGAGCGAATGGTCATTTGGTCTCAGACAAAATAAATATTCCATTTGAAAGAACAGACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                            <  1:1186115/223‑1 (MQ=255)
                                                                           taaCAATGACCAAAGCACCCTTAATATTAACGGAAAAAAACCTTTAAAACCCGCACGGTTTTCCGTCACCTGGCAGCAACACCTTAGTCCTTAAAAAAATAATGAAGACGTTGGTGATCTGTTCTACGTCAACACCAGTATTCTTTTTGGAAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                            <  2:860923/236‑1 (MQ=255)
                                                                            ttaactgaccccacacatcTTTATTTTAAATGGAAAAAAGTGTTCAAACACCGCCAGTGTTTCCGAAAATTGGCACATTCAACTTAGTCTTTAAAAAATAAATAAAGGCGATGGTATCTTTTTCTACTTAAACATCATTATTCCTTTTGGAAAACCAACCCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                            <  1:608411/220‑1 (MQ=255)
                                                                            gACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgagat                                                                                                                                                          >  2:216473/1‑239 (MQ=255)
                                                                                         gCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAtgtg                                                                                                                                             >  2:908895/1‑239 (MQ=255)
                                                                                             ccTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCTTAAAACCCACTATGGTACTGGAGTCATGGGGGTCGTCTCAGATGAGATTTCCTACTATGTGAATg                                                                                                                                         >  2:53383/1‑239 (MQ=255)
                                                                                             ccTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACAATGTGAAt                                                                                                                                          >  2:642375/1‑238 (MQ=255)
                                                                                                  aaTATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGa                                                                                                                                     >  2:22938/1‑239 (MQ=255)
                                                                                                   atatTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTCCTGGTGTCATGGGAGTCGTCTCCGATGGAATGTCCTACTATGTGAATGCTGa                                                                                                                                     >  1:58808/1‑238 (MQ=255)
                                                                                                   atatTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATACTAAAACCCACTCTGGTACTGGGGTCATGGGGTACGTCTCAGACGAGAAGTCcaactcatgaaattcttaa                                                                                                                                    >  1:402782/1‑225 (MQ=255)
                                                                                                    attaaaCGGGAACAAACCTTTCAAACCCGCCACGGTTCCCACCAACTGGAAGCATAAATTTAGTCCTTAAAAAAAAAAGGAAGACGTGTGGACTCTTTTCTCCGTCAACAATATTTTTCCATTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAAc                                                                                                                                   <  2:58808/235‑1 (MQ=255)
                                                                                                        aaCTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCAt                                                                                                                                                                               >  2:235906/1‑191 (MQ=255)
                                                                                                         aCTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACgaga                                                                                                                               >  2:426017/1‑238 (MQ=255)
                                                                                                         aCTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCACCTATGGAACTGGGGTCATGGGAGTCGTCTCAGATAAGATGTCCAACTAATTGAATGCTGAAGGAAAt                                                                                                                              >  1:395402/1‑239 (MQ=255)
                                                                                                         aCTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGCCAACTGGCAGCATCCACTTAGCGCTTCAGAAAATAATCACGCCGATGTTGATCTGTTCTCCCTCAACATCAGAATTCCATTTGGCAGATCACACACCTCCACACTGTATACCCGCCATGAAGATCCTCAAAACCACTATGCGACTGGAGCCATAGAAGCCGCACCAGGTGAGCTGTCCAACACTGTGAATGCCGAACgaga                                                                                                                               >  1:524255/1‑238 (MQ=255)
                                                                                                               aaCAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacg                                                                                                                          >  2:1304065/1‑238 (MQ=255)
                                                                                                                     acgggcgtcaaaaccccccagtttttCCCTCAACGGGAAACAAAAACTTGGGGCTTAAAAAAAAAATGAAGACGAGGGTGACTTTTTCTCCGTCAACATCAGTATTCCTTTTGGCAAATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacgacg                                                                                                                       <  2:479655/216‑1 (MQ=255)
                                                                                                                         tCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTa                                                                                                                                                                                         >  2:806495/1‑165 (MQ=255)
                                                                                                                         tCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGTCATGGGAGTCGGCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACAACGAACGGGaa                                                                                                               >  2:693840/1‑239 (MQ=255)
                                                                                                                          cAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGaa                                                                                                               >  2:93530/1‑238 (MQ=255)
                                                                                                                                      tgggTTTCCGAAAACTGCAACCATCAACTTTGTCCTTAAAAAAATAATGAAGACGGGGGTGATCTGTTCTACGTCAACATCATTATTCCATTTGAAAGAAAAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGaa                                                                                                    <  2:882427/237‑1 (MQ=255)
                                                                                                                                                  gtaaaCTGGCAGCATCAACTTAGGGCTTCAGAAAATAATGAAGGCGATGGTGATCTGTTCTACGTCAACATCAGTATCCTTTTGGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTAt                                                                                           <  2:395402/236‑1 (MQ=255)
                                                                                                                                                       cTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGt                                                                                       >  2:331092/1‑238 (MQ=255)
                                                                                                                                                               aTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGGGAAACGGGCTTGAAAGGcgagacaagttcacatccc                                                                              >  1:486888/1‑221 (MQ=255)
                                                                                                                                                                        gttctaaaaaaaaTAAGGAAGACGATTTGTTCCTTTTCTACGTAAACAACAGTTTTCATTTTGACATATAAAACCCCGCCACACTGTTCACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACAcc                                                                     <  1:542467/230‑1 (MQ=255)
                                                                                                                                                                        gTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTacac                                                                      >  2:803293/1‑238 (MQ=255)
                                                                                                                                                                         tGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGATTCGTCTCAGATGAGATGTCCTACTATGGGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTAGAACGGCAGTAACAATTCCAAACCCCATTaaaacc                                                                    >  1:67031/1‑234 (MQ=255)
                                                                                                                                                                              caaaaatTAATAAAGCCGGTGGTGACTGTTCTCACGTCAAAATCAGCATTCCATTTGGCAGTACAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTc                                                               <  2:67031/236‑1 (MQ=255)
                                                                                                                                                                              cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAgcc                                                               >  2:423918/1‑237 (MQ=255)
                                                                                                                                                                              cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTc                                                               >  2:1141075/1‑239 (MQ=255)
                                                                                                                                                                              cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTAATGGGAGTCGTCTCATATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGCCGAACGGGAAACAAGCTTGCACGGCAGGAACAGATCCAATTGCCAATAAACCCaaaga                                                               >  1:633374/1‑236 (MQ=255)
                                                                                                                                                                              cAGAAAATAAGGAAGACGATCGTGATCTGTTCTACGGCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTTGTGTCATGGCAGTCGTCTCAGATTAGATGTCCTACTATGTGAATGCTGAACGAGATACCGACGAACGGAAAACGGACTTGACCGTCACGATCAATACCAATCTCCATTAAACACAAgaa                                                               >  1:1002266/1‑237 (MQ=255)
                                                                                                                                                                                     taatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGTCATGGGAGTCGTCTCAGATGAGCAGACCTACTATGAGAATGCTGAACGAGATCACGACTAACGTGAAACAGCCTAGAACAGCGATACCCATTCCAATCTCCAGAACACCAAAGACAgaatag                                                        >  1:361669/1‑234 (MQ=255)
                                                                                                                                                                                                  tGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGGGAAACGAGCTTGAACGGCATAATCAGTTCCAACCCCCATTCAAACCAAATCAGCCTTGCCGCAgagacaa                                            >  1:500012/1‑233 (MQ=255)
                                                                                                                                                                                                  tGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGAACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGCGAATGCTGAACGAGATCACGACGACCGGACAACGAGCTTGAACGGAAGTATCAGTTCCAATCTCCATTACAacccgtctcttatacacatttcccac                                             >  2:753434/1‑215 (MQ=255)
                                                                                                                                                                                                         tCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGACCGGCAGTATCAGTTCCAACCTCCATTACACCCAAGTCAGCCTTGCCGCCAGAGAAACCGGCa                                      >  1:529087/1‑239 (MQ=255)
                                                                                                                                                                                                           cTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAg                                     >  2:1157862/1‑239 (MQ=255)
                                                                                                                                                                                                            ttggtttcacataaacaacaacgttCCATTTGGACAAGCAAAGCCCCGGCACACAGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCgg                                        <  2:1223657/215‑1 (MQ=255)
                                                                                                                                                                                                               ttCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGGGCAACGGGCAGGGa                                  >  2:33024/1‑239 (MQ=255)
                                                                                                                                                                                                               ttCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTAAACCCAAGGCAGCCTTGCCGCAGGAGCAAGCGGCAgcaa                                  >  2:1312014/1‑236 (MQ=255)
                                                                                                                                                                                                               ttCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAAAGGCAGTATACGTTCCAATCTCCATTACACCCAAATCACCCTTGCCGCCGGAACAAGCGGCAGGGa                                  >  1:961965/1‑239 (MQ=255)
                                                                                                                                                                                                                     cctaaacatatgtagtctttttGGCAGAAAAAACCCCGCCACACTGTATACCCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGcc                             <  2:402782/220‑1 (MQ=255)
                                                                                                                                                                                                                     cGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGcc                             >  2:70320/1‑239 (MQ=255)
                                                                                                                                                                                                                          aCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAAc                                                                                                             >  2:562226/1‑154 (MQ=255)
                                                                                                                                                                                                                          aCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACt                        >  2:40738/1‑239 (MQ=255)
                                                                                                                                                                                                                          aCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTAAACCCAAGCCAGCCTTTCCCCAGGAGCAAGCGGCatggtaaaccgaact                        >  2:609602/1‑225 (MQ=255)
                                                                                                                                                                                                                                 tATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGa                                                >  2:247053/1‑208 (MQ=255)
                                                                                                                                                                                                                                      tttttGGCAAAACAAACCCCGCCACACTTTATAACCGCCATGAAGGTCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACg            <  1:70320/237‑1 (MQ=255)
                                                                                                                                                                                                                                      cATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGCACTTACAACGGTCCg            >  2:1276/1‑239 (MQ=255)
                                                                                                                                                                                                                                        tttGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGAt          >  2:54587/1‑239 (MQ=255)
                                                                                                                                                                                                                                         ttGGCAGATCAAACACCGCCACACTGTATGCTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAg                                               >  2:446694/1‑201 (MQ=255)
                                                                                                                                                                                                                                           ggCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctac                                                                                                                                                   >  2:393597/1‑99 (MQ=255)
                                                                                                                                                                                                                                           gccAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctac                                                                                                                                                   <  1:393597/97‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAgggg  >  2:1065919/1‑237 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
CTTGGTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTG  >  NC_000913/3364176‑3364641

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