Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,0800GA24.1% 87.1 / 10.3 54Q89Q (CAG→CAAyraJputative fimbrial usher protein YraJ
Reads supporting (aligned to +/- strand):  ref base G (26/15);  new base A (0/13);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGTCTCCC  >  NC_000913/3288830‑3289330
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           
accaGGGACATAAACTCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATAGGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAg                                                                                                                                                                                                                                                             <  1:490013/251‑1 (MQ=255)
  caGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGaa                                                                                                                                                                                                                                                           <  1:112425/251‑1 (MQ=255)
                  tacCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATAGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACACAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAACAATATGCACAACTACCAGGAACATATCAGGATGATATCTGGCTGAATAAACAGAAGGTTTCACAGAAAAAAATTACATTTAcc                                                                                                                                                                                                                                           >  2:469352/2‑251 (MQ=255)
                  cacCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGGGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGAATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAcc                                                                                                                                                                                                                                           >  2:399430/1‑251 (MQ=255)
                                           caAACGTTTGTTGCCGAGGGGGCACACGGGCAGCGTCATGCTTTGCACTACCGCTAACGCCGAATGGTATTATTTCGACCCCAATAGGCTGAAAACAAAAAAAAGGGTTATGCAAAAAACCGATCTGTCGTGTTTATCAAAAAAATATGCACAACGCCAAGGAAGAAAACAGGTTGATATGTGGCTGAATAAAAAGAAGGTTGCACAGAAAAAAATATCAGTTACCGCCAATGCAGAGCAACTTCTGCAGc                                                                                                                                                                                                                 <  2:83003/251‑1 (MQ=255)
                                             aaCGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGcc                                                                                                                                                                                                                >  1:7317/1‑250 (MQ=255)
                                                     ttGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTAc                                                                                                                                                                                                       >  1:444357/1‑251 (MQ=255)
                                                          gATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGt                                                                                                                                                                                                    >  1:219462/1‑249 (MQ=255)
                                                                      tagCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGAAACTTATCAGGTTGATATCTGGCTGTATAAACAGCAGGATTCACAGAAAAAAATTACATTTACCGCCAATGAAAAGCAACTTCTGCAGTCACAGTTTACGGTAGAACAACTACg                                                                                                                                                                                        >  2:385044/3‑249 (MQ=255)
                                                                          aacgtcATGCTTTGCACTACCGCTAACGCCGACGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGGATGCAAACCACCGATCTGTCGCGTTTTTCAAACAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGACACAAATTACATTTACCGCGAATGCAGAGCAACTTCTGCAGCCACAGTTTAAGGAAGAAAAACTCCGTgggcg                                                                                                                                                                                  >  2:185302/6‑247 (MQ=255)
                                                                                           tACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATg                                                                                                                                                                 >  2:340979/1‑251 (MQ=255)
                                                                                            aCCGCTACCGCCGAAGAGTATTATTCCGCCCCCATTATGCTGGAACCCAAAAAAAGTGGTATGCAAACAACCGATCTGTCGCTTTTTAAAAAAAAATATGCACAACTACCAGGAACTTATCAGGTGGATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGAt                                                                                                                                                                  <  1:33927/249‑1 (MQ=255)
                                                                                            aCCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGAAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAATAAAAAGAAGTTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATg                                                                                                                                                                 <  1:142465/250‑1 (MQ=255)
                                                                                                 aaaCGCCGAAGAGTATTATTTCGACCCCATTAGGCTGAAAACCACAAAAAGTGTTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATcc                                                                                                                                                           <  1:340979/250‑1 (MQ=255)
                                                                                                    cGCCGAAGAGTATTATTTCGCCCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTCCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCgg                                                                                                                                                        <  1:170238/251‑1 (MQ=255)
                                                                                                     gCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAACAAGTGGCATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAAAAAAACTACGTGCGCTGGGAATTAAGGTGGATGAAATCCAGgc                                                                                                                                                       >  2:233079/1‑251 (MQ=255)
                                                                                                        gAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgct                                                                                                                                                    >  2:192590/1‑251 (MQ=255)
                                                                                                         aaGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATAAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACAACGGGAGCTGGGTATTAAGGTGGAGGAAATCCAGGCgctg                                                                                                                                                   >  2:71959/1‑251 (MQ=255)
                                                                                                              taatattTCGACCCCATTATGCTGGAAACCACAAAAAGTGTTATGAAAACCACCGATCTGGCGCGTTTTTAAAAAAAATATGCACAACTACCAGAATCTTATCAGGTTAATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctggctg                                                                                                                                               <  1:30611/247‑1 (MQ=255)
                                                                                                                            cATTAGGCTGTAACCCAAAAAAAGTGTTAGGCAAAAACCCGACCTGTCGCTTTTTAAAAAAAAAAATGCAAAACTACCAGGACCTTATCAGGTTGATATCTGGCTAAAAAAAAAGAAGGTTTCAAAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATa                                                                                                                                <  1:111841/251‑1 (MQ=255)
                                                                                                                             ataaggcgggAAACCAAAAAAAGTGGAAGGCAAACAACCGATTTGGCGGGTTTTTAAAAAAAATATGCCCAACTACCAGGAAATTAGCAGGTTAATATTTGGCTAAATAAAAAGAAGGTTCCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAg                                                                                                                               <  1:116168/243‑1 (MQ=255)
                                                                                                                               tatgcggaaaaCCACAAAAAGGTGTATGCAAACACCCGTTCTGTCGCTTTTTAAAAAAAAATAGGAACAACAACCAGGAACTTATCAGGTTAATATCTGGCGAAATAAAAAGAAGGTTTCAAAAAAAAAAATTACTTTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGc                                                                                                                              <  1:140466/242‑1 (MQ=255)
                                                                                                                                     ggAAACCACAAAAAGTGGTATGCAAACAACCGATCTGGCGCGTTTTTAAAAAAAATATGCACACCTACAAGGAACTTATCAGGTTGATATCTGGCTGAAAAAAAAGAAGTTTCCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCa                                                                                                                       <  1:42931/251‑1 (MQ=255)
                                                                                                                                              aaaaaGTGGTATGCAAACACCCAATCGGTCGCTGTTTTAAAAAAAAAATGCACAACTACAAGGACCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGTTTACACAAAAAAAAATTAAATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTg                                                                                                              <  1:541071/251‑1 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGa                                                                                                                                                                                                                                                            >  1:34740/1‑105 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACGTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGa                                                                                                                                                                                                                                                            <  2:34740/105‑1 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTAAAAAAAATATGCACAACTACCAGGAACTTAAAAGGTTGATATCTGGCTGAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACaa                                                                                                         <  1:71959/251‑1 (MQ=255)
                                                                                                                                                    gggTAGGCAAACACCCGACCTTCCGCTTTTTTCAAAAAAAAATGCACAACTACCAGGAACTTATCAGTTTGATATCTGGCTGAAAAAAAAGAAGTTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACaaa                                                                                                        <  1:527806/250‑1 (MQ=255)
                                                                                                                                                      gTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATc                                                                                                      <  1:457618/251‑1 (MQ=255)
                                                                                                                                                                     atatgtcgcagtttacAAAAAAAAATGAACAACTACCAAGAACGAAACAGGTTTATATCTGGATGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTAACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAgat                                                                                       <  2:397891/236‑3 (MQ=255)
                                                                                                                                                                         gTCGCGTTTTAAAAAAAAATATGCACAACTACCAGGAACTTATCAGTTTGATATCTGGCTGAATAAAAAGAAGGTTCCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAgctg                                                                                      <  1:248355/248‑1 (MQ=255)
                                                                                                                                                                               tttttCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATATGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGTTAGAACAACTACGAGAGCTGGGTATTAAGGTgg                                                                                                                                                                    >  2:419327/1‑164 (MQ=255)
                                                                                                                                                                               tttttAAAAAAAATAGGCACAACTACCAGGAACTTATCAGGTGGATATCTGGCTAAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTgg                                                                                                                                                                    <  1:419327/164‑1 (MQ=255)
                                                                                                                                                                               tttaaaaaaaaaaTATACAAAACCACCAGAAACTTATAAGGTTGAATTCTGGATAAATAAAAAGAAGGGTGCACAAAAAAAAATTCAATTACCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAAtt                                                                               <  2:219462/243‑1 (MQ=255)
                                                                                                                                                                                 tttCAAAAAAATATGCACAACTACCAGGAACTTATCAGATTGATATCTGGCTAAATAAAAAGAAGGTTTCACAGAAAAAAATTAATTTTACCGCCAATGCAGAGCAACTTATGCAGCAACAGTTTACGGTAGAACAACTACGTGAGCTGGTTATTAAGGTGGATGACATCCCGGCACTGGCTGAAAAAGATGAGGATAGCGTGCTAAACTCCCTTGAACAAATCATTCCCGGTACAGGTACGGAATTTGAt                                                                           >  2:366192/1‑251 (MQ=255)
                                                                                                                                                                                         aaaTATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAatc                                                                    >  1:422161/1‑250 (MQ=255)
                                                                                                                                                                                            tatGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGACAATGAAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAACTCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCACCTCGCTTGAACAAAGCATTACCGGTACAGCTGCTGAATTTGATTTCAatcatcc                                                                >  1:355648/1‑250 (MQ=255)
                                                                                                                                                                                             atGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAg                                                               <  1:255685/251‑1 (MQ=255)
                                                                                                                                                                                                  caACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAACAAAATTAAATTTACAGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGGGATCAACTCGCTTGAACAAATCATTCCCGGAACAGCTGCTGAATTTGATTTCAATCATCAGCAACt                                                          >  2:60321/1‑251 (MQ=255)
                                                                                                                                                                                                  aaaCTCCCAGAAATTTATCGGGTTGAAATCTGGCTGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACt                                                          <  1:60321/250‑1 (MQ=255)
                                                                                                                                                                                                    aCTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAAAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTATGCAGCAACAGTTTACGGTAGAACAACTACGTAAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTAAAAAAGATGACGATAGCGTGATCAAATCGCTTGAACAAATCATTCCCGGTACAGATGCTGAATTTGATTTCAATCATCAGCACCtt                                                         >  2:180804/1‑250 (MQ=255)
                                                                                                                                                                                                        ccaggacattacctggTTGATATCTGGATAAATAAAAAGAAGTTTTCACAAAAAAAAATTCTATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAAttt                                                    <  1:541233/237‑1 (MQ=255)
                                                                                                                                                                                                         cAGGAACTTATCAGGTTGGTATCAGGCTGAATAAAAAGAAGGTTTCACCGAAAAAAATTACACTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTTAGCTGGGTATTATGGTGGATGCAATCCCGGCGCTGGCTGAAAATGATGGCGATAGCGTGATCACATCGATTAAACAACTCATTCCCGGTACAGCTGCTGATTTTGCTTTCAATCATCAGCAACTTAGTTTg                                                   >  2:534362/1‑251 (MQ=255)
                                                                                                                                                                                                           gaaatttATCAGGTTGATATCTGGCTGAAAAAAAAGAAGGTTTCACAAAAAAAAATTACTTTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAg                                                 <  1:239282/246‑1 (MQ=255)
                                                                                                                                                                                                              aCTTATCGGGTTGATATCTGGCTGAATAAAAAAAAGGTTTCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCAt                                              <  1:349985/251‑1 (MQ=255)
                                                                                                                                                                                                                     aGGTTGATATCTGGCTGAATAATAAGAAGGTTTCCCAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAAAAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGATGAATTTGATTTCAATCAGCCGCAAATTAATTTAAGCATTCCCCa                                        >  2:387724/1‑250 (MQ=255)
                                                                                                                                                                                                                         ggATAGCTGGATGAATAAAAAGAAGGTTTCAAAGAAAAAAATTATATTTCCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAAtt                                    <  1:387724/249‑1 (MQ=255)
                                                                                                                                                                                                                               cTGGCTGAATAAAACGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGCACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGCTGACGATAGc                                                                                                                              >  2:516418/1‑154 (MQ=255)
                                                                                                                                                                                                                               cTGGCTGAATAAAAAAAAGGTTTCACAAAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGc                                                                                                                              <  1:516418/154‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTAcc                          >  1:88047/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGa                                                                                                                                     <  1:153868/141‑1 (MQ=255)
                                                                                                                                                                                                                                     gAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGa                                                                                                                                     >  2:153868/1‑141 (MQ=255)
                                                                                                                                                                                                                                       aacataagaagGTTTCACAGCAAAAAATTACATTTACCGCAAATGAAGAGCAACTTATGCAGCAAAAGTTTACGGTAGAAAACCTACGTGAGCTGTGTATTAAGGTGGATGAAATCCCGGAGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACACATCATTCCCGGTAAAGATGCTGCATTTGATTTCAATCATCTGCAACTTAATTTGAGCATTCCCCAACTTGCACTTTACCGTGAt                     >  2:571557/6‑251 (MQ=255)
                                                                                                                                                                                                                                                    ttCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCATCTTAATTTGAGCATTCCACAAATTGCACTGTACCGTGATGCAAGAGGTTACg        >  1:88892/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cacaGAAAAAGATTACATTTACCGCCTATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACCACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGAAGATAGCGTGATCAACTCGCTTGAACAAATCATTCACGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAAATTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGtc      >  2:512440/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          gAAAAACATTACATTTACCGCCAATGCAGGGCAACTTCTGCAGCCAAAGTTTACGGTAGAACAACTACGTGAGCTGGGCATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTAGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATCTGATTTCAATCATCATCAAATTAATTTGAGCATTCCCCAAATTGCACTGTAACGTACTGCAAGAGGATACGtacccc  >  2:181383/1‑246 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           
ACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGTCTCCC  >  NC_000913/3288830‑3289330

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: