Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,781 2 bp→AT coding (239‑240/597 nt) nrfG → putative formate‑dependent nitrite reductase complex subunit NrfG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,7810GA96.2% 167.5 / NA 54G80D (GGC→GAT) nrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/24);  minor base C (2/0);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,293,7820CT100.0% 178.5 / NA 54G80D (GGC→GAT) nrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/25);  total (29/25)

GCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTA  >  NC_000913/4293540‑4294021
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                 
gCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTa                                                                                                                                                                                                                                           <  1:301574/251‑1 (MQ=255)
   aTGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCt                                                                                                                                                                                                                                        <  1:192866/251‑1 (MQ=255)
     gaacctgCCCAAATTACCGGTGACAATGCTGGAAACCCTCACGATATTGATGGTAATCCTCTGTGTCGCCAGATATCTGTAAAGCCCAAAATGACAGGCGGCACGTGCGGCGTATAAGCGCCAGCGCGATCCACTCCATCAATTTGCCAGCCCGCCAACACCGGAAGCGCAGATTCCCGCATTGCAGGTTAAAATCCGCGCTAATCCACAAAACAGCGAACAGGGGGCGTTACTGGCTGAGTACTATcgat                                                                                                                                                                                                                                      >  2:24672/7‑248 (MQ=255)
        aCAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGc                                                                                                                                                                                                                                   >  1:309042/1‑251 (MQ=255)
               aaaaTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTAGGCGTTACTGGATGAGTACTATCTGTGGCAAAACg                                                                                                                                                                                                                             >  2:194698/1‑250 (MQ=255)
                    aCCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTAc                                                                                                                                                                                                                        <  1:26163/250‑1 (MQ=255)
                              aTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGtt                                                                                                                                                                                                             >  1:32553/1‑251 (MQ=255)
                              aTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAACCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCGTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTCATGGATGAGTACTATCTGTGGAAAAACGATTACAGCCATTctgt                                                                                                                                                                                                             >  2:159241/1‑248 (MQ=255)
                                         acTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTCTCTGTTCAGCCCAACATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCCGTTTGCCAGCCAGCAAACCCCGGAAGCGCCGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTAcc                                                                                                                                                                                                  >  1:190279/1‑251 (MQ=255)
                                         acTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTAcc                                                                                                                                                                                                  >  1:398364/1‑251 (MQ=255)
                                                     gATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAAGCGTTg                                                                                                                                                                                       <  1:114352/250‑1 (MQ=255)
                                                             ttCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGCGCGGAGTATCAGCGTCAGAGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGAAAAACGAATACAGCAATTCGTTGCTGGCGGACCGTCAGGCGTTGCAACTaagg                                                                                                                                                                              >  2:217785/1‑248 (MQ=255)
                                                              tctcTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCCACTACgtg                                                                                                                                                                             >  2:79327/1‑251 (MQ=255)
                                                              tctcTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtg                                                                                                                                                                             >  1:147572/1‑251 (MQ=255)
                                                                   gtgtCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggt                                                                                                                                                                           <  1:290599/248‑1 (MQ=255)
                                                                   gtggcgggaGTTATCTGTTAAGCCCAAAATGGCAGGGGGTACGTGCGGAGTAGCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAAGCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggtg                                                                                                                                                                          <  2:408201/241‑1 (MQ=255)
                                                                           aGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGa                                                                                                                                                               <  1:378118/251‑1 (MQ=255)
                                                                              tATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACt                                                                                                                                                            >  1:92308/1‑251 (MQ=255)
                                                                                  tGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCCCAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTAt                                                                                                                                                                                                                                          >  1:385561/1‑170 (MQ=255)
                                                                                  tGTTAAGCCCAAAAGTGCAGGCGGTGCGTGCGGAGTATCAGCGTCAGCGCGATCCGCGACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTAt                                                                                                                                                                                                                                          <  2:385561/170‑1 (MQ=255)
                                                                                   gTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTa                                                                                                                                                         <  2:234393/249‑1 (MQ=255)
                                                                                      aaGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTGGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATgcg                                                                                                                                                     <  1:264768/250‑1 (MQ=255)
                                                                                      aaGCCCAAAATGGCAGGCGGTACGGGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTAGCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACGACGTGGTGAGAACGC‑TGAACTTTATgcgg                                                                                                                                                    <  2:374475/251‑1 (MQ=255)
                                                                                              aaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGAAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACCGCGAACCGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAAAGC‑TGAACTGTATGCGGCGCTGGCg                                                                                                                                            >  2:257900/1‑251 (MQ=255)
                                                                                                         gTACGGGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGAAAGCAAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTACTCCACAAAACAGCGAACAGTGGGCGTTCCTGGCTGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTAGCAACTACGTGGGGAGACCGC‑TGAACTGTATGCGGCGATGTCGACGGTGtttt                                                                                                                                  >  2:344552/1‑250 (MQ=255)
                                                                                                           aCGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTAGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                               <  1:194698/251‑1 (MQ=255)
                                                                                                            cGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCCGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAACTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                               >  1:404850/1‑250 (MQ=255)
                                                                                                                cGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTAcc                                                                                                                            >  2:150773/1‑249 (MQ=255)
                                                                                                                    gTATCAGCGTCAGCGCGATCCGCTACATCAGTTTACCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TTAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgcca                                                                                                                      <  1:257829/251‑1 (MQ=255)
                                                                                                                        cagcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TTAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagcc                                                                                                                   <  1:362526/250‑1 (MQ=255)
                                                                                                                                gcgcGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATg                                                                                                           >  1:406391/1‑250 (MQ=255)
                                                                                                                                      tCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGccc                                                                                                   >  1:316437/1‑251 (MQ=255)
                                                                                                                                             cATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTcg                                                                                            >  1:169006/1‑251 (MQ=255)
                                                                                                                                             cATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGCGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTcg                                                                                            >  1:48208/1‑251 (MQ=255)
                                                                                                                                                     tGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGACGTTACTGGATGAGTACTATCTATGGAAAAACGATTACAGCAATTCGATGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATg                                                                                     >  1:359566/1‑250 (MQ=255)
                                                                                                                                                        cagccagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATGCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCg                                                                                 <  1:287837/251‑1 (MQ=255)
                                                                                                                                                         agccagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGa                                                                                <  2:309042/251‑1 (MQ=255)
                                                                                                                                                              gcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAg                                                                           <  2:349952/251‑1 (MQ=255)
                                                                                                                                                                       ggAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCCACTACGTGGTGAGAACGC‑TGAACTGTACGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCgcgc                                                                  >  1:306265/1‑251 (MQ=255)
                                                                                                                                                                                         aTTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATc                                                 >  1:381023/1‑250 (MQ=255)
                                                                                                                                                                                               ggataacctCCGCGCTAATCCACACAACAGCGAACAGTGGGCGCTACTGGATGAGCACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGAGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGccc                                          >  1:160462/9‑251 (MQ=255)
                                                                                                                                                                                                 aTAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGGTTACAGCGATTCGTTGCTGGCGTACCGTCAGGCGGTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGGGCTTTATTACCGGGTCAGCCAGGAT‑TGTACCGCCCAGAGTCGCGCAAGGATCGATAAGTGCCTTGCGCTGGGCAGTAAAGAAAGCAGCGCCCtt                                        <  2:160462/251‑2 (MQ=255)
                                                                                                                                                                                                    aaaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtg                                     <  1:8622/251‑1 (MQ=255)
                                                                                                                                                                                                       tCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgctg                                  >  1:24247/1‑251 (MQ=255)
                                                                                                                                                                                                        ccGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgctgc                                 >  1:409246/1‑251 (MQ=255)
                                                                                                                                                                                                             cTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCt                            >  1:74758/1‑251 (MQ=255)
                                                                                                                                                                                                                aTCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTAGCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑GGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTTATGCTGCTGTCTTcc                         <  2:248164/251‑1 (MQ=255)
                                                                                                                                                                                                                    acaAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATg                     <  1:235533/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaCAGGGAACTGTGGGCTTTACTGGTTGAGTCATATCTGTGGCAAAACGATTACAGCAATTCGGTGCTGGGGTACCGTCAGGCGTTGCACCCTCGTGGTGAGAACGCCTGA‑CTCTATGCGGCGCTGGCGACGGTGCTTGATTACCAGGCCAGCCAGTATATG‑ACCGCCCAGTCTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCGCCGCCCTGATGCTGCGGGCTTCCGATGCg                   <  2:80690/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGttt                >  1:98884/1‑251 (MQ=255)
                                                                                                                                                                                                                                aCAGTGGTCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAgg         <  2:406391/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTa  <  1:79327/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCta                                                                      >  2:374983/1‑180 (MQ=255)
                                                                                                                                                                                                                                         gTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTc                                                                      <  1:374983/181‑1 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                 
GCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGC‑TGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATG‑ACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTA  >  NC_000913/4293540‑4294021

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