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→CAA) yraJputative 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.

ACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGTCTCCC  >  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)
                                           caAACGTTTGTTGCCGAGGGGGCACACGGGCAGCGTCATGCTTTGCACTACCGCTAACGCCGAATGGTATTATTTCGACCCCAATAGGCTGAAAACAAAAAAAAGGGTTATGCAAAAAACCGATCTGTCGTGTTTATCAAAAAAATATGCACAACGCCAAGGAAGAAAACAGGTTGATATGTGGCTGAATAAAAAGAAGGTTGCACAGAAAAAAATATCAGTTACC‑GCCAATGCAGAGCAACTTCTGCAGc                                                                                                                                                                                                                 <  2:83003/251‑1 (MQ=255)
                                             aaCGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGcc                                                                                                                                                                                                                >  1:7317/1‑250 (MQ=255)
                                                     ttGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTAc                                                                                                                                                                                                       >  1:444357/1‑251 (MQ=255)
                                                          gATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGt                                                                                                                                                                                                    >  1:219462/1‑249 (MQ=255)
                                                                      tagCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGAAACTTATCAGGTTGATATCTGGCTGTATAAACAGCAGGATTCACAGAAAAAAATTACATTTACC‑GCCAATGAAAAGCAACTTCTGCAGTCACAGTTTACGGTAGAACAACTACg                                                                                                                                                                                        >  2:385044/3‑249 (MQ=255)
                                                                          aacgtcATGCTTTGCACTACCGCTAACGCCGACGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGGATGCAAACCACCGATCTGTCGCGTTTTTCAAACAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGACACAAATTACATTTACC‑GCGAATGCAGAGCAACTTCTGCAGCCACAGTTTAAGGAAGAAAAACTCCGTgggcg                                                                                                                                                                                  >  2:185302/6‑247 (MQ=255)
                                                                                           tACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATg                                                                                                                                                                 >  2:340979/1‑251 (MQ=255)
                                                                                            aCCGCTACCGCCGAAGAGTATTATTCCGCCCCCATTATGCTGGAACCCAAAAAAAGTGGTATGCAAACAACCGATCTGTCGCTTTTTAAAAAAAAATATGCACAACTACCAGGAACTTATCAGGTGGATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGAt                                                                                                                                                                  <  1:33927/249‑1 (MQ=255)
                                                                                            aCCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGAAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAATAAAAAGAAGTTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATg                                                                                                                                                                 <  1:142465/250‑1 (MQ=255)
                                                                                                 aaaCGCCGAAGAGTATTATTTCGACCCCATTAGGCTGAAAACCACAAAAAGTGTTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATcc                                                                                                                                                           <  1:340979/250‑1 (MQ=255)
                                                                                                    cGCCGAAGAGTATTATTTCGCCCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTCCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCgg                                                                                                                                                        <  1:170238/251‑1 (MQ=255)
                                                                                                     gCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAACAAGTGGCATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAAAAAAACTACGTGCGCTGGGAATTAAGGTGGATGAAATCCAGgc                                                                                                                                                       >  2:233079/1‑251 (MQ=255)
                                                                                                        gAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgct                                                                                                                                                    >  2:192590/1‑251 (MQ=255)
                                                                                                         aaGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATAAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACAACGGGAGCTGGGTATTAAGGTGGAGGAAATCCAGGCgctg                                                                                                                                                   >  2:71959/1‑251 (MQ=255)
                                                                                                              taatattTCGACCCCATTATGCTGGAAACCACAAAAAGTGTTATGAAAACCACCGATCTGGCGCGTTTTTAAAAAAAATATGCACAACTACCAGAATCTTATCAGGTTAATATCTGGCTAAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCgctggctg                                                                                                                                               <  1:30611/247‑1 (MQ=255)
                                                                                                                            cATTAGGCTGTAACCCAAAAAAAGTGTTAGGCAAAAACCCGACCTGTCGCTTTTTAAAAAAAAAAATGCAAAACTACCAGGACCTTATCAGGTTGATATCTGGCTAAAAAAAAAGAAGGTTTCAAAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATa                                                                                                                                <  1:111841/251‑1 (MQ=255)
                                                                                                                             ataaggcgggAAACCAAAAAAAGTGGAAGGCAAACAACCGATTTGGCGGGTTTTTAAAAAAAATATGCCCAACTACCAGGAAATTAGCAGGTTAATATTTGGCTAAATAAAAAGAAGGTTCCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAg                                                                                                                               <  1:116168/243‑1 (MQ=255)
                                                                                                                               tatgcggaaaaCCACAAAAAGGTGTATGCAAACACCCGTTCTGTCGCTTTTTAAAAAAAAATAGGAACAACAACCAGGAACTTATCAGGTTAATATCTGGCGAAATAAAAAGAAGGTTTCAAAAAAAAAAATTACTTTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGc                                                                                                                              <  1:140466/242‑1 (MQ=255)
                                                                                                                                     ggAAACCACAAAAAGTGGTATGCAAACAACCGATCTGGCGCGTTTTTAAAAAAAATATGCACACCTACAAGGAACTTATCAGGTTGATATCTGGCTGAAAAAAAAGAAGTTTCCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCa                                                                                                                       <  1:42931/251‑1 (MQ=255)
                                                                                                                                              aaaaaGTGGTATGCAAACACCCAATCGGTCGCTGTTTTAAAAAAAAAATGCACAACTACAAGGACCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGTTTACACAAAAAAAAATTAAATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTg                                                                                                              <  1:541071/251‑1 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGa                                                                                                                                                                                                                                                            >  1:34740/1‑105 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACGTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGa                                                                                                                                                                                                                                                            <  2:34740/105‑1 (MQ=255)
                                                                                                                                                   gTGGTATGCAAACAACCGATCTGTCGCGTTTTTAAAAAAAATATGCACAACTACCAGGAACTTAAAAGGTTGATATCTGGCTGAAAAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACaa                                                                                                         <  1:71959/251‑1 (MQ=255)
                                                                                                                                                    gggTAGGCAAACACCCGACCTTCCGCTTTTTTCAAAAAAAAATGCACAACTACCAGGAACTTATCAGTTTGATATCTGGCTGAAAAAAAAGAAGTTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACaaa                                                                                                        <  1:527806/250‑1 (MQ=255)
                                                                                                                                                      gTATGAAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATc                                                                                                      <  1:457618/251‑1 (MQ=255)
                                                                                                                                                                     atatgtcgcagtttacAAAAAAAAATGAACAACTACCAAGAACGAAACAGGTTTATATCTGGATGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTAACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAgat                                                                                       <  2:397891/236‑3 (MQ=255)
                                                                                                                                                                         gTCGCGTTTTAAAAAAAAATATGCACAACTACCAGGAACTTATCAGTTTGATATCTGGCTGAATAAAAAGAAGGTTCCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAgctg                                                                                      <  1:248355/248‑1 (MQ=255)
                                                                                                                                                                               tttttCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATATGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGTTAGAACAACTACGAGAGCTGGGTATTAAGGTgg                                                                                                                                                                    >  2:419327/1‑164 (MQ=255)
                                                                                                                                                                               tttttAAAAAAAATAGGCACAACTACCAGGAACTTATCAGGTGGATATCTGGCTAAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTgg                                                                                                                                                                    <  1:419327/164‑1 (MQ=255)
                                                                                                                                                                               tttaaaaaaaaaaTATACAAAACCACCAGAAACTTATAAGGTTGAATTCTGGATAAATAAAAAGAAGGGTGCACAAAAAAAAATT‑CAATTACCCGCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAAtt                                                                               <  2:219462/243‑1 (MQ=255)
                                                                                                                                                                                 tttCAAAAAAATATGCACAACTACCAGGAACTTATCAGATTGATATCTGGCTAAATAAAAAGAAGGTTTCACAGAAAAAAATTAATTTTACC‑GCCAATGCAGAGCAACTTATGCAGCAACAGTTTACGGTAGAACAACTACGTGAGCTGGTTATTAAGGTGGATGACATCCCGGCACTGGCTGAAAAAGATGAGGATAGCGTGCTAAACTCCCTTGAACAAATCATTCCCGGTACAGGTACGGAATTTGAt                                                                           >  2:366192/1‑251 (MQ=255)
                                                                                                                                                                                         aaaTATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAatc                                                                    >  1:422161/1‑250 (MQ=255)
                                                                                                                                                                                            tatGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GACAATGAAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAACTCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCACCTCGCTTGAACAAAGCATTACCGGTACAGCTGCTGAATTTGATTTCAatcatcc                                                                >  1:355648/1‑250 (MQ=255)
                                                                                                                                                                                             atGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTAAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAg                                                               <  1:255685/251‑1 (MQ=255)
                                                                                                                                                                                                  caACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAACAAAATTAAATTTACA‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGGGATCAACTCGCTTGAACAAATCATTCCCGGAACAGCTGCTGAATTTGATTTCAATCATCAGCAACt                                                          >  2:60321/1‑251 (MQ=255)
                                                                                                                                                                                                  aaaCTCCCAGAAATTTATCGGGTTGAAATCTGGCTGAATAAAAAGAAGGTTTCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACt                                                          <  1:60321/250‑1 (MQ=255)
                                                                                                                                                                                                    aCTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAAAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTATGCAGCAACAGTTTACGGTAGAACAACTACGTAAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTAAAAAAGATGACGATAGCGTGATCAAATCGCTTGAACAAATCATTCCCGGTACAGATGCTGAATTTGATTTCAATCATCAGCACCtt                                                         >  2:180804/1‑250 (MQ=255)
                                                                                                                                                                                                        ccaggacattacctggTTGATATCTGGATAAATAAAAAGAAGTTTTCACAAAAAAAAATTCTATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAAttt                                                    <  1:541233/237‑1 (MQ=255)
                                                                                                                                                                                                         cAGGAACTTATCAGGTTGGTATCAGGCTGAATAAAAAGAAGGTTTCACCGAAAAAAATTACACTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTTAGCTGGGTATTATGGTGGATGCAATCCCGGCGCTGGCTGAAAATGATGGCGATAGCGTGATCACATCGATTAAACAACTCATTCCCGGTACAGCTGCTGATTTTGCTTTCAATCATCAGCAACTTAGTTTg                                                   >  2:534362/1‑251 (MQ=255)
                                                                                                                                                                                                           gaaatttATCAGGTTGATATCTGGCTGAAAAAAAAGAAGGTTTCACAAAAAAAAATTACTTTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAg                                                 <  1:239282/246‑1 (MQ=255)
                                                                                                                                                                                                              aCTTATCGGGTTGATATCTGGCTGAATAAAAAAAAGGTTTCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCAt                                              <  1:349985/251‑1 (MQ=255)
                                                                                                                                                                                                                     aGGTTGATATCTGGCTGAATAATAAGAAGGTTTCCCAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAAAAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGATGAATTTGATTTCAATCAGCCGCAAATTAATTTAAGCATTCCCCa                                        >  2:387724/1‑250 (MQ=255)
                                                                                                                                                                                                                         ggATAGCTGGATGAATAAAAAGAAGGTTTCAAAGAAAAAAATTATATTTCCC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAAtt                                    <  1:387724/249‑1 (MQ=255)
                                                                                                                                                                                                                               cTGGCTGAATAAAACGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGCACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGCTGACGATAGc                                                                                                                              >  2:516418/1‑154 (MQ=255)
                                                                                                                                                                                                                               cTGGCTGAATAAAAAAAAGGTTTCACAAAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGc                                                                                                                              <  1:516418/154‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTAcc                          >  1:88047/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGa                                                                                                                                     <  1:153868/141‑1 (MQ=255)
                                                                                                                                                                                                                                     gAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGa                                                                                                                                     >  2:153868/1‑141 (MQ=255)
                                                                                                                                                                                                                                       aacataagaagGTTTCACAGCAAAAAATTACATTTACC‑GCAAATGAAGAGCAACTTATGCAGCAAAAGTTTACGGTAGAAAACCTACGTGAGCTGTGTATTAAGGTGGATGAAATCCCGGAGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACACATCATTCCCGGTAAAGATGCTGCATTTGATTTCAATCATCTGCAACTTAATTTGAGCATTCCCCAACTTGCACTTTACCGTGAt                     >  2:571557/6‑251 (MQ=255)
                                                                                                                                                                                                                                                    ttCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCATCTTAATTTGAGCATTCCACAAATTGCACTGTACCGTGATGCAAGAGGTTACg        >  1:88892/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cacaGAAAAAGATTACATTTACC‑GCCTATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACCACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGAAGATAGCGTGATCAACTCGCTTGAACAAATCATTCACGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAAATTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGtc      >  2:512440/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          gAAAAACATTACATTTACC‑GCCAATGCAGGGCAACTTCTGCAGCCAAAGTTTACGGTAGAACAACTACGTGAGCTGGGCATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTAGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATCTGATTTCAATCATCATCAAATTAATTTGAGCATTCCCCAAATTGCACTGTAACGTACTGCAAGAGGATACGtacccc  >  2:181383/1‑246 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           
ACCAGGGACATAAACGCACACCGAAACAGTTGGCGCTCATCATCAAACGTTGTTTGCCGATGGTGCTCACTGGCAGCGGCATGCTTTGCACTACCGCTAACGCCGAAGAGTATTATTTCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACC‑GCCAATGCAGAGCAACTTCTGCAGCCACAGTTTACGGTAGAACAACTACGTGAGCTGGGTATTAAGGTGGATGAAATCCCGGCGCTGGCTGAAAAAGATGACGATAGCGTGATCAACTCGCTTGAACAAATCATTCCCGGTACAGCTGCTGAATTTGATTTCAATCATCAGCAACTTAATTTGAGCATTCCCCAAATTGCACTGTACCGTGATGCAAGAGGTTACGTCTCCC  >  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: ‑