Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,660,102 C→T N678N (AAC→AAT ynfE → putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,1020CT95.3% 122.1 / ‑2.2 43N678N (AAC→AATynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (23/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

TGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCC  >  NC_000913/1659891‑1660331
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
tGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGt                                                                                                                                                                                                >  1:495055/1‑251 (MQ=255)
         ataAGCGTAAAGATCCCACTGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCTGTAAGATTGAAATTTATTCCAGCTGGTTGGCGGAAAGTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACTTTTGAATGCTGGAATTCCCCTGAACGTAGAACTTGCCCACTGCATCTGGTTGGTTTCCAtt                                                                                                                                                                                         <  2:490919/249‑1 (MQ=255)
         ataAGCGTAAAGATCCCAATGGCCAATTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTGAAAAAGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCAGAAATTGCCCGTACCTGGTAACTGGAAAAAGATGAAGTGATAAGCCCAGTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATGCCCCTGAACGTAGAACCTCCCCACTGCAACTGTTTGGTTTCCATTa                                                                                                                                                                                        <  2:273109/250‑1 (MQ=255)
              cGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTCCAAATc                                                                                                                                                                                  >  2:199800/1‑251 (MQ=255)
                         cAATGGCCAGTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAAACCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTTCCAGTTTATGCCTCAACCTTTGAAGGCTGGATCTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTactcact                                                                                                                                                                       <  2:495055/251‑1 (MQ=255)
                            tGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGAAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAACAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACAGCAACTGTTTGGTTTCCACTACAAATCCCGTACTCACTCgc                                                                                                                                                                    >  2:236678/1‑250 (MQ=255)
                             tgCCATATTGTCGCCTACATAGCAATTCGTGACGCCCCAGAAGCAACTCCACTTAAAAAGACTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTAGCGGAAATCGCCCGTACCTGGGAACTGGAAAAAGGTGAAGTGATAAACCCATTGCCAGTTTATGCCTCAAACTTTGAAGGCTGGGATTACTCTGAACGACGAACCTTACCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGAc                                                                                                                                                                   >  2:255201/2‑251 (MQ=255)
                                cATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTgg                                                                                                                                                                                                 >  2:335939/1‑218 (MQ=255)
                                cATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTgg                                                                                                                                                                                                 <  1:335939/218‑1 (MQ=255)
                                             aCAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTg                                                                                                                                                     <  1:252830/249‑1 (MQ=255)
                                                   cATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTg                                                                                                                                             >  1:453883/1‑251 (MQ=255)
                                                                   cGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTc                                                                                                                             >  1:430543/1‑251 (MQ=255)
                                                                              cACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGa                                                                                                                                                                                                                               >  2:476464/1‑142 (MQ=255)
                                                                              cACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGa                                                                                                                                                                                                                               <  1:476464/142‑1 (MQ=255)
                                                                                    aaaCGCCTTCCGGTAAGATTGAAATTTATTACAGCAGGCTGGCGGAACTTGCCCGTACCTGGGAACTGGAAAAAGTTGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCCACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAc                                                                                                            >  1:533172/1‑251 (MQ=255)
                                                                                             ccGGTAAGATTGAAATTTAATCCAGCAGGCTGGCGGAAATTGCCCGAACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCATTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGACCCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGa                                                                                                    <  2:361664/250‑1 (MQ=255)
                                                                                                        gAAATTTATTACAGCAGGCTGGCGGGAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTGCTAACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTCGCTCAACCCTATAGTTTCGCAGAAACg                                                                                        >  2:466419/1‑251 (MQ=255)
                                                                                                               agtCCAGCAGGCTGGCGGAAATGGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCGCAACATTTGAAGGCGGGAATTCCCCTGAACGGAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGGCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGt                                                                                       <  2:29378/243‑1 (MQ=255)
                                                                                                               aTTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGt                                                                                       >  1:29378/1‑245 (MQ=255)
                                                                                                                     gcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAAc                                                                           <  1:249745/251‑1 (MQ=255)
                                                                                                                      cagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGGTAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGCGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACg                                                                          >  2:175049/1‑251 (MQ=255)
                                                                                                                      cagGCTGGCGGAAATTGCCCGTACCTGAGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCAACCTACGGCAATAtt                                                                                                                                                      >  1:155935/1‑175 (MQ=255)
                                                                                                                        ggCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGGGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTGCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGAGATTGCCAACGGc                                                                        <  2:430543/251‑1 (MQ=255)
                                                                                                                               ggAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACATATCCCGTACTCACTCGACCTACGGCAATATTGATCTCATGAAGGCTGCATGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATgg                                                                 >  2:147967/1‑251 (MQ=255)
                                                                                                                                      gCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGgt                                                          <  1:498880/251‑1 (MQ=255)
                                                                                                                                      gCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGgt                                                          <  1:393399/251‑1 (MQ=255)
                                                                                                                                         cGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                       >  1:563873/1‑251 (MQ=255)
                                                                                                                                              cTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTGTATGCCTCAACCTTTGAAGGCTGGCATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAcc                                                  >  1:321387/1‑251 (MQ=255)
                                                                                                                                                   aaCTGGAAAAAGCTGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGt                                             >  1:383690/1‑251 (MQ=255)
                                                                                                                                                        gAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCg                                         >  1:496269/1‑250 (MQ=255)
                                                                                                                                                           aaaGATGCAGTGATAAACCCATTGACAGTTTATGCCTGAACCTTTGAAGGATGGAATACCCCTGAACGTAGAACCTTCCCACTGCAAATGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGAATTATTTATCTCCTGAAGGCTGCCTGCCGACAGGCGGAGTGGAGCAAACCTATAGATGCGCAGAAACGTGGACTTGCCATCGTCGATATGGTGCGGGTGATTAACACTCGTGGCGCAGt                                     >  2:558431/1‑251 (MQ=255)
                                                                                                                                                            acgATGAAGCGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATCGATGCGCAGACACGTGGGATTGCCAACGGCGATATGGTGCGGGAGTTTAACCATCGTGGCGCAGtt                                    >  2:310201/3‑251 (MQ=255)
                                                                                                                                                             aGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATgg                                                                 <  1:500039/221‑1 (MQ=255)
                                                                                                                                                               aTGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCgg                                 >  1:183326/1‑251 (MQ=255)
                                                                                                                                                                     gggaaatcccatgtccAGTTTATCCCCCACACTTTGAAGTCTGGAATTACCCTGAACGTATTACTCTCCCCCTCCAACTGTTTGGTTTCCAATACCAACCCCGTACTCCCACGAGCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGCCCGGAGGGGCGGATCAACCCTATAGAAGCGCAGAAACGTGGTATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTCGCGAAGTTCGGCTagga                           <  2:533172/237‑4 (MQ=255)
                                                                                                                                                                      gATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGGTTAACCATCGTGGCGAAGTTCGGCTACCa                           <  1:199800/250‑1 (MQ=255)
                                                                                                                                                                             ccATTGCCATTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCgg                                 <  2:235418/237‑1 (MQ=255)
                                                                                                                                                                             ccATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCgg                                 >  1:235418/1‑237 (MQ=255)
                                                                                                                                                                             ccATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGt                   >  1:76387/1‑251 (MQ=255)
                                                                                                                                                                               aTTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTaa                 <  1:484567/251‑1 (MQ=255)
                                                                                                                                                                                     aGTTTATGCCTCATCCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAAACCCTTACTCACTCGACCTACGGCAATAGCGATCTCCTTAAGGCGGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGAAGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCACCGTGGCGAAGTTCGGCTACCAGCCAAAGTAAcaccac           <  2:76387/251‑1 (MQ=255)
                                                                                                                                                                                            gCCTCAACCTTTGAAGGCTGGAATTCCCCTGACCGTAGAACCTTCCCACTGCAAATGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGGGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGAGCGGGTGTTTAACCATCGTGGCGCAGTTCGGCTACCAGCCAAAGTAACACCACGTATtct    >  1:605209/1‑251 (MQ=255)
                                                                                                                                                                                                taacttattaatgCTGGAATTCCCCTGAACGTCGAACTTACCCACTGCAACTGTTTGTTTTCCATTACAAATCCCGTACTCACTCTACCTACGGCAATATTGATTTCCTGAAGGCTGCCCGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGTGTGTTTATCCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTactcgca  <  2:370645/237‑7 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
TGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCC  >  NC_000913/1659891‑1660331

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