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,7810GA90.9% 114.4 / ‑5.9 45G80D (GGC→GAT) nrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (17/23);  minor base C (3/1);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,293,7820CT100.0% 136.7 / NA 45G80D (GGC→GATnrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/24);  total (21/24)

CCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAAC  >  NC_000913/4293553‑4294019
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
ccAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAAc                                                                                                                                                                                                                          <  1:301385/251‑1 (MQ=255)
     aTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCTGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGAtt                                                                                                                                                                                                                      >  1:199951/1‑250 (MQ=255)
          ggTGAAAATGCTGACAACACTCACGATGTTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCACAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGTGAGTACTTTCTGTGGCAAAACGATTACAga                                                                                                                                                                                                                 >  2:63471/1‑249 (MQ=255)
             gAAAATGCTGACAGAACTCCCGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGTAGTATAAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAAtt                                                                                                                                                                                                             <  1:194351/251‑1 (MQ=255)
              aaaaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTc                                                                                                                                                                                                            <  1:231882/251‑1 (MQ=255)
               aaaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCg                                                                                                                                                                                                           <  1:26922/251‑1 (MQ=255)
                         aacaCTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAACCGCTTACAGCAATTCGTTGCTGGCGt                                                                                                                                                                                                 >  1:413560/1‑251 (MQ=255)
                          acacTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTa                                                                                                                                                                                                <  1:147082/251‑1 (MQ=255)
                                  gATAGTGATGGTATTTCCCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGGTCGGAGTATCATCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAgg                                                                                                                                                                                        <  2:253147/251‑1 (MQ=255)
                                    tatTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGCTGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCg                                                                                                                                                                                      >  2:88769/1‑251 (MQ=255)
                                            gTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAACATGTCATTCGTTACTTTCGGAGTATCATCGTCAGCGCGCTCCGCTCCCTCCGTTTGCCAGCCCGCAAACCCCGGAAGCGCAGCTTCAGGCATTTCAGGATCAAATCCGCGCTAATCCCCAAAACATCGAACATTGGGCTTTACTGGATTGTTACTATCTTTGTCAAACCACTTACCGCCATTCGTTGCTGGCGTACCTTCCGGCGTTGCacct                                                                                                                                                                              >  2:127477/1‑248 (MQ=255)
                                                    ctGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGTAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggtg                                                                                                                                                                      >  1:384091/1‑251 (MQ=255)
                                                      gcgtcgtccgTTATCTGTTAAGCACAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCTCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCTTTACTGGCTGAGTCCTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggtg                                                                                                                                                                      <  1:251343/240‑1 (MQ=255)
                                                           ggCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGc                                                                                                                                                               >  2:300894/1‑251 (MQ=255)
                                                             cAGTTATCTTTTAAGTCCAAAAGGGCAGGGGGTACGTGCGGAGTATCATCGTCAGCGCGATCCGTTCCAGCATTTTGCCATCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCATGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGGGACTGGATGAGGACTAGCTGTGGCAAAACGAGTACAGCAATTCGGTGCTGGCGTACCCTCAGGCGTTTCAACTACGTGGTGAGAACGCTg                                                                                                                                                             <  2:400613/251‑1 (MQ=255)
                                                                         aaGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTTGGCGTTACTGGATGAGTACTATCTGTTGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATgcg                                                                                                                                                  <  1:264019/250‑1 (MQ=255)
                                                                         aaGCCCAAAATGGCAGGCGATACGTGCGGAGTATCAGCGTAAGCGCGATCCGCTACATCAGTTTGCCAGCCCGCAAACCCCGGAAGAGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCACTTCGTAGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGGGAACGCTGTACTGTATgcgg                                                                                                                                                 >  2:42904/1‑251 (MQ=255)
                                                                                aaaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGc                                                                                                                                          <  1:220880/251‑1 (MQ=255)
                                                                                    ggCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAACCAGCGAACAGTGGGCGTTAATGGCTGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACg                                                                                                                                      >  1:212679/1‑251 (MQ=255)
                                                                                    ggCAGGCGGAACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAAAAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGCTTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGGGCTGGCGACg                                                                                                                                      >  2:389191/1‑251 (MQ=255)
                                                                                     gCAGGAGGTACGGGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAAGCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCATCTACGTGGTGAGAACGCTGAACTGTATTCGGCGCTTGCGACgg                                                                                                                                     <  2:414860/251‑1 (MQ=255)
                                                                                       aGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAAACCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGCTGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGAGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACgagg                                                                                                                                   >  2:221350/1‑248 (MQ=255)
                                                                                            gTACGTGCTGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTtta                                                                                                                              <  2:199951/251‑1 (MQ=255)
                                                                                                gTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCAAGCACACACCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAACACAGCGAACAGAGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTAc                                                                                                                          >  1:233075/1‑251 (MQ=255)
                                                                                                   cGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAg                                                                                                                       <  1:325340/251‑1 (MQ=255)
                                                                                                      ggTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTGCTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgcc                                                                                                                    <  1:174705/250‑1 (MQ=255)
                                                                                                         aTCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAAACCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagc                                                                                                                 <  2:332768/251‑1 (MQ=255)
                                                                                                                       gATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAAAGCTGAAATGTATGCGGCGCTGGAGACGGTGCTTTATTACCAGGACAGCCAGCATATGACCGc                                                                                                   >  2:180274/1‑251 (MQ=255)
                                                                                                                        aTCCGCAACATCAGTTTCCCAGCCAGCAAACCCCGGCACCGAAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAGTCCAAAAAACAGCTCACAAAGAGCGATAATGGATGGATACTATCTATGGCAAAACGATTAAAGCAATTCGTTGCGAGCGTACCGTCAGAGGTTGCAACTACGTAGAGAGAAAGCTGAGCTGTATGAGGCGCTGACAGCGGCGCTTTAGTACCAagccagccgcaagataaccccc                                                                                                  >  2:281934/1‑237 (MQ=255)
                                                                                                                               aCATCAGGTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGAc                                                                                             <  1:63471/249‑1 (MQ=255)
                                                                                                                                     gTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAAGCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCaa                                                                                     <  2:285814/251‑1 (MQ=255)
                                                                                                                                                   aaaCCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGcc                                                                       >  1:330568/1‑251 (MQ=255)
                                                                                                                                                     cccccGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCt                                                                     <  2:384091/250‑1 (MQ=255)
                                                                                                                                                                               gCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCAcc                                            <  1:180274/250‑1 (MQ=255)
                                                                                                                                                                                   gATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATAAGTACTATCTGTGGCAACACGTTTACAGCAATTCGTTGCTTGCGTACCGTCAGTCGTTGCAACTACGTGGTGAGAACTCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGCAATCACCGCCCt                                       >  1:276615/1‑251 (MQ=255)
                                                                                                                                                                                    aTAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGAAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTg                                      >  1:329397/1‑251 (MQ=255)
                                                                                                                                                                                       aaaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGAGGCAAAACGATTACAGCAATTCGTTGCTGGAGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACAGCCCAGACTCGCGCAATGATCGCCAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCGGAtg                                   >  2:22500/1‑251 (MQ=255)
                                                                                                                                                                                        aatcctctcTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGTTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGTACAGTAATGAAATCACCGCCCTGAtgc                                  <  1:149828/243‑1 (MQ=255)
                                                                                                                                                                                         aTCCGCACTAATCCACAACACCGAGACCAGTGGGCGTTACTGGATGCGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTAGTGAGAACGCTGAACTGTATGCGGCGCTGGCGCCGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCCATGATCGACAAAGCCCTCGCGCTGGACCGTAATGAACTCCCCGCCATGAtgct                                 >  2:56486/1‑251 (MQ=255)
                                                                                                                                                                                                  aaTCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAACACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTTGCTTc                        >  2:389140/1‑251 (MQ=255)
                                                                                                                                                                                                          aaaCAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGt                <  1:10844/251‑1 (MQ=255)
                                                                                                                                                                                                                 gAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGc          >  1:370831/1‑250 (MQ=255)
                                                                                                                                                                                                                  aaCAGTGGGCGTTACGGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCTCGCTGGACAGTAATGAAAGCACCGCCCGGATGCTGCTGGCGTCCGAGGCGTGTATGCAg        <  2:255022/251‑1 (MQ=255)
                                                                                                                                                                                                                        gggCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAAc  <  1:375960/251‑1 (MQ=255)
                                                                                                                                                                                                                        gggCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAAc  <  1:300894/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
CCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAAC  >  NC_000913/4293553‑4294019

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