Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,540,284 T→G 26.7% intergenic (‑782/‑674) nanC ← / → fimB N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,540,2840TG26.7% 69.3 / 3.1 45intergenic (‑782/‑674)nanC/fimBN‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base G (11/2);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCT  >  NC_000913/4540061‑4540502
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
gtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTACATTCTGAAGACGATCTGGGGAGGCTTGTTGATTGTGGTGttt                                                                                                                                                                                                                           >  1:63259/1‑228 (MQ=255)
      aaaaaCAAAAGATTTTAGTTTAACTATTTGTTTTAAAAATAATTTATTAAGAGTCAAAACAGGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTTTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGga                                                                                                                                                                                                                     <  1:603114/228‑1 (MQ=255)
                   cTTAGTTTAACTATTTGTTTTAAAACTAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGACGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACtt                                                                                                                                                                                                        <  2:63259/228‑1 (MQ=255)
                                        ataaataaTTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTTGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTTTTTGCCGTATCg                                                                                                                                                                                      >  2:428837/1‑225 (MQ=255)
                                        ataaataaTTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTTTTGATGTTGGTGGTTTTTGGATGAAGGTTCACTTTGTTTGTTTGcctaattg                                                                                                                                                                                      >  2:247020/1‑219 (MQ=255)
                                                          gTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTTATGCTGATGTTTCACTTAGTTTGTTTGCCGTATTGCCCTCGGAATGGCTGTGGttt                                                                                                                                                               >  2:594883/1‑227 (MQ=255)
                                                          gTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTAGATCTGGAGCGGCTTGTTTTTGTTTGTGTTTTCTGGATGCTTTTTCCCTTTGTTTGTTTTCCGTATTGCCCGGGGAAAGGTTGTGtttt                                                                                                                                                               >  1:117409/1‑227 (MQ=255)
                                                                     gggAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTTAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTTGTGTTGGGGAAGGtttt                                                                                                                                                   >  2:253898/1‑227 (MQ=255)
                                                                     gggAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCACATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTCAGGATGATGGTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGGATGGCTGTGCTTGTGGAAGGgtt                                                                                                                                                    >  2:267328/1‑227 (MQ=255)
                                                                         gCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCAATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGCGAGGCTTGTTGATGTTGGGGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTTTCGACCGGCGAATGGTTGGGGTTGAGGAAAGGTTAATTc                                                                                                                                               >  2:197909/1‑228 (MQ=255)
                                                                            ttGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTa                                                                                                                                            <  2:448735/228‑1 (MQ=255)
                                                                             tGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGTTGATGTTTTACTTAGCTTGTTTGCCGTCTAGCCCGGCGGATGGGTGTGATTGGGGAACGTTTAAGTCGTAg                                                                                                                                           >  2:502513/1‑228 (MQ=255)
                                                                                    tAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATACTTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTTTTGTTGTGGTTTTCTGGGTGGTGTTTCACTTAGGTTGTTTGCCTGTTGGCCCGGTGCATGGCGGTCTTGTGTGAAGGTTTTATTCGTTGTGAccca                                                                                                                                    >  1:266421/1‑226 (MQ=255)
                                                                                            tGAGCTTGGTGAAAATCAGAGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGATTGTGGCTTGTAACTGTTCACTTCTGAAGTCGATCTGTAGAGTCTTGTTTTGGTTGGTGTTTTCAGGATGATTGTTCACTTTGTTTTTTTGCGGTATTGGCGCGGGAACGGGTGTGTTTAGGGGATGTTTTATGCGTACAGGCCCAAGCTatat                                                                                                                            >  1:633654/1‑228 (MQ=255)
                                                                                                cTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCCCTTCTGCAGTCGATCTGGAGAGGCTTGTTGATGTTTGGGTTTTTCGGGTGGTGTTTTACTTAGGTTGTTTGGCGTATTGCCCGGGGAATGGGCTGGGTTTGGGGAGGTTTTTGTCGTATTGACCAAAGGTATATTTAc                                                                                                                        >  1:570440/1‑228 (MQ=255)
                                                                                                cTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTTAGGAAGGGTTAAGTCGTAATGACCAAAGCCTTATTTAc                                                                                                                        >  2:199046/1‑228 (MQ=255)
                                                                                                             aGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTTACGTTGGGGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGATTTGATTGTGGAAGGGTTTAATCGTTTTAGCCCAGACCATATTTACCAAccgattgaga                                                                                                           >  2:679431/1‑219 (MQ=255)
                                                                                                                         ccATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGGTGATGTTTCACTTAGTTTGTTTGCCGTATGTCCCGGCGAATTGCTGTGTTTGTGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTGGATGAAAAAatca                                                                                                >  2:656349/1‑227 (MQ=255)
                                                                                                                               gggtGAGCATAATTCTCATCATGGAATATGTGTCCTGGTTTGTGGATTTTAACTGGGCACTTGGGAAGGTGAGCTGGAGAGGCGTGGTGATGTTGGTGTTTGCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTg                                                                                         <  1:502513/225‑1 (MQ=255)
                                                                                                                                    aGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGTGAGGCTTGTTTATGTTGGGGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGTTTGAGGAAGGTTTTAGTCGTATTGGCCAAATCTATATTTACCAACGAATTTAAATGAAAAAAACATTTCCTGCGt                                                                                      >  2:105185/1‑226 (MQ=255)
                                                                                                                                    aGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAg                                                                                                                                                 <  1:206262/167‑1 (MQ=255)
                                                                                                                                    aGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAg                                                                                                                                                 >  2:206262/1‑167 (MQ=255)
                                                                                                                                      cATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTTACTTAGTTTGTTTGCCGTATTGCCCGGCGAATGGTTGTGATTTAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGGAAAAATCGTCTCCTGCGTTccc                                                                                  >  2:645835/1‑228 (MQ=255)
                                                                                                                                          tttCTCATTGTGAAATATGTTTCCTGGTTTGTGGCTTGTTACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCat                                                                              <  1:199046/227‑1 (MQ=255)
                                                                                                                                          gaggtgtataagggacaggTTTCCTGGTTTGTGTCTTGTAACTGGTAACTGCTGAAGTCGATCTGGAGAGGCTTGGTGGTGTTGGGGTTTTCAGGATGATGTTTAACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTTACCAAAGCTATATTTACCAACGAATGTAGAGGAAAAAATCATCTCCTGCGTTCCCCCat                                                                              <  2:73554/210‑1 (MQ=255)
                                                                                                                                                atcatGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTTGGGAAGGTTTTAGTCGTATGGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATctct                                                                        >  2:77730/1‑228 (MQ=255)
                                                                                                                                                 tcatGAAATATGTTGCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAatc                                                                                                 <  2:160027/202‑1 (MQ=255)
                                                                                                                                                   gtGAAATTTGTTTCCTGGTTTGTGGCTGGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTGGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTa                                                                       <  2:307316/225‑1 (MQ=255)
                                                                                                                                                   atGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGGGTTTTCCGGGTGGTGTTTCCCTTCGTTTGTTTGCCCTATCTCTCGGGCAATAGCTGTGGTTGGTGGAGGTTTATGTCGTATTGACCAAAGCTATATTTAACAAACAATTTATTTGAAAAAATACCTTCCTGAGTTCCCCCATAtcttcatg                                                                     >  1:561250/1‑223 (MQ=255)
                                                                                                                                                         tatGTTTGTGGGTTTGGTGTTTTTAACTGGTCACTTCGGAAGTTGATCTGGAGAGGCTTTTTGAGTTTGGTGTTTTCAGGATTATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATaaaa                                                               <  1:645835/228‑1 (MQ=255)
                                                                                                                                                                 cTGGTTTGTTGCTTGTAACTGGTCACTTCGGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATaaaa                                                               <  2:298877/220‑1 (MQ=255)
                                                                                                                                                                               gTAACTGGTCACTTCTGAAGTCGATCTGGATAGGCTTTGGTATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTAACGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGc                                         <  2:297186/228‑1 (MQ=255)
                                                                                                                                                                                       tCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCt                                  <  2:117409/227‑1 (MQ=255)
                                                                                                                                                                                                     gATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTTAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGTTGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAACTCACCAAAcacacgaa                   >  1:95627/1‑222 (MQ=255)
                                                                                                                                                                                                     gATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGGAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCg                                                                                       >  2:125380/1‑160 (MQ=255)
                                                                                                                                                                                                      aTCTGGAGAGGGTTGTTGATGTTGGTGTTTGCAGGATGATGGTTCACTTAGTTTGTTTGCCGAATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAGAATCATCTCCTGCGTTCCCCCATATCTCTAGGATaaaa                                                               <  2:338345/183‑1 (MQ=255)
                                                                                                                                                                                                      aTCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGtt                                                                                     <  2:657159/161‑1 (MQ=255)
                                                                                                                                                                                                      aTCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGGTTCACTTAGTTTGTTTGCCGGATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGt                                                                                      <  1:125380/160‑1 (MQ=255)
                                                                                                                                                                                                       tCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGaat                 <  2:95627/228‑1 (MQ=255)
                                                                                                                                                                                                               ggCTTGTTGATGTTGGTGTTTTCAGGATGATGGTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTTAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAACATCATCTCCTGCGTTTCCCCATATCTCTAGGATAAAAAGGAAAGTAAAAATATCATTGCGTAAGCTGTCGAATAAGCAGGAAAAATcacaa         >  1:698140/1‑224 (MQ=255)
                                                                                                                                                                                                                gCTTGTTGATGTTGGTGTTTTCAGGATGATGGTTCACTTAGTTTGTTTGCCGTATCGCGCGGCGAATGGCTGGGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATTATCTCCAGCGTTCCCCCATATCTCTAGGAATAAAAGGAAATTAACACTCCCATTGCGTAAGCTGTCGAATCAACATGGATAATCGCAAg        >  1:691955/1‑228 (MQ=255)
                                                                                                                                                                                                                gCTTGTTGATGTTGGGGTTTTCAGGATGATGTTTCACGTAGGTTGTTTGCCGTATCGCCCGTCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTTACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAg        <  1:428837/228‑1 (MQ=255)
                                                                                                                                                                                                                gCTTGGTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTATTGACCAAAGCTATATTTACCAACGACTGTAGGTGAAAAAATCATATCCTACGTTGCCCCATATCTCTAGGATAAAAAGGAATGTAAAAATCTCATTGCGTAAGCTGCCGAATcaaaagggaaaaaagccag        >  1:339851/1‑211 (MQ=255)
                                                                                                                                                                                                                       gATGTGTGTGTTTTCAGGGTTATTTTTCACTTAGTTTGTGTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCCTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCt  <  1:644394/227‑1 (MQ=255)
                                                                                                                                                                                                                          gTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAAAAATTTCATGGCGTAAGCTGACGACTAAGCAGGAATAATCGCAAGGGAcc   >  2:644394/1‑223 (MQ=255)
                                                                                                                                                                                                                              gtgtTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGc                                                                                        >  2:594825/1‑134 (MQ=255)
                                                                                                                                                                                                                               tgtTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCg                                                                                       <  1:594825/134‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
GTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCT  >  NC_000913/4540061‑4540502

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