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.

CGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGCCGTGG  >  NC_000913/4351517‑4351986
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
cgcgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAAc                                                                                                                                                                                                                                   <  1:148107/250‑1 (MQ=255)
          cACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCa                                                                                                                                                                                                                       <  2:571246/251‑1 (MQ=255)
              gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAAATTCAAGGGGGGAAAGATGGTTTGCTTTTTAGGGGGAGAGAACTCTGGCAAGGGaaatggcagggc                                                                                                                                                                                                                    >  1:78044/1‑242 (MQ=255)
              gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGc                                                                                                                                                                                                                   <  2:132917/251‑1 (MQ=255)
                        aaGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGGAGAGAATCTGGCAACGGAAAATGCCAGTGCTGTATGa                                                                                                                                                                                                            >  1:524868/1‑248 (MQ=255)
                          gTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGTTACTCAAGATCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCGCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGAtt                                                                                                                                                                                                      <  2:548320/251‑1 (MQ=255)
                                          accaATAGCACCGAGAACAGTACCGCACTGACAATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTtat                                                                                                                                                                                    <  2:365403/251‑1 (MQ=255)
                                             aaTAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGc                                                                                                                                                                                 <  1:23217/251‑1 (MQ=255)
                                                 gCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAAGGGGGTAGAATGGTTTGCGATATATGGGGAGGGAACATGGCAACGGGAAACTCCAGTGCCGTATGAAGATTCCCGGGGCACGCTTATAGCGAtt                                                                                                                                                                             >  2:96656/1‑250 (MQ=255)
                                                   aCCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAAt                                                                                                                                                                           <  2:81548/251‑1 (MQ=255)
                                                   aCCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAAGTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCATCGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAAt                                                                                                                                                                           <  2:188086/251‑1 (MQ=255)
                                                                 cgctcTGACCATTAAGATATCTGTTGTACTCAATGTCCACGGACGGTTTCGTAACATGCGGTAGGTCAATTAATGTCTCATCAGCGTCCTGGGGTGACAAATTTCTTAAGCGTTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGAGTCCGGGGCTATGCTTATAGCTATAATCATACTGATgagag                                                                                                                                                             <  2:180196/247‑1 (MQ=255)
                                                                 cGCACTGACCATTAAGATCACTGTGGTACTCAATTTCAGCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTAGGCTTATAGCGATAATCATACTGATgagag                                                                                                                                                             <  2:513859/251‑1 (MQ=255)
                                                                       gACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGAGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAATCAAAGGGGGGAAGAATGATTTGGTATATAGGGGGAGAAGCCCTGGCAAAGGGAAATGCCAGAGCTGTAAGCAGAATCCGGGGGTTATGCTATAGGGAGAATTATACTGGTGGGAGAggg                                                                                                                                                         >  1:37290/1‑249 (MQ=255)
                                                                        aCCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAAGGAGGAAGAATGATTGGCTATATAGGGGGAGAGACTCTGGGAACGGAAACATCCAGGGCTGGATGAAGGTTCCGGGGGTACGCTTATAGCGGTAATTCTAATGTTGAGAGAgggg                                                                                                                                                        >  2:49304/1‑248 (MQ=255)
                                                                           aTTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTc                                                                                                                                                   <  1:187012/251‑1 (MQ=255)
                                                                           aTTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTc                                                                                                                                                   <  1:578296/251‑1 (MQ=255)
                                                                                      tgtgGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAGTTCAAGGGAGTAAGAATGATTTGGTTTATAGGGGGAGGGACCCTGGCAACGGGAAATTACCAGGCTGTATGGAGGTTCCAGGGCTATGGTTATAGCGATAAATATACTGGTGGGAGAGGGAAGGGGATGGATccggg                                                                                                                                        >  2:438467/1‑247 (MQ=255)
                                                                                       ttgggaCTCAATGTCATCGGACGTTTGCGTAACATGCGGTAGTGCAATGAATGTCTCAGCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGGAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGgc                                                                                                                                        <  2:139748/245‑1 (MQ=255)
                                                                                           tACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCCTCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTTAAGGGAGTAAGAATGATTGTCTTTATAGGGGGAGGGACTACTGCACCGGAAAATGCCCGGGCTGTATGGAGAGTCCGCGGCTAAGCTTTTAACGATAATCATTCTTTTTAGAGAGGGGGGGTCATGGGTCTGGCGCtc                                                                                                                                    >  2:78265/1‑249 (MQ=255)
                                                                                                aaTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTTATGAGGCGGGTAATTCAAAGTGAGTAAGAATGTTTGGCTATATAGGGGAAGAGTCTCTGTCAACGGAAACTTCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTAGAGCGATACGCATGCTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTAATGGAc                                                                                                                              <  2:48825/251‑1 (MQ=255)
                                                                                                 aTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGGGTAAGAAAGATTTGCTTTTTTGGGGGAGAGACTCTGGCAACGGAAACTGCCCGTGCTGTATTAAGATTTCGGGGGTTTGCTTTTATCGATAATCATACTTTTGAGAGAGGGAGGGTCATGGTTCAGGCGCTACt                                                                                                                                  >  2:157377/1‑246 (MQ=255)
                                                                                                 aTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACt                                                                                                                                  <  1:157377/246‑1 (MQ=255)
                                                                                                         ggACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAATTAAAGGGGGGAAGAATGATTGGGTTTATAGGGGGAGAGACTCCGGCAACGGGAAATGCCCGTGCTGTATTAAGATTTCGGGGCTATGCTTTTAGAGGTAATCATACTGATGAGAGAGGGAAGGGCATTGGTCAGGGGCTACTGGACGGGGGTTAt                                                                                                                    >  1:287906/1‑251 (MQ=255)
                                                                                                          gACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTTGGACAAATTTCTTAAGCATTATCGCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAGCTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGTGTTATc                                                                                                                   <  2:419943/251‑1 (MQ=255)
                                                                                                           aCGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAAGGGGGAAGAATGATTGGCTTTTTAGGGGGAGAGAATTTGGCAACGGAAACTGCCCGGGCTGGTTGAAGATTTCCGGGCCATGCTTATAACGGGAATTATTATGATGAGAGAGGGAAGGGCATGGGTAAGGGGCTCCTCGAGGGGGGGTATCg                                                                                                                  >  2:159706/1‑251 (MQ=255)
                                                                                                                 gagTAACAGGCGGTAGGGCAATGAAGGACTCATCAGCTTCCTTTTGTGACAAATTGCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGTCGATATAGGGGAGGAGACTCTGGCAACGGAAGCTGCCAGTGCTGTATGAAGAGTCTGGGGCTATGCTTATAGCGATAATCATACTGATGATAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTta                                                                                                            <  2:310539/249‑1 (MQ=255)
                                                                                                                        aTGCGGTCGGGCAATGACTGTCTCTACATCTTCCTTGTGTGACAAATTGCTTAAGCATTATCTCAGATGAGGCGGGTAATTTAAAGGGATTAAGAATGATTGGTTATATAGGGGAAGAGACTCTGGCATCGGAAACTGCCAGTGCTGTATGAAGAGTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGATAGGTAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATAcctgc                                                                                                     <  2:127488/251‑4 (MQ=255)
                                                                                                                                       gaatgTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCCACGGGAGTAAGAATGATTGGGTATATCAGGGGAGAGACTCTTGTAACGGAAAATGCCCGAGGTGTATGAAGGATTCGGGGCTATGCTTATAGCGATAATAATACTGATGGGAGAGGGAAGGGCATGGATCAGGCGCTACTGGACGGGGGGTATCGCTGTTTTACCGGCGAAAAGATCGATGTc                                                                                      >  1:379443/1‑251 (MQ=255)
                                                                                                                                       gaatgTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGt                                                                                       <  1:196603/250‑1 (MQ=255)
                                                                                                                                               catcaGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTc                                                                                <  2:524868/249‑1 (MQ=255)
                                                                                                                                                           tgtgtgACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGatat                                                                   <  1:98293/250‑1 (MQ=255)
                                                                                                                                                             tgtgACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAAGGAGGAAGAATGGTTGGCTATCTAGGGGGAGGGAATCTGGCCACCGAAAATGCCAGTGGTGTATGAAGATTCCGGGGCTATGCTTTTAGCGGTAATCATACTGATGAGAGAGGGAAGGGCATGGGTCAGGCGCTACTGGACGGGGGGTATCGCTGTTATAACGGCGCAAAAATCGATGTCTATTTTAAACCTGCGATATGTc                                                                >  1:531309/1‑251 (MQ=255)
                                                                                                                                                                gACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAAATTAAAGGGAGTAAGAATGGTTTGCTATATAGGGGGAGAGACTCTGGCAACGGAAAATGCCAGTGCTGTATGAAGAATCCCGGGCTATGCTTATAGCGATAATCATACTGATGGGAGAGGGGAGGGCATGGATCAGGCGCTACTTGACGGGGGGTATCGCTGTTATACCGGCGAAAAGATCGGTGTCTATTTCAACACTGCGATATGGCAGc                                                             >  1:1828/1‑251 (MQ=255)
                                                                                                                                                                  cAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGAGTAAAAATAATTTGCTATATAGGGGGAGAGACTCTGGCAACGGAAACTGCCCGTGCTGTCTGAAGATTTCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGGTTACGCTGTTTTTCCTGCGAAAAGAACGATGTCTATTTCAACCCCTCGATATGTCAGCAt                                                           >  1:336225/1‑251 (MQ=255)
                                                                                                                                                                  cAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGTCTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATATCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATGTCAACACGGCGATATGTCAGCAt                                                           <  1:449450/251‑1 (MQ=255)
                                                                                                                                                                  aaaaTTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCAt                                                           <  2:200066/250‑1 (MQ=255)
                                                                                                                                                                           tAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGGAAGAATGATTGGCTATATAGGGGGAGAGACTCTGGCAACGGAAACTGCCCGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGGTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAAt                                                  >  1:524901/1‑251 (MQ=255)
                                                                                                                                                                               cATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCg                                              <  2:161579/251‑1 (MQ=255)
                                                                                                                                                                                         gatgaGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTTCCAGTGCTGTATGAAGAGTCCGGGGCGATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCATCATTCTGTCAATTGCGTACGTGGc                                      <  2:493048/249‑1 (MQ=255)
                                                                                                                                                                                                ggggTAAGTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAg                             <  2:311225/250‑1 (MQ=255)
                                                                                                                                                                                                  ggTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGtt                           <  1:143255/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGGCTCTGGCAACGGAAGCTGCCAGTGCTGTATGAAGATGCCGGGGCTATTCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAAt                    <  2:504787/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAAt                    >  1:259327/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAAGGGAGTAAGAATGATTGGCTATAGAGGGGACGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGATGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGGTATCGCTGTTATGCCGGCGAAAAGATCTATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTaa                     <  1:539354/250‑1 (MQ=255)
                                                                                                                                                                                                                    tAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACg            >  1:36976/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGa           <  1:7254/250‑1 (MQ=255)
                                                                                                                                                                                                                        aaTGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGc       >  1:305627/1‑251 (MQ=255)
                                                                                                                                                                                                                              aggCTAGATAAGGGAAGAGACTCTGGCTACGTAAACTGCCAGTGCTGTATGAAGAATCCGGGGCGATGCTTATAGCGTTAATCATACTGGGGAGAGAGGGAAGGTCATGGATCAGGCGCTATTGGACGGGGGTTATCGCTGTTATACCGGCGAAACGATCGATGTCTATTTCAACACTGCGATATGTCATCATTCTGGCAATTGCTTGCTGGGCAACGGCAAGTTATTTAATCTCAAACGaaatcagtgc  <  2:219250/249‑8 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
CGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGGAAGGTCATGGATCAGGCGCTACTGGACGGGGGTTATCGCTGTTATACCGGCGAAAAGATCGATGTCTATTTCAACACTGCGATATGTCAGCATTCTGGCAATTGCGTACGTGGCAACGGCAAGTTATTTAATCTCAAACGAAAGCCGTGG  >  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: