Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,9320AT23.1% 120.1 / 10.5 52Y42F (TAT→TTT) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (8/32);  new base T (12/0);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 6.10e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAA  >  NC_000913/4548706‑4549165
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
aGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATG‑TTTTTGTAAAACCTTg                                                                                                                                                                                                                                  >  1:1706757/1‑239 (MQ=255)
  gTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGGCGCCCCCAAG‑TTTTTTT‑aaaaccttcg                                                                                                                                                                                                                                >  1:1485058/1‑231 (MQ=255)
  gTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCAGCGCCAAAG‑TTTTTGT‑AAAACTTgcg                                                                                                                                                                                                                                >  1:1701342/1‑238 (MQ=255)
         cATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAGCGCCAAAG‑TTTATGT‑AAACCTTGCG‑gcccccgc                                                                                                                                                                                                                       >  1:114043/1‑234 (MQ=255)
                  tAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCCTC‑GTTATTT‑AAACCTTTCT‑GCCGTCGTCGAGGTggg                                                                                                                                                                                                              >  2:1091305/1‑239 (MQ=255)
                                                aGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑GTTATTT‑AAAACTTGCG‑CCCCTCGCGAAAGTGGGGCAAAAACTGGTCGGGGGTCtttttgcgcc                                                                                                                                                                                >  1:366646/1‑232 (MQ=255)
                                                  gCAGTGATTAGCATCACCTATACCTACAGCCGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAGCGCCCAAG‑TTTTTGT‑AAAACCTGCG‑CCCCTCCGGGATGTGGGGGCAAACCTGGGGGGGGGTCTTTTC‑ACGCaaa                                                                                                                                                                             >  1:156198/1‑239 (MQ=255)
                                                       attgcgCACCACCTATACCAAGAGCTACACACGAAGGTGTGTTTGAAAGAAAACAGTTTATTCCCTTT‑TTTTCTGTATTGTTG‑TGGGGTGTGTGGGTAAA‑TCCCTGGTCATTCGCCAGA‑AAAT‑CCGCCAATGGTACCGCTGTCGCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑ACTGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCttt                                                                                                                                                                        <  2:1701342/234‑1 (MQ=255)
                                                                       aCCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGGAGCGCCCAAG‑TTTTTTT‑AAACCTTTCG‑CCCCCCGCGAATGTGGGGCAAAAACTgg                                                                                                                                                                                                   >  1:574433/1‑197 (MQ=255)
                                                                        ccaacgacagAACCCGAAAAGAGGATTGTAGTGAAACGTGTTTTCACCCGT‑TTTGCTGTACTCCTGA‑GGGGTG‑GTTGGTAAATTGCTTGGTCTTGCGCC‑GT‑AAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCt                                                                                                                                                                                                     <  2:605007/185‑1 (MQ=255)
                                                                                  aaCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATG‑TTTATTT‑AAACATTTCG‑CCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGag                                                                                                                                                                                           >  1:605007/1‑193 (MQ=255)
                                                                                   aCCCGGAGGTATGATTGTAATGAAGCGATTTA‑TACCCTGTTTTGCGGGCCGTCTGAGGGGCTG‑GTGGTTAAG‑TGCCTGGCCTTTCGCATAA‑AAAA‑CCGCCAATGGTATCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTa                                                                                                                                            <  2:1644302/239‑1 (MQ=255)
                                                                                              tgtttGTAATGAAACGAGTTATTACCCTT‑TTTGCTTTTCTGGTGGTGGGCGG‑GTGGGTAAA‑TGCCTGGTCATTCGCTTGA‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGt                                                                                                                                  <  1:1273342/235‑1 (MQ=255)
                                                                                                    tAATGAAACGATTTTTTCCCTTT‑TTTGTTTTCCTGCTGAGGGGCTG‑GTCGGTAAA‑TGCCTGGTCTTCCGCCTTT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTcaca                                                                                                                              <  1:104324/236‑1 (MQ=255)
                                                                                                      aTGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCCGCGCCCCTG‑TTTTTTT‑AAAcc                                                                                                                                                                                                                                     >  1:1817097/1‑133 (MQ=255)
                                                                                                                 ttttccccTT‑TTGGCTTTACTGCTGAGGGGTGG‑TTCGGAAAA‑TGCGTGCTCTTTCCCCTAA‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCt                                                                                                                                                                                                                                    <  2:138163/118‑1 (MQ=255)
                                                                                                                            tttGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGGGCGCCCCAG‑TTTTTGT‑AAACCCTTCT‑CCcgtc                                                                                                                                                                                                                         >  1:138163/1‑123 (MQ=255)
                                                                                                                                   tACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCTACCCCAACG‑TTTATGT‑AAACCTTTTG‑CCcctcgc                                                                                                                                                                                                                       >  1:798540/1‑117 (MQ=255)
                                                                                                                                    cctgtgtgggggCTG‑GTCGGTAAA‑TCCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTAt                                                                                            <  1:1600827/230‑1 (MQ=255)
                                                                                                                                        ctgATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCAGCGCCCATG‑TTTTTTT‑aaa                                                                                                                                                                                                                                       >  1:946873/1‑98 (MQ=255)
                                                                                                                                        ctgATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCCCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCTGCGCCCAAT‑TTTTTGT‑AAACCCTCCG‑CGcgtcgggt                                                                                                                                                                                                                     >  1:776482/1‑112 (MQ=255)
                                                                                                                                        cggaggggCTG‑CTCGTTAAA‑TGCCTGGTCATTCGCTTGA‑AAAC‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGa                                                                                                                                                             <  1:1897135/165‑1 (MQ=255)
                                                                                                                                                  g‑gTGGGTAAA‑TGACGGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                               <  1:2089219/238‑1 (MQ=255)
                                                                                                                                                  g‑gTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTTGCGGTGGCAGCACCAATG‑TTTTTTT‑ATACCTTTTG‑CCCCCCTTGTATGTGGGGGAAAACCCGGTTGGGGGTCTTTTG‑TCGCACAACCTTTTCCCTAAAGAATATCCCGGAAACATTACAGACTCATTCGCACTGCAACAAAGCGCGGCGTCTTGCGGCGAGTTGTCTTAttttttttgg                                                                              >  2:2089219/1‑234 (MQ=255)
                                                                                                                                                  g‑gTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGGTGGAGCCCCAACGTTTTTTGA‑AAACCTT‑CG‑CCCGGCGTGGATGTGGGGGGAAACCTGGTCCTGGAGCTTTTT‑ACGCAAAACTTTTTCCCTAACGCATATACggga                                                                                                                                                   >  2:1897135/1‑168 (MQ=255)
                                                                                                                                                    tacGGTAAA‑TGCGAGGTATTACACCTAA‑AAAA‑CCTCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATaa                                                                                                                                                                <  1:2060575/154‑1 (MQ=255)
                                                                                                                                                     tgggTAAT‑TCCCTGGTCATCCCCCTGA‑AAAC‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:1007938/236‑1 (MQ=255)
                                                                                                                                                     tgggTAAA‑TGCCTGGTCTTTCCCCTAT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:1091305/236‑1 (MQ=255)
                                                                                                                                                     tCGGTAAT‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:1900161/238‑1 (MQ=255)
                                                                                                                                                     tCGGAAAT‑TGCTTGTTAATTCGCTTGT‑AAAC‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:1121818/238‑1 (MQ=255)
                                                                                                                                                     gtggTAAG‑TGCCTGTCCATCCCCTTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:1487/236‑1 (MQ=255)
                                                                                                                                                            aatggcctgTCCTTCCCCCGT‑AAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTGGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaa                                                                       <  1:656195/231‑1 (MQ=255)
                                                                                                                                                            aaatccTGGTTATTCCCCCTTAAAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                        <  1:1512636/234‑1 (MQ=255)
                                                                                                                                                               gCCTGGTCTTTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCa                                                                                                                                                                                                           <  1:1409220/104‑1 (MQ=255)
                                                                                                                                                                    gTCATTCCCCTGT‑AAAA‑CCGCCAATGGTTCCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGa                                                                                                                                                    <  1:301453/153‑1 (MQ=255)
                                                                                                                                                                     tCATTCCCCTTA‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:147160/223‑1 (MQ=255)
                                                                                                                                                                      tgttCGCCTTT‑AAAC‑CCGCCAATGGTCCCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                            <  1:78951/237‑1 (MQ=255)
                                                                                                                                                                         tCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCCACG‑TTTATTT‑AAAACCTTCC‑CCCCTCGTGGATGTGGGGGAAAAACCGGTc                                                                                                                                                                                                 >  2:1409220/1‑104 (MQ=255)
                                                                                                                                                                            cccttaaaaa‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                       <  1:440649/233‑1 (MQ=255)
                                                                                                                                                                              tGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCGGCGGCAAAG‑TTTTTGT‑TAACCCTGCGCCCCCCCGGGAAGGGGGGGCAAAACC‑GGGCGGGGGGCTTTTG‑GCGCAAAACTTTTTTCCTAAAGAGTATCCGGGAAACATTata                                                                                                                                          >  2:301453/1‑151 (MQ=255)
                                                                                                                                                                                t‑AAAA‑CCGCCAATTGTACCGCTATCCCTATTGGCGGTGGCAGCTCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                  <  1:523182/239‑1 (MQ=255)
                                                                                                                                                                                  aaac‑CCGCCAATGGTCCCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                  <  1:322323/238‑1 (MQ=255)
                                                                                                                                                                                   aaa‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                <  1:74411/239‑1 (MQ=255)
                                                                                                                                                                                               ggTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACa                                                                                                                                          <  1:694595/138‑1 (MQ=255)
                                                                                                                                                                                                         tCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                           <  1:582246/239‑1 (MQ=255)
                                                                                                                                                                                                         tCCCTATTGGCGGTGGCAGCGCCAAAG‑TTTATGT‑TAAACCTGGG‑CCCGTCGGGAAGGTGGGGCCAAACCCGGTCGGGGGGTTTTCGCAGCCAACTTTTTTCCCCTAA‑GAATATCCCGGGACCCTTACAGACCATGTcc                                                                                                                                >  2:694595/1‑137 (MQ=255)
                                                                                                                                                                                                             tATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:959243/185‑1 (MQ=255)
                                                                                                                                                                                                             tATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                            <  1:1270839/202‑1 (MQ=255)
                                                                                                                                                                                                                        gCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat            <  1:683211/239‑1 (MQ=255)
                                                                                                                                                                                                                            cgcCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGAGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGTCTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGCGTTTTCTAATTTTTCCCGGACCGTAACATACAATTGCCGGAGCTTACCATTTCCTAACCCCAGCGAAAACCCACTCGGTGTTTATaac         >  2:683211/1‑237 (MQ=255)
                                                                                                                                                                                                                             gcCAATG‑TTTATGT‑AAACCGTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat            <  1:794983/234‑1 (MQ=255)
                                                                                                                                                                                                                                   tggttATGT‑AACCCTTCGG‑ACCGTCGTGAATGTGGGGCAAAACCTGGTCTGGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa  <  2:366646/236‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
AGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTG‑TTTGCTGTACTGCTGATGGGCTG‑GTCGGTAAA‑TGCCTGGTCATTCGCCTGT‑AAAA‑CCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATG‑TTTATGT‑AAACCTTGCG‑CCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCG‑ACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAA  >  NC_000913/4548706‑4549165

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