Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,8230AC41.3% 54.2 / 12.7 48T6P (ACC→CCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base A (21/6);  new base C (1/18);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 7.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTA  >  NC_000913/4548577‑4549061
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ataAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTAccc                                                                                                                                                                                                                                               <  1:409333/250‑1 (MQ=255)
             gaacgcgCAACTCATCCAGTTAGCGCTACAGGATGACAGTGGCAACACATTGAATCCTGGCGCAACCAAACAAGTTCAGGCGGATGATTCCTCACAATCAGCGCACTTACCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGCGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTActg                                                                                                                                                                                                                                 >  2:256162/4‑251 (MQ=255)
                        aCATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGGCTg                                                                                                                                                                                                                     <  1:201392/251‑1 (MQ=255)
                            ccAGTTGGAGATACAGGATGAAAGGGGCAACACATTGAATACTGACGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGAACTTCACGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTAACCTGTTTGCTGGACTGCTG‑ATGGGCTGGt                                                                                                                                                                                                                   >  2:264256/1‑249 (MQ=255)
                                             aTGACAGTGGCTACACATTGAATTCTGGCGCAACCAAACCAGTTCAGGTGGATGATTCATCACCATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGACACCATTCAGGCCGTGATTAGCATCACCTATCCCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTCTTCACATGTTTGCTGTACTGCTG‑CTGGGCTGGCCGGTAAATGCCTGG‑TCAtt                                                                                                                                                                                               >  2:171853/1‑251 (MQ=255)
                                             aTGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTGCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCAGTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGGAAATGCCTGG‑TCAtt                                                                                                                                                                                               <  2:405943/251‑1 (MQ=255)
                                                       cAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTaa                                                                                                                                                                                      >  2:201910/1‑250 (MQ=255)
                                                         acacATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGACACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAAc                                                                                                                                                                                   >  1:294583/1‑251 (MQ=255)
                                                           acaTTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTAATTAGCATCACCTATACCTACAGCTGAACCC‑GCAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTCAATGCCTGG‑TCATTCGCCTGTAAAACCg                                                                                                                                                                                 >  1:248507/1‑251 (MQ=255)
                                                                       ggcgaaacaaaaccatttCAGGGGGATGATTCCCCAAAATCAGCGCACTCCCCGTTCCAGATAAGAGAATTGACAGTAAAGGGCGGGGCCCCTCAGGGACCCATTCAGGCTGTGATTAGCATCACCTATACCTCCAGCTGACCCC‑GAAGAGATGATGGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑CCATCCGCCTGTAAAACCGCCAATGGTAc                                                                                                                                                                       <  2:122958/233‑1 (MQ=255)
                                                                           aaacaaaaaCAGTTCAGGTGGATGATCCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATCCCTACAGCTGACCCC‑GAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTAt                                                                                                                                                                 <  1:178870/246‑1 (MQ=255)
                                                                                 aaaCAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTAt                                                                                                                                                           >  1:130578/1‑251 (MQ=255)
                                                                                 aaaCAATTCAGGTGGATGATTCCCCAAAAACAGCGCACTCCCCTTTTCAGGTAAGAGCATGGACTGAAAATGGCGGAGCCACTCAGGGACCCATTCAGTCAGTGATTAGCATCCCCTATACCTACAGTTGAACCC‑GAAGAGAGGGTTGTAATGAAACGATTTATTCCCCTGTTTGCTGTACTGCTG‑AGGGGCTGGCCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTAt                                                                                                                                                           <  2:86121/251‑1 (MQ=255)
                                                                                     aGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGc                                                                                                                                                       >  1:179053/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:284849/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:101434/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATAATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGTGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGCGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:318691/1‑251 (MQ=255)
                                                                                          ggggggATGATTCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAAGGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCCCCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGATTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCTTGG‑TCTTTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTgg                                                                                                                                                  <  1:92215/247‑1 (MQ=255)
                                                                                          aGGTGGATGACTCCTCACAATCAGCGAACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGAAGCCACTCAGGGAACCCTTCAGGCAGTGATTAGCATCACCTATATCTACCGCTGAACCC‑GAAGAGATGATTGTAAGGAACCGAGTAATAACCAGATTTGCTGGACTGCTG‑ATGGTCGGGTCGGTAAATGCTTGG‑TCATTCGCCTGTAACACCACCAATGACACCGCTATCCCTATTGACGGTgg                                                                                                                                                  >  2:224332/1‑251 (MQ=255)
                                                                                            gTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGc                                                                                                                                                 >  1:250125/1‑250 (MQ=255)
                                                                                             tGGATGATTCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTCCAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTggtcc                                                                                                                                               >  1:42293/1‑248 (MQ=255)
                                                                                                atgatTCATCACAATCAGCGCCCTTCCCGTTACAGGTCAGCGCATTGACAGTAAATGGCGGAGCCACTCAGTGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAATGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCACTGGTACCGCTATCCCTATTGGCGGTGGccgacc                                                                                                                                            >  1:322444/1‑246 (MQ=255)
                                                                                                  gatTCCTCACAATCAGCGCCCTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGAAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTACATGCCTGG‑TCATTCGCCTGTAAAACCGCCACTGGTACCGCTATCCCTCTTGGCGGttggcgctcc                                                                                                                                           >  2:95263/1‑241 (MQ=255)
                                                                                                  gatTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCGGCGccc                                                                                                                                          >  1:77574/1‑250 (MQ=255)
                                                                                                                  cgcATTTCCCGTTACAGGTAAGAGCATTGAAAGAAAATGGAGGAGCCACTCAGGGAACCATCCAGGAATTGATTAGAATCACTTATACCTACAGCTGACCCC‑GAAGAGATGATGGTAATGAACCGAGTTATTACCCTGTTTGCTGTACTTCTG‑ATGGGCTGGTCGGTAATGGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCt                                                                                                                          <  1:284059/251‑1 (MQ=255)
                                                                                                                    cACTTCCCGTTACAGGTCAGAGCATTGAAAGTAAATGGCGGAGCCACTCAGGGACCCTTTCAGGCAGTGATTAGCATCCCCTATACCTACAGCTGACCCA‑GAGGAGATGATTGTAATAAAACGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCTTTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTg                                                                                                                        <  2:109714/251‑1 (MQ=255)
                                                                                                                              tacagtccAGAGCATTGACAGTAAAGGGCGGAGCCACCAAGGGACCCATTCAGGCAGTGATTAGCATCACTTATACCTACAGCTGACCCC‑GAAGGGATGATTGTAATAAAACGAGTTATTCCCCTGTTGGCTGTACTGTTG‑ATGGGCTGGTCGGTAAATGCTTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCcgtcg                                                                                                               <  1:117343/243‑1 (MQ=255)
                                                                                                                                  ggTCAGAGCATTGCCAGTAACTGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTCCAGCTGAACAC‑GACGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCACTGTTTAAGTAAACCCTGCGCCCCTCGTgttt                                                                                                          >  2:73562/1‑248 (MQ=255)
                                                                                                                                     cAGATCATGGACATTAAAGGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATCCCTCCAGCTGAACCC‑GAAGTGATGATTGTAATGAACCGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGTTAATTGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCTCCCGTCGTGAAtgt                                                                                                        <  1:118933/250‑1 (MQ=255)
                                                                                                                                           aTTGACAGTAAAGGGAGGCGCCTCCCAGGAAACAATGCAGGCGGTGATTACCATCCCCTATCCCTCCAG‑TGGACCCTAAATAGATGATTGTAAGGAAACGATTAATCTCCCTGTTCGCGTTACGGCTA‑AGGGTCTGGTCGGAAATTGCCTGG‑TCTTTCGCCAGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAGCCTTGCGCCCGTCGTGAATGTGGGGCaa                                                                                                 <  1:171853/251‑1 (MQ=255)
                                                                                                                                              gACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAATGGTTTTGTAAACCTTTCTCCCGCCGGGGATGTggggggaaa                                                                                               >  2:346328/1‑245 (MQ=255)
                                                                                                                                                   aaaaTGGCGGAGCCACTCAGGGAACCATTCAGGCAGTAATAAGCATCACCTATCCCTACAGCGGAACCC‑GAAGAGATGATTGTAAGGAAACGAGTTATTCCCCTGTTTGTTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCTTGG‑TCATTCCCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGGGAATGTGGGGCAAAACCTGGt                                                                                         <  1:324419/250‑1 (MQ=255)
                                                                                                                                                    aaaTGGAGGAGCCACTCAGGGAACCATGCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGGTTGGAATGAAACGAGTTATTCCCTTGTTTGCTGTCCTGCGG‑ATGGGCTGGTCGTTAAATGCGGGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTc                                                                                        <  1:61098/251‑1 (MQ=255)
                                                                                                                                                                cccTCAGGGAACCATCCAGGAAGTGATTAGCATCCCCTATCCCTACAGCTTAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTGGCTTTCCTGCTG‑ATGGGCTGTTCGGTAAATGCCTGT‑TCATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATTTGGGGCAAAACCTGGTCGTGGATCTTTCg                                                                            <  1:346328/249‑1 (MQ=255)
                                                                                                                                                                    cagggatctaatCAGGCAGAGATGAGCATCGCCTATCCCTACATCTAACCCC‑GAAGAGAGGATGGTAATGAAACGAGTTATTCCCCTGTTTGGTGTACTGTGG‑AGGGTGGGGTCGGAAAATGC‑TGGCTCATTCGCCGGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCGTGCTCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGAc                                                                          <  1:257740/238‑1 (MQ=255)
                                                                                                                                                                                   gcagtgttagacctCACCTATTCCTACAGCTAAACCC‑GAAGAGATTATTGAAATGAAAGGGTTTATTCCCCTGTTTGCTGTACTGCTG‑AGGGGCGGGTCGGTAAAGGCCTGT‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATa                                                         <  2:130578/238‑1 (MQ=255)
                                                                                                                                                                                             gaatccccTATACCTACAGCTGAACCC‑AAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTTTACTGCGT‑AGGGGCGGGTCGGTAAATCCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATcc                                               <  1:131444/245‑1 (MQ=255)
                                                                                                                                                                                                  cctgttaatatcatgtGAACCCTGAAGGGATGATGTTAATGCACCGATTTATACCCCTGTTTTCTGTACTTCGG‑AGGGCCTGGTCTGAAAACTCCGGG‑TCTTCCGCCTTTAAAACCGCCAATGGTACCGCTATCCCTACTGGCGGTGGCAGCGCCAATGTTTATTTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACTCAAATCTGTTGCCATAACGATTAGCCGGaa                                           <  2:294583/236‑1 (MQ=255)
                                                                                                                                                                                                      atcccTACAGCTGACCCC‑GAAGAGAGGATTGAAATGAAACGTGTTATTCCCCTGTTTGCTGTTC‑GCCGTAGGGGCTGTTCGGTAAAGCCCTGG‑CCATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGGGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAGCGTTTGCCATAACGATTATCCTGAAAcc                                        <  2:161573/246‑1 (MQ=255)
                                                                                                                                                                                                      ataCCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCCAGTTTATTTAAAACTTGCGCCCGTCGTGGATGTGGGGCAAAAACTGGTCGTGGCTCTTTTTCCGCAAATATTTTTCCCGAACAATTATCCCGTAACCAt                                      >  2:248872/1‑251 (MQ=255)
                                                                                                                                                                                                        cccTACAGCTGACCCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCGGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTa                                    <  1:393754/250‑1 (MQ=255)
                                                                                                                                                                                                         ccTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGTTAACCTTTCGCCCGTCGTGAATGTGGGGCCAAACCTGGGCGTGGGTCTTTCCACGCAAATCTTTTTCCATAACGATTATCCGGAAACCATTAc                                   >  1:55933/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCCC‑GAAGAGATGATTTTAATAAAACGATTTTTTTCCCTGTTTGCGTTCCTGTTA‑AGGGGCTGTTCGGAAAATTCTTGG‑TCTTCCCCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAGCCATTACAGACTATGTc                         <  1:351177/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaCCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAATGTTTTTGTAAAACTTTCCCCCCTCGGGAATGTTGGGCAAAACCTGGTCGTGGATCTTTTCACGCCAAACTTTTTCCCTAACGATTTTCCCGAAAcccataaaaactatttc                         >  1:282577/1‑236 (MQ=255)
                                                                                                                                                                                                                     ccccgaatagagtatTGAAATAAAACGAGTTATTACCGTGTTTGCTGTACGGCGG‑ATGGGCTGTTGGGTAAAGTCCTGG‑TCATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTca                        <  2:211418/238‑1 (MQ=255)
                                                                                                                                                                                                                           agagATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGTTG‑AGGGGCTGGTCGGTAAAGGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGc                  <  1:201910/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gAAACGAGTTATTCCCTTGTTTGCTTTACTGCTG‑AGGGGCTGGGCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCt    <  1:131409/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCGAGTTATTCCCCTGTTTGCTGTACTGCTG‑ATGGTCTGGGCGGTAAATCCCTGG‑TTATCCCCTTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa  <  2:146927/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCC‑GAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTG‑ATGGGCTGGTCGGTAAATGCCTGG‑TCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTA  >  NC_000913/4548577‑4549061

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