Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,958 G→A D64N (GAC→AAC)  yahH → putative uncharacterized protein YahH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,9580GA86.5% 75.9 / ‑4.5 37D64N (GAC→AAC) yahHputative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base A (11/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAG  >  NC_000913/339713‑340202
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
tCTTGGCTGCTGCGACAATTCCGGCGACTGACATACATGCTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACAGGCGCCCGAATATGTAGGCCTGATAAGCGTGGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTTGCCCCC‑TTATTCCGCCTCCGAATGGCGCACGAATCTGCAGGCCGGCTGAACGTAGCGCAT‑CAGGCATTTGTGCCCATTGACGGATGCAGCGTGGACGcc                                                                                                                                                                                                                                                  >  2:37508/1‑251 (MQ=2)
       tGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTGTCACCATTGCCGGTTGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCCAGAATCCTTAGGCCTGATAAGCGTAGCGCCT‑CCGGCATTTGTCACCACTGCCGGATGCGGCGTGGaacccctttccc                                                                                                                                                                                                                                           >  2:47332/1‑240 (MQ=17)
              aaCGGTACGGCGTCTGGCCTACATGTTCATGCCGGATGCGGCGTGATCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATGTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CATGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAGCCGGC‑CTACg                                                                                                                                                                                                                                   <  1:282909/251‑1 (MQ=14)
               aCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGC‑CTACg                                                                                                                                                                                                                                   >  1:97133/1‑250 (MQ=255)
                 gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGt                                                                                                                                                                                                                                <  1:260671/251‑1 (MQ=255)
                    aCGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGAGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACCCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCAT‑CAGGAATTTGTCACCATTGCCGGATGCGGCGTGAACGACTTATCCGGC‑CTACGAGTGgc                                                                                                                                                                                                                             >  2:442757/1‑251 (MQ=14)
                               ccTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAttgt                                                                                                                                                                                                                  <  1:192464/251‑4 (MQ=255)
                                 tACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGACTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGACAGCC‑TTATCCGTCCTACGAATGGCGCAAGAATTTGTAGGCCTGATAAGCGGAGCGCATTCAG‑CATTTGTCACCATTGACTGATGCGGCGGGGACCCCTTCTCCggtccacgcgtggggggagcatgtt                                                                                                                                                                                                                   >  1:235567/1‑226 (MQ=14)
                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAACGTAGCGCCT‑CAGGCATTTGTCCCCCTTTCCGGATTCGGCGGTGACCCCTTTTCCCGC‑CTCCGTGTGGCGCGGGAAT‑TTTAGGCCTGATaacc                                                                                                                                                                                                    >  1:236759/1‑249 (MQ=11)
                                                     ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCAt                                                                                                                                                                                           <  2:441825/251‑1 (MQ=9)
                                                        gtgtaCGCCTTATCCGGCCTTAGAATGGCGCAAGAATCTGTAGGCCCGAAAAGCGTTGCGGATCCGGCATTTGTCACCCTTGGCGGATGCGGCGGGAACGCCCTTA‑CCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                                                                        <  2:445827/247‑1 (MQ=11)
                                                              gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCAGGATGCGGCGTGAACACA‑TTATACGGCCTACGAATGGCGCAAGCATCTGTAGGCCTGATAAGCGTAGAGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGAGTGAACGCCTTATACGGC‑CTACGCGGGGCGCGAGAGT‑TGTAAGCCTGATAAGCAGAGCGCATCAGGCaga                                                                                                                                                                                   >  1:267241/1‑248 (MQ=11)
                                                                        gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTCCGAGTGGCGGGAGGAT‑TGTAGGCCGGATAAGCGTAccgcctccggcggactggcgttggg                                                                                                                                                                        >  1:191240/1‑230 (MQ=11)
                                                                                                                                              aTGCGGCGTGCACGCC‑TTCTCAGTCTTAGGAATGTAGAAAGAATCTGGAGGCGGGATAAGCGGAGCTGGT‑CAGGCAGTTGTCTTCACTGCCGGATGGGGCGTGAACGCCTGATCCGGC‑CTACGAATGGCGCAAGAATCGGTAGGCCGGATAAGCGTAGCGCAGCAGGCAGTTGGCTCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGAGTGGCGCGATAATTGTAGGTCTGATAAGCGTAgc                                                                                                   <  2:364572/251‑1 (MQ=2)
                                                                                                                                                                tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                                    <  1:47332/249‑1 (MQ=2)
                                                                                                                                                                tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTg                                                                                    <  2:267701/249‑1 (MQ=2)
                                                                                                                                                                      gCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGCGTTGGt                                                                            >  1:46632/1‑251 (MQ=11)
                                                                                                                                                                       ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTc                                                                           <  1:325921/251‑1 (MQ=255)
                                                                                                                                                                         tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCAt                                                                         >  1:199900/1‑251 (MQ=255)
                                                                                                                                                                           cGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                       <  1:89022/251‑1 (MQ=255)
                                                                                                                                                                                       gaagcTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCct                                                           <  2:339191/247‑1 (MQ=255)
                                                                                                                                                                                         aTCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGGGAGAATTGTAGGGCTGATAAGCGGTGCGCATCAGGCAGTCTGGCGTTGGTCATAAccccctcccccccc                                                         >  1:407299/1‑241 (MQ=11)
                                                                                                                                                                                             gTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTa                                                      <  1:263158/250‑1 (MQ=255)
                                                                                                                                                                                                 gCCTGATAAGCGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGa                                                 <  2:222080/251‑1 (MQ=255)
                                                                                                                                                                                                           cGTAGCGCAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCCGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaac                                        >  1:254881/1‑250 (MQ=255)
                                                                                                                                                                                                                 gcAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGAGTGGCGCGAGTATTGTTGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGGGTTGGTCATAACCCCATCACCCCCTGTAGCGGACATaacaactttaa                                   >  1:35059/1‑244 (MQ=14)
                                                                                                                                                                                                                  cAT‑CAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAaga                                <  1:445370/251‑1 (MQ=255)
                                                                                                                                                                                                                      cAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTGCGAATGGCGCAAGAATCTTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGGTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCt                             <  1:434886/251‑1 (MQ=17)
                                                                                                                                                                                                                          cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg                         <  1:73791/251‑1 (MQ=255)
                                                                                                                                                                                                                          cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGGGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg                         >  1:348467/1‑251 (MQ=255)
                                                                                                                                                                                                                             ttgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGt                       <  2:50824/246‑1 (MQ=255)
                                                                                                                                                                                                                              tgtcaccatTGCCGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTc                      <  1:173061/247‑1 (MQ=255)
                                                                                                                                                                                                                                         ccGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACGAGTGGCGCGAGAATTGTAGGCCCGATAAGCGTAGCGCCTCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGAAATAACAACATGAAAGAGGTTGTGGTCGGTGCCATtggg          >  2:134932/1‑250 (MQ=18)
                                                                                                                                                                                                                                          cGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCTCCATTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCGACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTGTGGCGTTGGTCATAAGCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTGGTGGTCGTTGCCAAtggtg         <  2:31426/251‑1 (MQ=18)
                                                                                                                                                                                                                                          cGGATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAAAAACATGAAAGAGCTTGTGGTCGTTGCCATttgtg         >  2:205054/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gATGCGGCGTGAACGCCTTATCCGGC‑CTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGc       <  2:302776/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggcgTGAACTCCTTATCCGGCTCT‑CGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGCGCGAGAATTGTATGCCTGTTAAGCGTAGCGCAGCAGGCAGTCTGGCGTTGGTCATAGCCCCATCCCCCGCTGGAGCGGACATAACAACATGAAAGAGCGTCTGGGCTTTGCCATTGGTGGCAACAg  <  2:438238/251‑1 (MQ=17)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
TCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCC‑TTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGC‑CTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAG  >  NC_000913/339713‑340202

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