Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 913,349 G→A L156L (CTG→TTG)  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_000913913,3490GA100.0% 255.7 / NA 78L156L (CTG→TTG) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/35);  total (43/35)

CGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGAC  >  NC_000913/913218‑913484
                                                                                                                                   |                                                                                                                                       
cgcgTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCAcagc                                                                                                                                  >  2:1734108/1‑139 (MQ=255)
    ttCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGAc                                                                                                                              >  2:679298/1‑139 (MQ=255)
                                            tATGGTACTGGGGATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                                   >  2:684046/1‑126 (MQ=255)
                                            tATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCTGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                                   <  1:684046/126‑1 (MQ=255)
                                                   cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGcc                                                                                                                       <  1:1410680/99‑1 (MQ=255)
                                                   cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGcc                                                                                                                       >  2:1410680/1‑99 (MQ=255)
                                                   cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAggg                                                                               >  1:815302/1‑139 (MQ=255)
                                                    tGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAgggg                                                                              >  2:1605601/1‑139 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             <  1:685998/90‑1 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             >  2:685998/1‑90 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATg                                                                                                                >  2:518412/1‑103 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATg                                                                                                                <  1:518412/103‑1 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                                                <  2:723608/119‑1 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                                                >  1:723608/1‑119 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                                  >  1:32885/1‑133 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                                  <  2:32885/133‑1 (MQ=255)
                                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                                            >  2:1353930/1‑139 (MQ=255)
                                                           aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                    >  2:1032760/1‑110 (MQ=255)
                                                           aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                    <  1:1032760/110‑1 (MQ=255)
                                                           aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                              <  1:872122/116‑1 (MQ=255)
                                                           aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                              >  2:872122/1‑116 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGc                                                                                                                        >  2:1023675/1‑80 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGc                                                                                                                        <  1:1023675/80‑1 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGtt                                                                                                              <  2:539610/90‑1 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGtt                                                                                                              >  1:539610/1‑90 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                                                                                               >  1:405412/1‑105 (MQ=255)
                                                                     gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                                                                                               <  2:405412/105‑1 (MQ=255)
                                                                              aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             <  2:140667/66‑1 (MQ=255)
                                                                              aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             >  1:140667/1‑66 (MQ=255)
                                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             <  2:402600/64‑1 (MQ=255)
                                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             >  1:402600/1‑64 (MQ=255)
                                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTc                                                  >  1:1528662/1‑139 (MQ=255)
                                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTc                                                  >  1:1345447/1‑139 (MQ=255)
                                                                                  agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                                             >  2:1178490/1‑62 (MQ=255)
                                                                                  agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAggg                                                                               >  1:735487/1‑108 (MQ=255)
                                                                                  agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAggg                                                                               <  2:735487/108‑1 (MQ=255)
                                                                                  agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGCCCATACAAGCACAGCATACg                                                                                                                             <  1:1178490/62‑1 (MQ=255)
                                                                                   gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc                                                                       >  1:766176/1‑115 (MQ=255)
                                                                                   gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTc                                                                       <  2:766176/115‑1 (MQ=255)
                                                                                   gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGcc                                               >  1:395046/1‑139 (MQ=255)
                                                                                        tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                              >  2:975399/1‑87 (MQ=255)
                                                                                        tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                              <  1:975399/87‑1 (MQ=255)
                                                                                        tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAt                                          <  1:1734108/139‑1 (MQ=255)
                                                                                        tGCGCGTGATCCATATAGGACGCCGCACCTTTCAGGCCATACAAGCACAGAAGACGCATGCCGAGAATATTTTCGCCAATCGCCGCTTTATCTTAGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAt                                          >  2:1651106/1‑139 (MQ=255)
                                                                                         gcgcGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAtc                                         >  2:1068178/1‑139 (MQ=255)
                                                                                            cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAAt                                                                                                     <  1:166411/76‑1 (MQ=255)
                                                                                            cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAAt                                                                                                     >  2:166411/1‑76 (MQ=255)
                                                                                                  ccATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                <  2:395046/139‑1 (MQ=255)
                                                                                                  ccATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                <  2:815302/139‑1 (MQ=255)
                                                                                                         ggCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCGGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGc                         >  2:1469788/1‑139 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                               >  1:933690/1‑36 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                               <  2:933690/36‑1 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgc                                                                                                  <  1:943530/65‑1 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgc                                                                                                  >  2:943530/1‑65 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                     <  2:702866/110‑1 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                     >  1:702866/1‑110 (MQ=255)
                                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCa                        >  2:1544566/1‑139 (MQ=255)
                                                                                                            cgccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAgg                      >  2:420872/1‑139 (MQ=255)
                                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                     >  2:871072/1‑107 (MQ=255)
                                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                     <  1:871072/107‑1 (MQ=255)
                                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                >  2:350417/1‑128 (MQ=255)
                                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                <  1:350417/128‑1 (MQ=255)
                                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt                     <  2:1345447/139‑1 (MQ=255)
                                                                                                              ccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                    <  2:1318686/59‑1 (MQ=255)
                                                                                                              ccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                                    >  1:1318686/1‑59 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATc                                                                                         >  1:1503834/1‑69 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATc                                                                                         <  2:1503834/69‑1 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTc                                                  >  2:1311915/1‑108 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTc                                                  <  1:1311915/108‑1 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                >  1:1513266/1‑126 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc                                <  2:1513266/126‑1 (MQ=255)
                                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCa                   >  1:817514/1‑139 (MQ=255)
                                                                                                                   ccTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCAAGATCGTCGCTCACCAGTTGCAGGTCAGCCa               <  1:679298/139‑1 (MQ=255)
                                                                                                                    cTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCAt              <  1:1605601/139‑1 (MQ=255)
                                                                                                                    cTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCAt              <  2:1528662/139‑1 (MQ=255)
                                                                                                                           gCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCt                                                               >  1:1466182/1‑83 (MQ=255)
                                                                                                                           gCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCt                                                               <  2:1466182/83‑1 (MQ=255)
                                                                                                                                aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGAc  >  1:1469213/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGAC  >  NC_000913/913218‑913484

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