Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 912,638 G→A 6.1% R393C (CGT→TGT)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913912,6380GA 6.1% 211.9 / 3.0 66R393C (CGT→TGT) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

AGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATCACCGTCC  >  NC_000913/912493‑912783
                                                                                                                                                 |                                                                                                                                                 
aGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGc                                                                                                                                                 >  4:89002/1‑148 (MQ=255)
     gTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCAc                                                                                                                                             <  3:112402/147‑1 (MQ=255)
      ttAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCAcc                                                                                                                                            <  4:4231/147‑1 (MQ=255)
         aaGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGa                                                                                                                                         >  4:219258/1‑147 (MQ=255)
             gATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCTCCGCGTTCGCCGTCACAGCCACCAAGCAGGAAGATATGAAGCAGTTTTTCACGGCTCACCAGatcaat                                                                                                                                    >  3:53080/1‑148 (MQ=255)
                 tttACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAg                                                                                                                                 <  4:209409/147‑1 (MQ=255)
                 tttACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGc                                                                                                                                <  3:219258/148‑1 (MQ=255)
                     ccACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTa                                                                                                                             <  3:89002/147‑1 (MQ=255)
                      cacaGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATc                                                                                                                           <  3:64954/148‑1 (MQ=255)
                      cacaGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATc                                                                                                                           <  4:159462/148‑1 (MQ=255)
                       acaGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTAt                                                                                                                            >  1:155544/1‑146 (MQ=255)
                         aGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATcagc                                                                                                                        <  1:65253/148‑1 (MQ=255)
                                  tCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAag                                                                                                               >  2:194345/1‑148 (MQ=255)
                                  tCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAag                                                                                                               >  6:224353/1‑148 (MQ=255)
                                     ggATCAGGCAGTCATCCGGCACGCTGTTGGCGAAATCGGTGAAGTAGTGTCGCTCGGCGCGTGCGCCGTCACAGCCACATAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCTGCGTATCAGCAGCGCCAagcag                                                                                                            <  3:32480/148‑1 (MQ=255)
                                     ggATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAagcag                                                                                                            <  7:48899/148‑1 (MQ=255)
                                        tCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGt                                                                                                         >  8:292899/1‑148 (MQ=255)
                                         cAGGCAGTCATCCGGCACGCTGGTGGCGAAATAGGTGAAGTAGTCGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGt                                                                                                         <  5:67912/147‑1 (MQ=255)
                                          aGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCt                                                                                                       <  1:194345/148‑1 (MQ=255)
                                            gCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTgg                                                                                                     >  2:62459/1‑148 (MQ=255)
                                             cAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTAt                                                                                                                            <  2:222952/124‑1 (MQ=255)
                                             cAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTAt                                                                                                                            >  1:222952/1‑124 (MQ=255)
                                             cAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGc                                                                                                    >  8:70929/1‑148 (MQ=255)
                                                   tCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAg                                                                                                                                          <  2:244346/104‑1 (MQ=255)
                                                   tCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAg                                                                                                                                          >  1:244346/1‑104 (MQ=255)
                                                   tCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCa                                                                                               >  7:178990/1‑147 (MQ=255)
                                                    ccGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCa                                                                                                                                           <  4:25333/102‑1 (MQ=255)
                                                    ccGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCa                                                                                                                                           >  3:25333/1‑102 (MQ=255)
                                                    ccGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaa                                                                                              >  1:231698/1‑147 (MQ=255)
                                                     cGGCACGCTGGGGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaa                                                                                             <  5:182680/147‑1 (MQ=255)
                                                        cACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAAccc                                                                                         <  5:64128/148‑1 (MQ=255)
                                                           gCTGGTGGCGACATCGGGGAAGGAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACg                                                                                      <  2:151378/148‑1 (MQ=255)
                                                                 ggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaa                                                                                             >  8:178760/1‑135 (MQ=255)
                                                                 ggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaa                                                                                             <  7:178760/135‑1 (MQ=255)
                                                                 ggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCAAAGCAGCGTCTGGCGACCAAAACCCACGGTGAt                                                                                 >  8:328310/1‑147 (MQ=255)
                                                                    gAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAgg                                                                             <  1:62459/148‑1 (MQ=255)
                                                                      aaTAGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACACCCACCAAGCAGGAAGATATGACGCAGTTTTTAACGGCTCACCAGCTAAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCCACCAAAACCCACGGTGATAAGGTg                                                                           >  7:134644/1‑148 (MQ=255)
                                                                          ggTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCGCAGCCCCCAAGCAGGAAGATATGCCGCAGTTTTTCACGGCTCACGAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCgg                                                                        >  4:50395/1‑147 (MQ=255)
                                                                            tGAAGTAGTCGCGCGCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAAt                                                                     <  2:229987/148‑1 (MQ=255)
                                                                               agtagtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTc                                                                  <  8:178990/148‑1 (MQ=255)
                                                                               agtagtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTc                                                                  <  7:328310/148‑1 (MQ=255)
                                                                                  agtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaa                                                                                             <  8:128507/118‑1 (MQ=255)
                                                                                  agtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaa                                                                                             >  7:128507/1‑118 (MQ=255)
                                                                                    tGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGCCGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTg                                                              >  1:181583/1‑147 (MQ=255)
                                                                                         cTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTAAggg                                                        <  4:53080/148‑1 (MQ=255)
                                                                                            gccgcGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGa                                                       <  6:240441/146‑1 (MQ=255)
                                                                                                cgTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTg                                                 <  7:292899/148‑1 (MQ=255)
                                                                                                cgTGCCCAGTCACAGCCACCAAGCAGGAAGATATGAAGCAGTTTTTCCCTGCTCACCAGATCAAACAGCGTATCAGCAGCGCCAAGAAGCGACTGGCGACCAACACCCACCGACAAAAGGTTCGGAATTTCGCAGTACGGGAAGCCTg                                                 >  8:43460/1‑148 (MQ=255)
                                                                                                    cgcCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCa                                              <  2:155544/147‑1 (MQ=255)
                                                                                                      ccGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGGATCCGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGGGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCt                                           >  1:53326/1‑148 (MQ=255)
                                                                                                         tCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGtt                                        >  5:180320/1‑148 (MQ=255)
                                                                                                                caccaAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGGGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTg                                 <  7:43460/148‑1 (MQ=255)
                                                                                                                     aGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTg                             <  1:6481/147‑1 (MQ=255)
                                                                                                                     aGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTg                             >  2:6481/1‑147 (MQ=255)
                                                                                                                       cAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATa                          <  5:224353/148‑1 (MQ=255)
                                                                                                                               tatGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCa                    >  4:243371/1‑146 (MQ=255)
                                                                                                                                    cGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAat             <  2:181583/148‑1 (MQ=255)
                                                                                                                                     gCAGTTTTTCACAGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGATTTTCGCTGTACGGGAAGCCTGCCATCt                                           <  5:251225/117‑1 (MQ=255)
                                                                                                                                      cAGTTTTTCACAGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCt                                           >  6:255510/1‑116 (MQ=255)
                                                                                                                                      cAGTTTTTCACAGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCt                                           <  5:255510/116‑1 (MQ=255)
                                                                                                                                      cAGTTTTTCACAGCTCACCAGATCAAACAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCAACGGTGATAAGGTGCGGAATTTCGCAGTAAGGGAAGCATGCCATCTg                                          >  6:251225/1‑117 (MQ=255)
                                                                                                                                       aGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAatcat          >  5:83247/1‑148 (MQ=255)
                                                                                                                                          ttttCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATCAc       >  7:112391/1‑148 (MQ=255)
                                                                                                                                               aCGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATCACCGTc   <  2:231698/147‑1 (MQ=255)
                                                                                                                                               aCGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATAACCGTcc  >  2:86286/1‑148 (MQ=255)
                                                                                                                                                cGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATCACCGTc   >  3:220886/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTGCCATCTGTTGCGCCTGGGTGATAACCGCAGAGAAATCATCACCGTCC  >  NC_000913/912493‑912783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑