Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,156,913 G→A G46D (GGT→GAT)  ycfH → putative metal‑dependent hydrolase YcfH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,156,9130GA97.4% 115.7 / NA 39G46D (GGT→GAT) ycfHputative metal‑dependent hydrolase YcfH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/19);  minor base C (1/0);  total (20/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 = 9.95e-01

GTTACAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGTGAACTGAACAAGCCGG  >  NC_000913/1156668‑1157155
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
gTTACAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGattatt                                                                                                                                                                                                                                               <  1:156458/251‑1 (MQ=255)
gTTACAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGattatt                                                                                                                                                                                                                                               <  1:203138/251‑1 (MQ=255)
   aCAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTAc                                                                                                                                                                                                                                            >  1:452570/1‑251 (MQ=255)
   aCAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGGCGCCACAACATTACCGGATTATTTa                                                                                                                                                                                                                                             <  2:553951/250‑1 (MQ=255)
      ggCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTCCCTCATCTTTAAGAGAGACATCATGTGTTTAGTCGACTCACACTGCCATCGCGATGGTCTGGATTATGAATCTTTGCATAAGTACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACata                                                                                                                                                                                                                                         <  1:363780/251‑1 (MQ=255)
             aCCGCGATCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTTCCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCGCAACATTACCGGATTATTTACATATGCGGGa                                                                                                                                                                                                                                  <  2:348069/251‑1 (MQ=255)
                gcgAGCTTCGCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACTTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCt                                                                                                                                                                                                                               <  2:454474/251‑1 (MQ=255)
                                     ttttACTGCGTATTGAGCATTACCTGCAACCTGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTATGCGAACGTGACAAcg                                                                                                                                                                                                           <  2:57101/250‑1 (MQ=255)
                                                        ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAgg                                                                                                                                                                                                                         <  2:410262/217‑1 (MQ=255)
                                                        ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAgg                                                                                                                                                                                                                         >  1:410262/1‑217 (MQ=255)
                                                               cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACCGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCAc                                                                                                                                                                                <  2:292822/251‑1 (MQ=255)
                                                                 aCCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTc                                                                                                                                                                              >  2:104393/1‑251 (MQ=255)
                                                                  ccGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTc                                                                                                                                                                              >  2:128966/1‑250 (MQ=255)
                                                                           gtgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAGTTTGTCTGGCAGTCGCCACAACATTACCGGATTAGTTACATATTCGGGATCTGGTAGGCGAACGTGACAACGTTGTATTTTCTGGTTGCGTGCATCCACGCAACCAGAATg                                                                                                                                                                    <  2:569412/251‑1 (MQ=255)
                                                                                      tCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGAt                                                                                                                                                         >  1:473801/1‑251 (MQ=255)
                                                                                               ttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCTAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTagaag                                                                                                                                                  <  1:81768/249‑1 (MQ=255)
                                                                                                      agaCAGCATGTTTTTAGTCGACGCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATGTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTGACGcc                                                                                                                                         <  2:565184/251‑1 (MQ=255)
                                                                                                        aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCCGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGt                                                                                                                                       >  1:99447/1‑251 (MQ=255)
                                                                                                                tttttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggcggc                                                                                                                               >  1:276707/1‑251 (MQ=255)
                                                                                                                   ttAGTCGACTCACACTGCCATCTCGATGGTCTGGCTTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggcggc                                                                                                                               >  2:35057/1‑248 (MQ=255)
                                                                                                                     aGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCagaag                                                                                                                          <  2:276707/251‑1 (MQ=255)
                                                                                                                          aCTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCCTAAGGACGTGGATGCCGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCAACGAt                                                                                                                                                         >  2:525192/1‑215 (MQ=255)
                                                                                                                          aCTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGAt                                                                                                                                                         <  1:525192/215‑1 (MQ=255)
                                                                                                                              acacTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTTGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTa                                                                                                                                                             <  1:230141/207‑1 (MQ=255)
                                                                                                                              acacTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTa                                                                                                                                                             >  2:230141/1‑207 (MQ=255)
                                                                                                                                              ggTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGCAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGAGCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGGGAAAc                                                                                                 >  2:299356/1‑251 (MQ=255)
                                                                                                                                                                    aTAAGGTCGTGGATGTCGTTCTGGCGAAAGCCGTCGCACGCGATGTGAAATTTTGTCTGGCAGGCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGTCAACGTCTTATTTTCTTGTGGCGTGCATCCACTCATCCTGAATGATCCCTACGATGTAGAAGATTTACGCCGTCGGGCGGCAGAAGAGGGTGTTGTAGCGCTGTGTGAAACCGTGCGGGattattatttcccg                                                                           <  2:579544/251‑6 (MQ=255)
                                                                                                                                                                    aTAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACg                                                                           <  1:104393/251‑1 (MQ=255)
                                                                                                                                                                                    cGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTAc                                                           >  1:68712/1‑251 (MQ=255)
                                                                                                                                                                                     gTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACg                                                          >  1:539954/1‑251 (MQ=255)
                                                                                                                                                                                               caGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGCTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGCTGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGGTTATTATTACACGCCGGAAACTAAAGTACGTCAGCAag                                                  >  2:78207/2‑249 (MQ=255)
                                                                                                                                                                                                    gccgcACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTAGTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCtt                                           <  1:325037/251‑1 (MQ=255)
                                                                                                                                                                                                        cACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTcatc                                       >  2:507092/1‑251 (MQ=255)
                                                                                                                                                                                                                           tCTGGCAGTCGCCCCAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCCTCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTcatccagatcggaaaagcgtcgg                    >  2:297431/1‑234 (MQ=255)
                                                                                                                                                                                                                           tCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGTATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTcatcc                                      <  1:297431/233‑1 (MQ=255)
                                                                                                                                                                                                                              ggCAGTCGCCACAACATTACCGGATTATTTACATAGGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGt                  <  2:449048/250‑1 (MQ=255)
                                                                                                                                                                                                                                cAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGtga                >  1:174224/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gCCACACCATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGCACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGCTTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGATGGGAGAAACCGGGCTGGGTTATTATTACACGCCGGAAACTAAAGTACGCCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGtgaactga           >  2:250150/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ccACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggc                                                                                                                                  <  2:194597/130‑1 (MQ=255)
                                                                                                                                                                                                                                      ccACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggc                                                                                                                                  >  1:194597/1‑130 (MQ=255)
                                                                                                                                                                                                                                             aTTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGTGAACTGAACAAGCCgg  >  1:301679/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
GTTACAGGCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGGGCTGGATTATTATTACACGCCGGAAACTAAAGTACGTCAGCAAGAGTCCTTCATCCATCATATCCAGATTGGTCGTGAACTGAACAAGCCGG  >  NC_000913/1156668‑1157155

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