Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,3730GT34.2% 53.1 / 11.7 39G33C (GGC→TGC) yedEputative selenium transporter YedE
Reads supporting (aligned to +/- strand):  ref base G (8/17);  new base T (12/1);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 5.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCA  >  NC_000913/2008144‑2008622
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                             
gCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTCTTTTTGCATTACTCGCACCtttttgg                                                                                                                                                                                                                                        >  1:513616/1‑248 (MQ=255)
                     gagaAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGTCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTGGGCATTACTGGCACCTTTTGGTCTGTCACGGGTGAATTTAccc                                                                                                                                                                                                                   <  1:359329/251‑1 (MQ=255)
                        aaaTTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAACATAAATTTGCCTGGTGTATGCTTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGAATTTTCTCTACCTTCTCTTTTTGCATTACTGGCACCTTTTGGGCTGTCACCGGGGAATTTACCCGtt                                                                                                                                                                                                                >  1:658689/1‑251 (MQ=255)
                         aaTTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCCGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTg                                                                                                                                                                                                               <  1:441077/251‑1 (MQ=255)
                               agcTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGggcgg                                                                                                                                                                                                         <  2:400407/249‑1 (MQ=255)
                                   aaCTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTCCTTTCTATTTTTGCCTTTCTGGCACCTTTTGGGCTGTCAAGGGGGAATTTACCCGTTGGGGCGGcc                                                                                                                                                                                                       >  1:350806/1‑249 (MQ=255)
                                       attattTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCCCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCt                                                                                                                                                                                                   <  1:235323/249‑1 (MQ=255)
                                               tgtgAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTGGGCTGTCCAGGGGGAATTTTACCGTTTGGGCGGCCCGCCCCTGCAACTg                                                                                                                                                                                         >  1:290366/1‑251 (MQ=255)
                                                  gAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTATTTTTGCATTAATGGCACCTTTTTGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCCGCTCCTGCCACTGttt                                                                                                                                                                                      >  1:576073/1‑250 (MQ=255)
                                                               gAAAATAAATTTGCATGGCGTATGCTTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCCAATTACTATTTTCGGATTACTGGCACCCTTTTGGCTGTCACCGGAGAATTTACCCGCTGGGGCGGCCAGCTCCTGCCCCTGTTCGGCGTACctgctg                                                                                                                                                                       >  1:635000/1‑251 (MQ=255)
                                                                  aataaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTGACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATGTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAag                                                                                                                                                                    <  1:680139/251‑1 (MQ=255)
                                                                     aaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAagag                                                                                                                                                                  <  1:167941/250‑1 (MQ=255)
                                                                       atTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATTGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTgg                                                                                                                                                               <  2:576073/251‑1 (MQ=255)
                                                                                  gtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTa                                                                                                                                                     <  1:71799/250‑1 (MQ=255)
                                                                                  gtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTTTCTAATTACTATTTTTGCATTACTTGCACCCTTTTGGGTTTCCCGGGGGGATTTACCCGTTGGGGGGGCCCGGTCTTCCTACTGTTCGGAGTccctggtgcagggggggggttctttc                                                                                                                                                       >  2:608639/1‑226 (MQ=255)
                                                                                    aTGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTTCTTACTATTTTGGCATTAATGGCACCTTTTTGGCTGTCACGGGGGAATTTTCCCGTTGGGGAGGCCCGCTCCTGCAACTGTTCGGCGGCCCTGCTGAAGGGTgggggttcttttaa                                                                                                                                                   >  1:261586/1‑240 (MQ=255)
                                                                                           gcgACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGAAATTTTTCTAATTTCTCTTTTGGCATTACTAGCACCTTTTTGGCTGTCACGGGGGAATTTTTCCGTTGGGGGGGCCCGCTCCTGCCACTGTTCGGCGTCCATTCTGAAGAGTGGGGGTACTTTTAAATTATccc                                                                                                                                           >  1:458175/1‑250 (MQ=255)
                                                                                                tctTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTTGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCCTGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTgg                                                                                                                                      >  1:660340/1‑251 (MQ=255)
                                                                                                   tGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGTGCCCCAGCCCTTCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACTGGTGAATTTACCCGTTGGGGCGGCTAGCTCCTGCAACTGTTCTGCGGCCTATCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGaa                                                                                                                                    <  1:669351/250‑1 (MQ=255)
                                                                                                     tttCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTAATTTAAAATTATCCCACTGGAAgg                                                                                                                                  >  1:598795/1‑250 (MQ=255)
                                                                                                       tCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGAt                                                                                                                                <  1:184203/250‑1 (MQ=255)
                                                                                                        cAAACCTTCGGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATGCTCTCGACTTACTATTTTGTCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTGTAAAATTATCCATCTGGAAGGATc                                                                                                                               <  2:388950/250‑1 (MQ=255)
                                                                                                          aaCCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTTCCTACTCTTTTTGCATTCCTGGCACCCTTTTGGCTGTCACCGGGGAATTTTCCCCCTTGGGCGGCCAGCCCCCGCAACTGTTTGGCGTCCATGCCTGAGAGGGGGGTTTCTTTTAAATTTTCCATCTGGAAGGGTCGcc                                                                                                                            >  2:657248/1‑251 (MQ=255)
                                                                                                               tCTGAGTTCAGAGGCTAAATTAATGTCATGGCACCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCAATCCCTGCGGTCATCGCGGCGGGCCTTCTCTCTTCTTACTATTTTTGCATTACTGGCACCTTTTTGGCTGTCCCGGGTGGATTTTCCCGTTGGGGCGGGCAGCTCCTGCAACTGGTCGGGGTTCATGCTGAAGAGGGGGGGTATTTTTAAATTTTCCATCTGGGAGGGTCGCCATTa                                                                                                                        >  2:433918/1‑250 (MQ=255)
                                                                                                                       cAGAGGCTACATTCATGTCATGGCGCCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTCGGGGCGGCCAGCTCCTGCTACTGTTCGGCGGCCATGCTGAAGAGTGGGGTGACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCAt                                                                                                                <  1:379305/250‑1 (MQ=255)
                                                                                                                        agagGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCg                                                                                                              <  1:603831/251‑1 (MQ=255)
                                                                                                                            gCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAATTTCTGGGCCCCCATCCCTGAGGTCATCGCGGCGGGGGTTTTTTCTCCTTCCCTTTTTGGCAATTATGGGACCTTTTTGGTTGTCACCGGGGAATTTACCCGTTGGGGCGGCAAGCTCCTGCAAACGTTCGGCGTCCATGCTGAAGGGGGGGGTTACTTTAAAATTATCCATCTGGGAGGGTCGCCATTAACCCGCCTTGACgg                                                                                                          >  1:39317/1‑251 (MQ=255)
                                                                                                                            gCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCGGCGGTCATCGCGGCGGGTATTCTCTCTACGTACTATTTTGGCATTACTGGCACCTTTGGGGCTGTGATGGGTg                                                                                                                                                                                                                            <  2:662341/139‑1 (MQ=255)
                                                                                                                            gCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCACTGCGGTCATCGCGGCGGGTATTCTCTCTACTTTCTATTTTGGGATTACTGGCACCTTTTTGGCTGTCACggggg                                                                                                                                                                                                                            >  1:662341/1‑137 (MQ=255)
                                                                                                                                               agcaATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTGGCACTAATGGCACCTTTTGGGCTGTCACGGGGGAATTTAACCGTTGGGGCGGCCCGCTCCTGCCACTGGTCGGCGGCCCTGCTGAAGAGGGGGGGTACTTTTAAAATATTCATCTGGAAGGATCGCCATTAACCCGCCTCGACGGGGTGATGATACTCGGGATg                                                                                       >  1:506538/1‑251 (MQ=255)
                                                                                                                                                       aaaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTACTTTTTTTGGATTTCTGGCACCTTTTTGGGTGTAACGGGTGAATTTACCCGTTGGGGCGGCCAGCCCCTGCAACTGTTTGGGGTCCATGCTGAAGAGGGGGGGTACTTTTAAATTATCCATATGGAAGGGTCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTggcgg                                                                               >  1:385873/1‑251 (MQ=255)
                                                                                                                                                        aaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCCTCGCGGCGGGCATTTTCTCTACTTAATATTTTGGGATTACTTGGACCTTTTTGGCTGTCACGGGGGGATTTACCCGTTTGGGCGGCCAGCTCCTGCCACCGTTTGGCGTCCCTTCTGAAGAGTGGGGGTACTTTTAAATTATCCATCTGGCAGGATAGCCATTAACCCCCAACGACGGGGTGCTGctcctccggctggtttgcgtg                                                                              >  2:94482/1‑236 (MQ=255)
                                                                                                                                                                  gTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTTCTATTTTTGCATTTCTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGGCAGCTCCTGCAACTGTTCGGGGTCCATGCTGAAGAGGGGGGGTTCTTTTAAATTATCCATCTGGGAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGGATGTTTTGCGGCCTCTTTGccc                                                                    >  1:222209/1‑250 (MQ=255)
                                                                                                                                                                                   cccccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCCCCCTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCcaac                                                   <  1:431909/251‑1 (MQ=255)
                                                                                                                                                                                       cagcccTGCGGTCATCGCGGCGGGTATTCTCTCTGCTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCGGCTTTGCCGCAGCGCTGTGGGCCAACAATg                                               <  2:634711/248‑1 (MQ=255)
                                                                                                                                                                                          cccTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACGATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTAGCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTc                                             <  1:58552/250‑1 (MQ=255)
                                                                                                                                                                                                        cggcgtGTCTTCTCTCTACTTTCTATTTTTGCATTACTGGCTCCTTTTGGGCTGTCACGGGTGAATTTTCCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAATAGTGGGGTTACTTTAAAATTATCCATCTGGAATGATCGCCATTAACCCGCATCGACGGTATGATGATCCTCGGTATGTTTGGTGGCTTCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGc                                <  1:723785/249‑1 (MQ=255)
                                                                                                                                                                                                         ggcggGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGAGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTCTGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCcgc                             >  1:260986/1‑251 (MQ=255)
                                                                                                                                                                                                                aTTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCg                      >  2:386560/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggCATTACTGGCACCTTTTGGGCTGTCCCGGGTGAATTTACCCGCTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAGCCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTAGCCGcg              >  2:54028/1‑237 (MQ=255)
                                                                                                                                                                                                                                     ggCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTGCCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCTCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCca  <  1:480885/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGCTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCa              <  1:54028/238‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                             
GCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCA  >  NC_000913/2008144‑2008622

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