Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,288,385 G→A V208I (GTT→ATT)  nrfA → cytochrome c552 nitrite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,288,3850GA87.2% 89.4 / 6.1 40V208I (GTT→ATT) nrfAcytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (16/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-02

GTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCA  >  NC_000913/4288158‑4288627
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
gTCTGACCCAGAAAGACGGGGAAGAGGGCTCCTCCCACGGTAAAGGGGCGCGCGGCGGCCCGGAAATCGTCAACAATTTAGGTTGGGCCGTTTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAg                                                                                                                                                                                                                             <  1:103291/251‑1 (MQ=255)
         aGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTc                                                                                                                                                                                                                    >  1:485965/1‑251 (MQ=255)
             agagggCGAAGATGTCTACTACAGCGGTAAAGGGGCGCGCGGCGGTCCGAAAACCGTAAACAACTGATGTTGGCCGGATTGCCTAAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAAGCGATGATTTGCGGTCAGTGCTATGTGGAGTATTACTTCGACg                                                                                                                                                                                                                <  1:378242/247‑1 (MQ=255)
                        aTGGCTACTTCACCGGAAAATAGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTCGGTTGTGCCGATTGCCCTAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACAGGAGTTAACCCTTTCCCGTCCGTATGCGGCCCGCGCGATGGAAACCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGTTGATTTGCGGTAAGTGACATGAGGAGCATAACTTCGACGGAGAAAACaa                                                                                                                                                                                                      >  2:252594/1‑250 (MQ=255)
                             tACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGt                                                                                                                                                                                                >  1:188726/1‑251 (MQ=255)
                                     gggTAAATGGGCGGGGGGCGGTCCGGAAATCGTCAACAACTTAGGTTGGGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTcc                                                                                                                                                                                        <  1:449221/250‑1 (MQ=255)
                                           gtGGGCGCGCGGCGCTCCGGAAATCGACAACAACTTATGTGGTGCCGATTCCCATACCACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTggg                                                                                                                                                                                  <  1:456837/250‑1 (MQ=255)
                                                             ggAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCg                                                                                                                                                                >  1:30059/1‑251 (MQ=255)
                                                                aaccgaaaacaatTTAGTTTGTGCGGATGGCCATACCACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGaaa                                                                                                                                                             <  1:199503/244‑1 (MQ=255)
                                                                              tAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTAt                                                                                                                                               >  1:364194/1‑251 (MQ=255)
                                                                               gggTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTAtt                                                                                                                                              >  1:110637/2‑251 (MQ=255)
                                                                               aGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAACGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGCAAATATGGAGCAGTAtt                                                                                                                                              >  1:110625/1‑251 (MQ=255)
                                                                                  ggtggcccgatggccaaacccccGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTAc                                                                                                                                            <  2:30059/230‑1 (MQ=255)
                                                                                    gtgCCGATTCCCAAAACCCCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGAc                                                                                                                                         <  2:91201/251‑1 (MQ=255)
                                                                                       ccGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGATTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGCGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGGTGACGGCATGCAAGTCgcaa                                                                                                                                                             >  2:120622/1‑225 (MQ=255)
                                                                                                             ctcCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACt                                                                                                                 >  1:248141/1‑251 (MQ=255)
                                                                                                               ccAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACACAATTGCCTTATCTGACTGGACTa                                                                                                                >  1:309092/1‑250 (MQ=255)
                                                                                                               ccAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGAc                                                                                                                  <  2:364194/248‑1 (MQ=255)
                                                                                                               ccAGAGTGCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGACCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGAc                                                                                                                  <  2:130235/248‑1 (MQ=255)
                                                                                                                 agagTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTTCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACt                                                                                                             <  1:115403/251‑1 (MQ=255)
                                                                                                                   agTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAACGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGATCAACTcc                                                                                                           >  2:441291/1‑251 (MQ=255)
                                                                                                                     ttCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCt                                                                                                         >  1:83668/1‑251 (MQ=255)
                                                                                                                       cGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTg                                                                                                        <  2:64584/250‑1 (MQ=255)
                                                                                                                       cGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCATTCTCTGACTGGACTAACTCCCTGt                                                                                                       >  2:240124/1‑251 (MQ=255)
                                                                                                                         ccAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCg                                                                                                     <  1:122808/251‑1 (MQ=255)
                                                                                                                          cAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCg                                                                                                     >  1:278562/1‑250 (MQ=255)
                                                                                                                                       gAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCGTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCt                                                                                       <  2:485965/251‑1 (MQ=255)
                                                                                                                                                             aTGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                                                 <  2:371490/251‑1 (MQ=255)
                                                                                                                                                                 ggCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGaaa                                                             <  1:320043/251‑1 (MQ=255)
                                                                                                                                                                 ggCTCGCGCGATGGAAGCCATTGGTAAAACTTTTGAGAAAGCAGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCAAGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAAACGAAAATATGGCGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTCACGAAAAAGCAAATGCTGAAAGCGCAGCACCCGGAATATGaa                                                              >  1:278225/1‑250 (MQ=255)
                                                                                                                                                                  gCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGAc                                                                                                                                                                                                                 <  2:92212/102‑1 (MQ=255)
                                                                                                                                                                  gCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGAc                                                                                                                                                                                                                 >  1:92212/1‑102 (MQ=255)
                                                                                                                                                                             ggAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGCTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCgg                                                 >  1:204050/1‑251 (MQ=255)
                                                                                                                                                                              gAAGCCATTGGTAAACCTTTTGAGAAATCCGGACGGTTCTACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTAGTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATGACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCggg                                                <  2:345042/251‑1 (MQ=255)
                                                                                                                                                                                 gCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCACAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAAAATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCTAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAt                                             >  2:381738/1‑251 (MQ=255)
                                                                                                                                                                                       ggTAAACCTTTGGAGAAAGCCGGAAGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGTACAGCGGGCATTCACg                                        <  2:29093/250‑1 (MQ=255)
                                                                                                                                                                                                             gACGTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTGAAAACAACGTGACCTGTATc                 >  1:96749/1‑251 (MQ=255)
                                                                                                                                                                                                                gTTTCGACCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGAc              >  1:187910/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCCAGCAATCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCat        <  2:204050/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaTCGATGATTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCa  <  1:324646/251‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
GTCTGATCCAGAAAGACGGCGAAGATGGCTACTTCCACGGTAAATGGGCGCGCGGCGGTCCGGAAATCGTCAACAACTTAGGTTGTGCCGATTGCCATAACACCGCCTCTCCAGAGTTCGCCAAAGGCAAACCGGAGTTAACCCTTTCCCGTCCGTATGCGGCTCGCGCGATGGAAGCCATTGGTAAACCTTTTGAGAAAGCCGGACGTTTCGACCAGCAATCGATGGTTTGCGGTCAGTGCCATGTGGAGTATTACTTCGACGGCAAAAACAAAGCGGTTAAATTCCCGTGGGATGACGGCATGAAAGTCGAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCTGTATCGACTGCCATATGCCA  >  NC_000913/4288158‑4288627

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