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,7810GA97.7% 137.7 / ‑4.8 44G80D (GGC→GAT) nrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/23);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
*NC_0009134,293,7820CT97.7% 148.7 / ‑4.2 44G80D (GGC→GATnrfGputative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/23);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGC  >  NC_000913/4293542‑4294011
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                     
gATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATc                                                                                                                                                                                                                             <  1:312352/251‑1 (MQ=255)
            cAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTGCATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACg                                                                                                                                                                                                                 <  2:280831/251‑1 (MQ=255)
              aaaTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTCAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCCTTGCAGGATAAAATCCGCGCTAATCCACAAAACCGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGAt                                                                                                                                                                                                               >  2:326492/1‑251 (MQ=255)
               aaTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGAtt                                                                                                                                                                                                              <  2:89813/251‑1 (MQ=255)
                aTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGAtt                                                                                                                                                                                                              >  2:256478/1‑250 (MQ=255)
                            aagCTGACTACGCTCAAGATATTGATGGTATTTCTCTGTGTCGGCAGTTACCTGTTAAGCCCAAAAGTGCTGGCGGTACGTGTGCAGTATCAGCGTCAGCGCGAGCCGCTTCATCAGTTTTCGATCCAGCAATCTCCGGAAGCGCAGCTTGATGCACTGCAGGATAAAATTCGCGCTATTCCACCAAACATCGAAGAGTTTGCGTTACTGCATGGGTATTATCCATGGCAAAGCGATTACAGCAATTCGtt                                                                                                                                                                                                 <  2:261136/249‑1 (MQ=255)
                                  acaacaCTAACGCTATTGATAGTAATTCTCTGTGTCCGCAGTTATCTGTTAACCCCAACATGGCAGGCGGTAAGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGAAGATTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATATGTGTCAAAACGATTCCAGCAATTCGTTGCTGGc                                                                                                                                                                                           >  2:387626/1‑251 (MQ=255)
                                       acTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTAcc                                                                                                                                                                                      <  1:288545/251‑1 (MQ=255)
                                          cACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGAAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCg                                                                                                                                                                                     <  2:282912/249‑1 (MQ=255)
                                           aCGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCa                                                                                                                                                                                  <  1:345840/251‑1 (MQ=255)
                                              atatTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGc                                                                                                                                                                               >  1:376337/1‑251 (MQ=255)
                                              atagtGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGc                                                                                                                                                                               <  1:317010/247‑1 (MQ=255)
                                                atTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCAGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTTGCGTACCGTCAGGCGt                                                                                                                                                                             >  2:356390/1‑251 (MQ=255)
                                                       gTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGATTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACt                                                                                                                                                                      <  1:117961/251‑1 (MQ=255)
                                                                 gtgtCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggtg                                                                                                                                                              >  2:129950/1‑249 (MQ=255)
                                                                     cGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACg                                                                                                                                                        >  1:73968/1‑251 (MQ=255)
                                                                         aGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTg                                                                                                                                                     >  2:51750/1‑250 (MQ=255)
                                                                                 gTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTAt                                                                                                                                             <  1:194482/250‑1 (MQ=255)
                                                                                   tAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATgcg                                                                                                                                          >  2:380147/1‑251 (MQ=255)
                                                                                            aaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCg                                                                                                                                 <  2:87006/251‑1 (MQ=255)
                                                                                            aaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCg                                                                                                                                 >  1:253212/1‑251 (MQ=255)
                                                                                                             gCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTAcc                                                                                                                 >  2:110497/1‑250 (MQ=255)
                                                                                                               ggAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAgg                                                                                                              <  2:27267/251‑1 (MQ=255)
                                                                                                               ggAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAgg                                                                                                              <  1:110497/251‑1 (MQ=255)
                                                                                                                 aGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgcc                                                                                                            >  2:130813/1‑251 (MQ=255)
                                                                                                                       agcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGc                                                                                                                                  <  1:241851/223‑1 (MQ=255)
                                                                                                                       agcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGc                                                                                                                                  >  2:241851/1‑223 (MQ=255)
                                                                                                                       agcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGATTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagccac                                                                                                      >  2:129167/1‑250 (MQ=255)
                                                                                                                         cgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCa                                                                                                    <  1:193970/251‑1 (MQ=255)
                                                                                                                                    tCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCaa                                                                                                                                                                                                       >  2:114827/1‑141 (MQ=255)
                                                                                                                                    tCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCaa                                                                                                                                                                                                       <  1:114827/141‑1 (MQ=255)
                                                                                                                                                 tttGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCatt                                                                            >  2:163273/1‑249 (MQ=255)
                                                                                                                                                                                 tCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAAt                                             <  1:212849/250‑1 (MQ=255)
                                                                                                                                                                                   aGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGaa                                          <  1:67666/251‑1 (MQ=255)
                                                                                                                                                                                   aGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGaa                                          <  1:257878/251‑1 (MQ=255)
                                                                                                                                                                                          gCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCAcc                                    <  2:200599/250‑1 (MQ=255)
                                                                                                                                                                                                   aaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgc                          >  1:97514/1‑251 (MQ=255)
                                                                                                                                                                                                             aaTCCACACAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTc                >  2:39468/1‑251 (MQ=255)
                                                                                                                                                                                                                ccACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCg              <  1:93742/250‑1 (MQ=255)
                                                                                                                                                                                                                 cacaAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTGCCAGGCCAGCCAGCATATGACCGCCCAGACTCGCTCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGAt            <  2:210775/251‑1 (MQ=255)
                                                                                                                                                                                                                  acaAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATg           >  1:413968/1‑251 (MQ=255)
                                                                                                                                                                                                                        cAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                    <  2:296984/140‑1 (MQ=255)
                                                                                                                                                                                                                        cAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                    >  1:296984/1‑140 (MQ=255)
                                                                                                                                                                                                                            gAACAGCGGGCGTTACTGGATGAGTACTATCGGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCCACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATCTGACCGCCCAGACTCGCGCAATGAACGACAACGCCCTCGCGCTGGACCGGACTGAAATAACCGCCCTGATGCTGCCGGCTTCCGTTGCGTTTATGc  >  2:397275/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                     
GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGC  >  NC_000913/4293542‑4294011

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