Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,7510AG28.6% 92.9 / 11.0 49intergenic (‑89/‑92)dcuS/yjdIsensor histidine kinase DcuS/PF06902 family protein YjdI
Reads supporting (aligned to +/- strand):  ref base A (4/31);  new base G (14/0);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 4.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGCCGTGG  >  NC_000913/4351517‑4351986
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
cgcgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAAc                                                                                                                                                                                                                                   <  1:148107/250‑1 (MQ=255)
          cACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCa                                                                                                                                                                                                                       <  2:571246/251‑1 (MQ=255)
              gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAAATTCAAG‑GGGGGAAAGA‑TGGTTTGCTTTTT‑AGGGGGAGAGAACT‑CTGGCAAGGG‑aaatggcagggc                                                                                                                                                                                                                    >  1:78044/1‑242 (MQ=255)
              gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGc                                                                                                                                                                                                                   <  2:132917/251‑1 (MQ=255)
                        aaGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGGAGAGA‑AT‑CTGGCAACGG‑AAAATGCCAGTGCTGTATGa                                                                                                                                                                                                            >  1:524868/1‑248 (MQ=255)
                          gTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGTTACTCAAGATCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCGCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑Att                                                                                                                                                                                                      <  2:548320/251‑1 (MQ=255)
                                          accaATAGCACCGAGAACAGTACCGCACTGACAATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTtat                                                                                                                                                                                    <  2:365403/251‑1 (MQ=255)
                                             aaTAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGc                                                                                                                                                                                 <  1:23217/251‑1 (MQ=255)
                                                 gCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑TTCAAAAGGGGGTA‑GAATGGTTTGCGATAT‑ATGGGGAGGGA‑AC‑ATGGCAACGGGAAACT‑CCAGTGCCGTATGAAG‑ATTCCCGGGG‑C‑ACGCTTATAGCGAtt                                                                                                                                                                             >  2:96656/1‑250 (MQ=255)
                                                   aCCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAAt                                                                                                                                                                           <  2:81548/251‑1 (MQ=255)
                                                   aCCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAAGTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCATCGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAAt                                                                                                                                                                           <  2:188086/251‑1 (MQ=255)
                                                                 cgctcTGACCATTAAGATATCTGTTGTACTCAATGTCCACGGACGGTTTCGTAACATGCGGTAGGTCAATTAATGTCTCATCAGCGTCCTGGGGTGACAAATTTCTTAAGCGTTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑AGTCC‑GGGG‑CTATGCTTATAGCTATAATCATACTGATgagag                                                                                                                                                             <  2:180196/247‑1 (MQ=255)
                                                                 cGCACTGACCATTAAGATCACTGTGGTACTCAATTTCAGCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTAGGCTTATAGCGATAATCATACTGATgagag                                                                                                                                                             <  2:513859/251‑1 (MQ=255)
                                                                       gACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGAGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑ATCAAA‑GGGGGGAAGAATGATTTGGTATAT‑AGGGGGAGAAG‑CC‑CTGGCAAAGG‑GAAATGCCAGAGCTGTAAGCAG‑AATCC‑GGGGGTTATGCT‑ATAGGGAGAATTATACTGGTGGGAGAggg                                                                                                                                                         >  1:37290/1‑249 (MQ=255)
                                                                        aCCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑AGGAGGAAGAATGATTGGCTATAT‑AGGGGGAGAGA‑CT‑CTGGGAACGG‑AAACATCCAGGGCTGGATGAAG‑GTTCC‑GGGG‑GTACGCTTATAGCGGTAATTCTAATGTTGAGAGAgggg                                                                                                                                                        >  2:49304/1‑248 (MQ=255)
                                                                           aTTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTc                                                                                                                                                   <  1:187012/251‑1 (MQ=255)
                                                                           aTTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTc                                                                                                                                                   <  1:578296/251‑1 (MQ=255)
                                                                                      tgtgGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAGTTCAA‑‑GGGAGTAAGAATGATTTGGTTTAT‑AGGGGGAGGGA‑CC‑CTGGCAACGGGAAATTACCAG‑GCTGTATGGAG‑GTTCC‑AGGG‑CTATGGTTATAGCGATAAATATACTGGTGGGAGAGGGAAGGGGATGGATccggg                                                                                                                                        >  2:438467/1‑247 (MQ=255)
                                                                                       ttgggaCTCAATGTCATCGGACGTTTGCGTAACATGCGGTAGTGCAATGAATGTCTCAGCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGGAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGgc                                                                                                                                        <  2:139748/245‑1 (MQ=255)
                                                                                           tACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCCTCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTTTAA‑GGGAGTAAGAATGATTGTCTTTAT‑AGGGGGAGGGA‑CTACTG‑CACCGG‑AAAATGCCCGGGCTGTATGGAG‑AGTCC‑GCGG‑CTAAGCTTTTAACGATAATCATTCTTTTTAGAGAGGGGGGGTCATGGGTCTGGCGCtc                                                                                                                                    >  2:78265/1‑249 (MQ=255)
                                                                                                aaTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTTATGAGGCGGGTAA‑TTCAAA‑GTGAGTAAGAATGTTTGGCTATAT‑AGGGGAAGAGT‑CT‑CTGTCAACGG‑AAACTTCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTAGAGCGATACGCATGCTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTAATGGAc                                                                                                                              <  2:48825/251‑1 (MQ=255)
                                                                                                 aTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTTAAA‑GGGGGTAAGAAAGATTTGCTTTTT‑TGGGGGAGAGA‑CT‑CTGGCAACGG‑AAACTGCCCGTGCTGTATTAAG‑ATTTC‑GGGG‑GTTTGCTTTTATCGATAATCATACTTTTGAGAGAGGGAGGGTCATGGTTCAGGCGCTACt                                                                                                                                  >  2:157377/1‑246 (MQ=255)
                                                                                                 aTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACt                                                                                                                                  <  1:157377/246‑1 (MQ=255)
                                                                                                         ggACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑ATTAAA‑GGGGGGAAGAATGATTGGGTTTAT‑AGGGGGAGAGA‑CT‑CCGGCAACGG‑GAAATGCCCGTGCTGTATTAAG‑ATTTC‑GGGG‑CTATGCTTTTAGAGGTAATCATACTGATGAGAGAGGGAAGGGCATTGGTCAGGGGCTACTGGACGGGGG‑TTAt                                                                                                                    >  1:287906/1‑251 (MQ=255)
                                                                                                          gACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTTGGACAAATTTCTTAAGCATTATCGCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAGCTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGTG‑TTATc                                                                                                                   <  2:419943/251‑1 (MQ=255)
                                                                                                           aCGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTTAAA‑AGGGGGAAGAATGATTGGCTTTTT‑AGGGGGAGAGA‑AT‑TTGGCAACGG‑AAACTGCCCGGGCTGGTTGAAG‑ATTTC‑CGGG‑CCATGCTTATAACGGGAATTATTATGATGAGAGAGGGAAGGGCATGGGTAAGGGGCTCCTCGAGGGGGG‑GTATCg                                                                                                                  >  2:159706/1‑251 (MQ=255)
                                                                                                                 gagTAACAGGCGGTAGGGCAATGAAGGACTCATCAGCTTCCTTTTGTGACAAATTGCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGTCGATAT‑AGGGGAGGAGA‑CT‑CTGGCAACGG‑AAGCTGCCAGTGCTGTATGAAG‑AGTCT‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGATAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTta                                                                                                            <  2:310539/249‑1 (MQ=255)
                                                                                                                        aTGCGGTCGGGCAATGACTGTCTCTACATCTTCCTTGTGTGACAAATTGCTTAAGCATTATCTCAGATGAGGCGGGTAA‑TTTAAA‑GGGATTAAGAATGATTGGTTATAT‑AGGGGAAGAGA‑CT‑CTGGCATCGG‑AAACTGCCAGTGCTGTATGAAG‑AGTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGATAGGTAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATAcctgc                                                                                                     <  2:127488/251‑4 (MQ=255)
                                                                                                                                       gaatgTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCCAC‑GGGAGTAAGAATGATTGGGTATATCAGGGGA‑GAGA‑CT‑CTTGTAACGG‑AAAATGCCCGAGGTGTATGAAGGATTC‑‑GGGG‑CTATGCTTATAGCGATAATAATACTGATGGGAGAGGGAAGGGCATGGATCAGGCGCTACTGGACGGGGG‑GTATCGCTGTTTTACCGGCGAAAAGATCGATGTc                                                                                      >  1:379443/1‑251 (MQ=255)
                                                                                                                                       gaatgTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGt                                                                                       <  1:196603/250‑1 (MQ=255)
                                                                                                                                               catcaGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTc                                                                                <  2:524868/249‑1 (MQ=255)
                                                                                                                                                           tgtgtgACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGatat                                                                   <  1:98293/250‑1 (MQ=255)
                                                                                                                                                             tgtgACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑TTCAAA‑AGGAGGAAGAATGGTTGGCTATCT‑AGGGGGAGGGA‑AT‑CTGGCCACCG‑AAAATGCCAGTGGTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTTTAGCGGTAATCATACTGATGAGAGAGGGAAGGGCATGGGTCAGGCGCTACTGGACGGGGG‑GTATCGCTGTTATAACGGCGCAAAAATCGATGTCTATTTTAAACCTGCGATATGTc                                                                >  1:531309/1‑251 (MQ=255)
                                                                                                                                                                gACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAA‑ATTAAA‑GGGAGTAAGAATGGTTTGCTATAT‑AGGGGGAGAGA‑CT‑CTGGCAACGG‑AAAATGCCAGTGCTGTATGAAG‑AATCC‑CGGG‑CTATGCTTATAGCGATAATCATACTGATGGGAGAGGGGAGGGCATGGATCAGGCGCTACTTGACGGGGG‑GTATCGCTGTTATACCGGCGAAAAGATCGGTGTCTATTTCAACACTGCGATATGGCAGc                                                             >  1:1828/1‑251 (MQ=255)
                                                                                                                                                                  cAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTTAAA‑GGGAGTAAAAATAATTTGCTATAT‑AGGGGGAGAGA‑CT‑CTGGCAACGG‑AAACTGCCCGTGCTGTCTGAAG‑ATTTC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGGTTA‑CGCTGTTTTTCCTGCGAAAAGAACGATGTCTATTTCAACCCCTCGATATGTCAGCAt                                                           >  1:336225/1‑251 (MQ=255)
                                                                                                                                                                  cAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGT‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATATCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATGTCAACACGGCGATATGTCAGCAt                                                           <  1:449450/251‑1 (MQ=255)
                                                                                                                                                                  aaaaTTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCAt                                                           <  2:200066/250‑1 (MQ=255)
                                                                                                                                                                           tAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGGAAGAATGATTGGCTATAT‑AGGGGGAGAGA‑CT‑CTGGCAACGG‑AAACTGCCCGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑GTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAAt                                                  >  1:524901/1‑251 (MQ=255)
                                                                                                                                                                               cATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCg                                              <  2:161579/251‑1 (MQ=255)
                                                                                                                                                                                         gatgaGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTTCCAGTGCTGTATGAAG‑AGTCC‑GGGG‑CGATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCATCATTCTGTCAATTGCGTACGTGGc                                      <  2:493048/249‑1 (MQ=255)
                                                                                                                                                                                                ggggTAA‑GTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAg                             <  2:311225/250‑1 (MQ=255)
                                                                                                                                                                                                  ggTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGtt                           <  1:143255/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGG‑CT‑CTGGCAACGG‑AAGCTGCCAGTGCTGTATGAAG‑ATGCC‑GGGG‑CTATTCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAAt                    <  2:504787/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAAt                    >  1:259327/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAA‑GGGAGTAAGAATGATTGGCTATAG‑AGGGGACGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGATGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑GTATCGCTGTTATGCCGGCGAAAAGATCTATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTaa                     <  1:539354/250‑1 (MQ=255)
                                                                                                                                                                                                                    tAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACg            >  1:36976/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGa           <  1:7254/250‑1 (MQ=255)
                                                                                                                                                                                                                        aaTGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGc       >  1:305627/1‑251 (MQ=255)
                                                                                                                                                                                                                              aggCTAGAT‑AAGGGAAGAGA‑CT‑CTGGCTACGT‑AAACTGCCAGTGCTGTATGAAG‑AATCC‑GGGG‑CGATGCTTATAGCGTTAATCATACTGGGGAGAGAGGGAAGGTCATGGATCAGGCGCTATTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAACGATCGATGTCTATTTCAACACTGCGATATGTCATCATTCTGGCAATTGCTTGCTGGGCAACGGCAAGTTATTTAATCTCAAACGaaatcagtgc  <  2:219250/249‑8 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
CGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAA‑TTCAAA‑GGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGA‑CT‑CTGGCAACGG‑AAACTGCCAGTGCTGTATGAAG‑ATTCC‑GGGG‑CTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGG‑TTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGCCGTGG  >  NC_000913/4351517‑4351986

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