Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,9320AT28.6% 57.6 / 10.7 35L412F (TTA→TTT) sfmDputative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand):  ref base A (8/17);  new base T (10/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 3.44e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTATCGTCGAATGGAGGGGTACGAATATGATTACGACGGTGAGC  >  NC_000913/560714‑561182
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
gATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGCAGGGCGTTTCAAATTTGCCCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGGTTACCACAGGAATTTACTGCCTACGGCGGGGAGCAATTATTTTCC‑AATTTCACCCCCCTTTTTTTT‑TGGGCTGGGGGGCacttccgggactgggg                                                                                                                                                                                                                              >  2:430148/1‑236 (MQ=255)
    gCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTCATAGTCAGCATTCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTCCCACAGGAATTTACTGCCTACGGCGGGACGGAATTTTCTGCC‑TATTACACCGCC‑TTTTTTTT‑AGGGCTTGGGGGCCAACTCggggacctgggcgcg                                                                                                                                                                                                                         >  1:203061/1‑239 (MQ=255)
             agAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCCATTATTTGCC‑AATTTCCCCCCC‑CTTTTTTT‑TGGGGTGGGGGGCAATCTCCGGGACCGGGGCGCCgtgttgccg                                                                                                                                                                                                                >  1:254082/1‑246 (MQ=255)
               aaTTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGa                                                                                                                                                                                                              <  2:254082/251‑1 (MQ=255)
                           ccGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTTTCTGCC‑AATTTCACCGCC‑CTTTTTTT‑TGGGCTGGGGCGCAATCTCGGGGACTGGGGCGCCGTGTCGCTGGGTGTAACGCA‑Tgc                                                                                                                                                                                                 >  1:191601/1‑251 (MQ=255)
                                   aaCAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGcagtc                                                                                                                                                                                         <  1:98378/251‑1 (MQ=255)
                                       gTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGAAGCAATTATCTGCC‑CAATAAACCGCC‑TTTTTTTT‑AGGGGTGGGGGGCAATCTCCGGAACTGGGGCGCAGGGTCGCTGGATGTAACGCA‑Tgcgcggcgtcaggt                                                                                                                                                                                     >  1:113354/1‑242 (MQ=255)
                                               tttACAACGCGAAGGGCGTTTCACATTTGACCTGACGGCGGGCGATTT‑CGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTTTCTGCC‑AATTACACCGCC‑TTTTTATT‑TGGGCTGGGGGGCCATCTCGGGGACTGGGGCGCAGTGTCGCTGGGTGTAAAGCC‑TTCGCGCCGTCCGTTAGCCCACGc                                                                                                                                                                            >  1:118814/1‑251 (MQ=255)
                                               tttACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGAc                                                                                                                                                                              <  2:191601/250‑1 (MQ=255)
                                               tttACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTGTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCTGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGAc                                                                                                                                                                              <  1:24485/250‑1 (MQ=255)
                                                    aaCGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCCCTCGGCGGTTTCCCACAGGCATTTCCTGCCTACGGCGGGACGCAATTATCCGCCCAATAAC‑CCGCC‑TTTTTTTT‑AGGGCTGGGGCGCAATCTCGGGGACTCGGGCGCAGGGTCGCTGGGTGTAACGCC‑TGCCCGCCGTCAGGTAGCCCACGCCacg                                                                                                                                                                        >  1:413322/1‑249 (MQ=255)
                                                       gcgAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAgtcgt                                                                                                                                                                     <  2:113354/251‑1 (MQ=255)
                                                         gAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGGCGCAATTATCTGCC‑CATTACACCGCC‑TTTTTTTTTAGGG‑TGGGGGGCAATCTCGGGAAATGGGGCGCAGTGTCGCTTGGTGTAACGCA‑TGCGCGCCGTCAGTTTGCCGGCGCCAgtcggc                                                                                                                                                                    >  1:229578/1‑248 (MQ=255)
                                                                                     cgggcgATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATTTGCC‑AATTTCACCGCC‑TTTTTATT‑AGGGCTGGGGGGCAATATCGGGAACTTGGGCGGAGTGTCGCTGGATGTAAAGCC‑TTCGGGCAGGCAGGTAGCCGACGGCCGTCGTCATGAGGGGGGTTCTATTCGCttacgctat                                                                                                                                       >  1:61645/1‑244 (MQ=255)
                                                                                              ttCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCg                                                                                                                              >  1:316332/1‑251 (MQ=255)
                                                                                               tCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCg                                                                                                                              <  1:39860/250‑1 (MQ=255)
                                                                                                   aGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTTTT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGGCAGTTAGcc                                                                                                                                                                                 >  1:388324/1‑195 (MQ=255)
                                                                                                   aGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGcc                                                                                                                                                                                 <  2:388324/195‑1 (MQ=255)
                                                                                                      ggTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGTGTACGGCAGTCGGCGGTTTACCACAGGAATTTTCTGCCTACGGCGGGACGCAATTATCGGCC‑AATTACACCGCC‑TTTTTAGT‑AGGGCTGGGGCGCAATCTCGTTAACTGGGGCGCAGTGGCGCTGGAGGTAACGCA‑TGCTCGCAGTCAGTTATCCGACGCCAGTCGTCATGAGGTGGATTCTATTCGCGTCCTCTATGCGAAATCGATGAacac                                                                                                                      <  1:415856/251‑1 (MQ=255)
                                                                                                         aaTAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCCCTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTAGCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCtt                                                                                                                   <  2:2923/251‑1 (MQ=255)
                                                                                                            aGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCTCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAGTCGATGAACACCTTCg                                                                                                                 <  1:395982/250‑1 (MQ=255)
                                                                                                                                  tttttCAGGGGACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAAGCGAGGAACACCTTCGGCACCAATGTTCAGTTAATg                                                                                             <  2:197908/248‑1 (MQ=255)
                                                                                                                                    tttCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACCGCCTACGGCCGGACGCAATTCTCTTCC‑AATTTCACCGCC‑TTTTAATT‑CCGGATGGGGCGCAATCTCGGGAACTGGGGCGCCGTGCCGCTGGGCGTTACGCCCTCGGGGCGTTCA‑TTAACCGACCCCAGTCGGCATGAGGGGGGTTCTATTGGCTTCCTCTATTGAAAATCGATGGACACCTTCGGCCCCAAATTTTAGTTTATGGGttt                                                                                        >  1:188512/1‑250 (MQ=255)
                                                                                                                                       cAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGAGGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTAcc                                                                                      <  1:430148/250‑1 (MQ=255)
                                                                                                                                           gTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCCCC‑TTTTTATT‑AGGGCTGGGGCGCAAACTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCCGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTAt                                                                                 >  1:105545/1‑251 (MQ=255)
                                                                                                                                               ggCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGCCGCAATTATCTGCC‑AATTACACCCCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGCCGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCg                                                                              >  1:387820/1‑250 (MQ=255)
                                                                                                                                               ggCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGAGGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGa                                                                             <  1:411826/251‑1 (MQ=255)
                                                                                                                                                  aCTCGGCGGTTCACCACAGGAATTTACTGCCTCCGTCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTCTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCGTCCTCTTTGTGAAATCGATGAACACCCTCGTCACCAACTTTCATGGAATGGTTTACCGCGATTCGacac                                                                          <  2:203061/251‑1 (MQ=255)
                                                                                                                                                      ggcggTTTACCACAGGAATTTACTGCCTACGGCGGGACGCCATTATCTGCC‑CATTACACCGCC‑TTTTTTTT‑AGGGCTGGGGCGCCATCTCGGGGACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCCCCAATTTTCCGTTAATGGGTTACCGCTATTCGTCACAAgg                                                                      >  2:152675/1‑251 (MQ=255)
                                                                                                                                                       gcggTTTACCACAGGAATTTACTGCCTACGGCGGGACAAAATAAAATTCC‑CATTACACCGCC‑CTTTTTTT‑CGGGCTGGGGCGCAATCTCGGGAACTGGCGCGCAGAGTCGCTGGATGTAACGCA‑CTCGCGCAGTCAGTTAACCGACGCCAGTCGTCCTGAGGGGGATTCTATTCGCTTCCTCTATGCGGTATCGAAGAACACCTTCGGAACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACaacgg                                                                     >  1:430463/1‑248 (MQ=255)
                                                                                                                                                                                   cggcggGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTAtc                                         <  2:357294/251‑1 (MQ=255)
                                                                                                                                                                                       gggACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTTATGATGTTGCGTAtcgtcg                                     <  1:126814/251‑1 (MQ=255)
                                                                                                                                                                                        ggACGCAATTATCTGCA‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTAtcgtcg                                     >  1:377912/1‑250 (MQ=255)
                                                                                                                                                                                                    cTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGATATATCGGCAACTGAAGCGCAGTGCCGCTGGACGTTCCGGA‑TGCACGCAGCCGGCTAGCTCACCCCCGTCGGCATGCGGGGGAATCTATTGGATTCATCTATGCGAACTCGATGAACACCTTCAGCACCACTTTTCAGCTGATCGGTTACAGATAAACGCCACAAGCTTTTTATACCCTTGTTGATGTTGCGTATCGTCGAAAGGAggggacc                        >  1:205315/1‑248 (MQ=255)
                                                                                                                                                                                                             tACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTATCGTCGAATGGAGGGGTACGAATATg                 <  1:239455/250‑1 (MQ=255)
                                                                                                                                                                                                                           tatt‑aGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGc                                              <  1:354450/208‑1 (MQ=255)
                                                                                                                                                                                                                           tatt‑aGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGc                                              >  2:354450/1‑208 (MQ=255)
                                                                                                                                                                                                                            att‑aGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTCTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGACAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTATCGTCGAATGGAGGGGTACGAATATGATTACGACGGTGAGc  >  1:60872/1‑251 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
GATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGGGCGATTTTCGTAGCGGTAATAGTCAGCAATCATCGCCTTTCTTTTTTCAGGGTACGGCACTCGGCGGTTTACCACAGGAATTTACTGCCTACGGCGGGACGCAATTATCTGCC‑AATTACACCGCC‑TTTTTATT‑AGGGCTGGGGCGCAATCTCGGGAACTGGGGCGCAGTGTCGCTGGATGTAACGCA‑TGCGCGCAGTCAGTTAGCCGACGCCAGTCGTCATGAGGGGGATTCTATTCGCTTCCTCTATGCGAAATCGATGAACACCTTCGGCACCAATTTTCAGTTAATGGGTTACCGCTATTCGACACAAGGTTTTTATACCCTTGATGATGTTGCGTATCGTCGAATGGAGGGGTACGAATATGATTACGACGGTGAGC  >  NC_000913/560714‑561182

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