Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,239 G→T A74E (GCG→GAG)  yaaA ← DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,2390GT89.8% 110.6 / ‑2.5 52A74E (GCG→GAG) yaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (29/15);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/6003‑6481
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
tgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGGGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCTTGGAATCt                                                                                                                                                                                                                                         >  2:609999/1‑251 (MQ=255)
  gTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTgg                                                                                                                                                                                                                                       >  1:580282/1‑251 (MQ=255)
    gctgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTTGCt                                                                                                                                                                                                                                     >  2:45848/3‑251 (MQ=255)
        atatCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCa                                                                                                                                                                                                                                <  1:78756/251‑1 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTc                                                                                                                                                                                                                             >  1:443623/1‑251 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTc                                                                                                                                                                                                                             >  1:83315/1‑251 (MQ=255)
           tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCAGGCT‑GCCCGTc                                                                                                                                                                                                                             >  2:385562/1‑251 (MQ=255)
             cccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCAt                                                                                                                                                                                                                           <  1:497737/250‑1 (MQ=255)
              ccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCa                                                                                                                                                                                                                            <  1:30001/249‑1 (MQ=255)
                     aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAagcg                                                                                                                                                                                                                    >  1:463931/1‑250 (MQ=255)
                       tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAagcgagc                                                                                                                                                                                                                 >  1:145433/1‑251 (MQ=255)
                         aTACAGATCTTTCCCTCGGGCATTCGTAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCTGGAAAGCATT‑CGCAAATGCGGTTGGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGTCTGC‑AAGCCGGTGGAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCGGGCT‑GCCAGTCATGAAAGCGAgcgg                                                                                                                                                                                                               <  2:618537/251‑1 (MQ=255)
                           aCAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATATGGCTGCCTTAAATCGAGCGGGCGCAGTCCGCCATCCACGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGCAATCGTCTTCGATGAAGGTTTCGGCCTGC‑ACGCCGGTGTAGCCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCCGTAAGGAAAGCGAgcggcg                                                                                                                                                                                                             >  2:305372/1‑251 (MQ=255)
                                   ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑g                                                                                                                                                                                                         <  1:49905/246‑1 (MQ=255)
                                    ccccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGGCATGAAAGCGAGCGGCGTT‑GATAccc                                                                                                                                                                                                   <  2:114642/250‑1 (MQ=255)
                                                           aTCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCAt                                                                                                                                                                            >  1:294486/1‑251 (MQ=255)
                                                            tCCCCATTTCCAGACGATAAGGATGCCTTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGAAAATGCTGTTGGGCACAATCGAAATCGTCTTCGATGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCCTGCAGTCTGGCT‑GCCAGTCCTGCAAGGGGGCGGCGTTTGA‑ACCCGCCAGtttttggcttatgagcat                                                                                                                                                                            >  2:436236/1‑235 (MQ=255)
                                                                cattgcaaGACGATAAGGCGGCATAAAATCGAGCGGGCGCAGTACCCCATACAAGCCGGAAAGTATT‑CGCAAATGCGGTTGGGCAAAAGCGAAATCGTCTGCGCGTAAGTTTTCGGCCTGC‑AAGCCGGTGTATACATCACCTTTAAACGCCATAATCGCCTGGCGGGCATTCTCCGTCGTGAAATCTGGCG‑GCCAGTCAGGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                                                                                                                                                       <  1:436236/244‑1 (MQ=255)
                                                                 aTTTCCAGACGATAAGAATGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGAAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTACGTAGTGTAATCTGGCT‑GCCAGTCATTAAAGCGCGCGGCATT‑GATACCCGCCAGTTTGACGCTGATGCGCATccgatt                                                                                                                                                                      >  2:44535/1‑246 (MQ=255)
                                                                         gcaaTAAGGCAGAATTAAAGCGAGCGGGCGCAGTACCCCAACCAGGCCGGAA‑GCCTTGCGAAAATGCTGTTGGGAAAAATCGAAATCTTCTTCGCTGAGGGTTCCGGCCGGA‑AAGCCGGTGTAGACATCACGTTTAACCGCCAGAATCGCCTGGCGGGCATTCCCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGGGAGGGGCGTT‑AATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCt                                                                                                                                                             <  2:145433/248‑1 (MQ=255)
                                                                                     cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCAACTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGAAGCGCATCAGCGTG‑CTCATCTGCGGCGGCGTc                                                                                                                                                  >  1:99905/1‑250 (MQ=255)
                                                                                          aaTCGAACGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACACCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGGCATGAAAGAGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAACCTGCGGAGGCGTCAGTTTc                                                                                                                                            >  1:291768/1‑251 (MQ=255)
                                                                                              gagtgggcgcAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCTGCGttt                                                                                                                                                                                                          <  2:636121/183‑2 (MQ=255)
                                                                                              gAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑g                                                                                                                                                                                                         >  1:636121/1‑187 (MQ=255)
                                                                                              aaGCGGACGCAGTACGCCATACAAAACGGAAAGCATT‑CGCAAATGCTGTTGGGCACAATCGAAATCGTCTTCGCTGAAGGTTTAGGCCTGC‑AAGCCGGTGTAGACATCACCTTTACACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCGTGAAATCTGGCG‑GCACGTCCAGAAAGCTCGCGGATTT‑GATCCCCGCAAGTTTGTGTCTTGTGCG‑ATAAGCATGGCTAATCTGCGGAGGGGTCAGTTTccacg                                                                                                                                        >  2:501103/2‑248 (MQ=255)
                                                                                                 ggggcgcAGTCCGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCgcgctt                                                                                                                                     <  1:45848/250‑3 (MQ=255)
                                                                                                              cctcaacccccctcttccattctgCAAATGTTGTTGGGCAAAACCAAAACCGTCTTCGCTGAAGGTTTCGGC‑TGCCAAGCCGTTGTAGACATCACCTTAAACCGCCAGAATCGCCTGGCGGGAATTCCCCGGCGTAAAATTTGGCT‑GCCAGTCATAAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                         <  2:386061/228‑1 (MQ=255)
                                                                                                               cATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgc                                                                                                                       <  1:659388/251‑1 (MQ=255)
                                                                                                                aTACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgct                                                                                                                      >  1:153083/1‑251 (MQ=255)
                                                                                                                        cGGAAAGCATT‑CGCAAATGATGTTGGGCAAAATCGAAATCGTCTTAGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCCGTCATGAAGGAGAACGGCGTT‑GATACCCGCCAATTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAAGTTCCGCGCTTCATGGATAAACGTCTGGGAATTg                                                                                                              >  2:676507/1‑251 (MQ=255)
                                                                                                                           aaaGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                           <  1:166533/251‑1 (MQ=255)
                                                                                                                            aaGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                           >  1:615531/1‑250 (MQ=255)
                                                                                                                             aGCATT‑CGCAACTGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCAGTGTAGACATCACCTTTAAACGCAAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGGTTACGCGCTTCATGGATCAACTGCTGGGAATTGTCTa                                                                                                          >  1:372901/1‑250 (MQ=255)
                                                                                                                                aTT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAATCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAAc                                                                                                        >  1:631213/1‑249 (MQ=255)
                                                                                                                                    agCAAATGCTGTTAGGCAAAATGGAAATAGTCTTCGCTGAAGGTCTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGACAAAATCGCATGGCGGGCATTCTACGGCGTGACATATGGCT‑GCAAGTCATGAAAGAGAGGGACGTT‑GATACCCACAATATTGACGCTGATGCGGAACAACGTG‑CATATATGAGGAGGAGTCGGGTTACTCGCTTCATGGCACAACGGCTGGGAATGTTCTAACAGCTc                                                                                                   >  2:652653/2‑251 (MQ=255)
                                                                                                                                     gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                         >  1:386061/1‑228 (MQ=255)
                                                                                                                                     gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTcc                                                                                                  >  1:396470/1‑251 (MQ=255)
                                                                                                                                        aaTGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACACCCCCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGttt                                                                                                                                                                                           <  1:374004/160‑1 (MQ=255)
                                                                                                                                          tGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGACAGCGAGCGGCGTT‑GATACCCGCCAGTTTTTCGCTGATGCGCATCCGAGTG‑CTACACTGCGGAGGCGTAAGTTTCAGCGCTTCCTGGATCAACTGATGGGAATTGTTTACCAGCTCCGGccg                                                                                             >  2:634000/1‑249 (MQ=255)
                                                                                                                                                    cAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcgc                                                                                   >  1:466424/1‑251 (MQ=255)
                                                                                                                                                     aaaaTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGcgcg                                                                                  >  1:494217/1‑251 (MQ=255)
                                                                                                                                                         ccgaaatcgTCTCAGCTGAAGGTTCCGGCCTGA‑AAGCCGGTGTAGACATCACTTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCgtggt                                                                              <  1:664721/250‑1 (MQ=255)
                                                                                                                                                            aaatcgTCGTCACTGAAGGTTTCGGACTGC‑AAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGGATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGGGCGGCGTT‑GATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTG‑CTAATATGCGGAGGCGTGAGTTTCCGAGCTTCAAGGATCAAATGAGGGGCATTATCTGACAGCTCGGGCAGCGTAAAGCGCGtggtggg                                                                           >  2:333223/1‑250 (MQ=255)
                                                                                                                                                                     tCGCGGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGCCCTCACCTTTAAACGCCAGAATCGCCTGGCGGCCATTCTCCGGCGTGAAATCGGGCT‑GCAAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACgg                                                                     <  1:373925/248‑1 (MQ=255)
                                                                                                                                                                       gCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCttt                                                                >  1:434204/1‑251 (MQ=255)
                                                                                                                                                                           atggTTGCGGCCTGC‑AAGCCGTCGGAGAAACCCCCTTAAAACGCCAGAATTGCGTGCCGGGCACTCGCCGGCGATAAGTATGGCT‑GCATGTCATGAAAGCGAGTGTGGTT‑GATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGGGGTGGTCAACGTGCTTTGGTa                                                            <  2:631213/249‑1 (MQ=255)
                                                                                                                                                                               tttCGGCCTGC‑AGGCCGGTGGAGACATCACCTTTAAACGCAAGAACCGCCTGGCGGGCATGCTCCGGCGTGAAATCGGGCT‑GCCGGTCATGAAAGCGAGGGGGGTT‑GATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATGGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATc                                                         <  1:305372/250‑1 (MQ=255)
                                                                                                                                                                                   cggctggaaaGCCGGTGTAGACATCACCTTAAAACGCCATAAGCCTCTGGCGGGCATTCTCCGGTGTGAAATCTGGCT‑GCCCGTCATGAAAGCGGGCGGGGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGc                                                     <  1:126040/243‑1 (MQ=255)
                                                                                                                                                                                              ccGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAgg                                          >  1:109737/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccgggCGGGCATTCTCTGGCGGTAAATCTGGCT‑GCCAGTCAGGAAAGCGAGCGGTGTT‑GATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAgg        <  1:634000/248‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGGGCATTCGCCGGCGTGAAGTC‑GGCTCGCCAGTCATGAAAGCGGGGGGCGTT‑TACACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCCGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGa       <  1:263979/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gggCATTCCCCGGCGGAAAATCTGGCT‑GCCAGTTATAAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt  <  1:518608/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
TGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/6003‑6481

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