Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,548,360 C→T D320D (GAC→GAT)  fdnG → formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,548,3600CT100.0% 259.4 / NA 79D320D (GAC→GAT) fdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/39);  total (40/39)

CGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACCTGCTG  >  NC_000913/1548231‑1548492
                                                                                                                                 |                                                                                                                                    
cGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTCGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTc                                                                                                                             >  2:1187594/1‑139 (MQ=255)
        aCATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGtt                                                                                                                              <  1:243827/130‑1 (MQ=255)
        aCATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGtt                                                                                                                              >  2:243827/1‑130 (MQ=255)
        aCATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACACCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                                     <  2:1237833/139‑1 (MQ=255)
                  ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                           >  2:216504/1‑139 (MQ=255)
                   cTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                          <  1:263214/139‑1 (MQ=255)
                   cTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                          <  1:1464843/139‑1 (MQ=255)
                      tCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacgacg                                                                                                                <  1:71521/130‑1 (MQ=255)
                      tCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacgacg                                                                                                                >  2:71521/1‑130 (MQ=255)
                        tGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                     <  1:32109/139‑1 (MQ=255)
                          gCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCa                                                                                                   >  2:944744/1‑139 (MQ=255)
                                tGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAACACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGAAGCTGAAAAACGTCAATACGa                                                                                             >  1:110431/1‑139 (MQ=255)
                                    cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                         <  1:396767/139‑1 (MQ=255)
                                    cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                         >  2:934711/1‑139 (MQ=255)
                                     tACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                        <  1:944744/139‑1 (MQ=255)
                                      aCCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATc                                                                                       >  2:366908/1‑139 (MQ=255)
                                      aCCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATc                                                                                       <  2:110431/139‑1 (MQ=255)
                                             cGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                                <  1:109278/139‑1 (MQ=255)
                                                     acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCa                                                                        >  1:199851/1‑139 (MQ=255)
                                                     acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCa                                                                        >  2:1063643/1‑139 (MQ=255)
                                                     acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCa                                                                        >  1:179169/1‑139 (MQ=255)
                                                              aCGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                                >  1:198828/1‑122 (MQ=255)
                                                              aCGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                                <  2:198828/122‑1 (MQ=255)
                                                              aCGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGa                                                               >  2:1537847/1‑139 (MQ=255)
                                                              aCGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGa                                                               <  1:366908/139‑1 (MQ=255)
                                                                gCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaaa                                                             >  2:847182/1‑139 (MQ=255)
                                                                 ccGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaaaa                                                            <  1:934711/139‑1 (MQ=255)
                                                                    aaTACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACgg                                                         >  1:1827675/1‑139 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGc                                                                                                               <  1:1526775/75‑1 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGc                                                                                                               >  2:1526775/1‑75 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCt                                                                                                              <  1:980494/76‑1 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCt                                                                                                              >  2:980494/1‑76 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                        <  1:1669761/98‑1 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                        >  2:1669761/1‑98 (MQ=255)
                                                                              gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGaaa                                               >  1:733912/1‑139 (MQ=255)
                                                                                 ttACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                           >  2:938046/1‑76 (MQ=255)
                                                                                 ttACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                           <  1:938046/76‑1 (MQ=255)
                                                                                  tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                                     >  2:269721/1‑65 (MQ=255)
                                                                                  tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTa                                                                                                                     <  1:269721/65‑1 (MQ=255)
                                                                                  tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                                           >  2:1292378/1‑139 (MQ=255)
                                                                                   acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTg                                                                                                                                >  2:1637401/1‑53 (MQ=255)
                                                                                   acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTg                                                                                                                                <  1:1637401/53‑1 (MQ=255)
                                                                                   acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGAt                                                                 <  2:683573/116‑1 (MQ=255)
                                                                                   acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGAt                                                                 >  1:683573/1‑116 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                           >  2:805543/1‑73 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaa                                                                                                           <  1:805543/73‑1 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGaaa                                               >  1:162161/1‑133 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGaaa                                               <  2:162161/133‑1 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACACTGAAAAACGTCAATACGATaaa                                                                                         <  1:240908/91‑1 (MQ=255)
                                                                                    cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACACTGAAAAACGTCAATACGATaaa                                                                                         >  2:240908/1‑91 (MQ=255)
                                                                                         aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                          >  2:747755/1‑69 (MQ=255)
                                                                                         aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaa                                                                                                          <  1:747755/69‑1 (MQ=255)
                                                                                         aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAac                                    >  1:1008226/1‑139 (MQ=255)
                                                                                         aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAac                                    >  1:1439179/1‑139 (MQ=255)
                                                                                         aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAac                                    >  2:1577808/1‑139 (MQ=255)
                                                                                          cGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                        <  2:1845591/86‑1 (MQ=255)
                                                                                          cGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                        >  1:1845591/1‑86 (MQ=255)
                                                                                              aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGc                                                                                                               <  1:739057/59‑1 (MQ=255)
                                                                                              aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGc                                                                                                               >  2:739057/1‑59 (MQ=255)
                                                                                              aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGa                               <  1:1063643/139‑1 (MQ=255)
                                                                                              aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGa                               <  2:199851/139‑1 (MQ=255)
                                                                                              aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGa                               <  1:216504/139‑1 (MQ=255)
                                                                                                ccTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCt                                                                                  <  2:154019/86‑1 (MQ=255)
                                                                                                ccTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCt                                                                                  >  1:154019/1‑86 (MQ=255)
                                                                                                        tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCACTACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                     <  1:1187594/139‑1 (MQ=255)
                                                                                                        tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                     <  1:737474/139‑1 (MQ=255)
                                                                                                            tgatgaTTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgc                      >  2:1235837/1‑134 (MQ=255)
                                                                                                            tgatgaTTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgc                      <  1:1235837/134‑1 (MQ=255)
                                                                                                              atgatTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                     <  2:354643/53‑1 (MQ=255)
                                                                                                              atgatTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                     >  1:354643/1‑53 (MQ=255)
                                                                                                               tgatTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCtgtgtg              <  2:1008226/139‑1 (MQ=255)
                                                                                                                    ttGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGaa         <  1:1577808/139‑1 (MQ=255)
                                                                                                                      gCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAAcc       >  1:168899/1‑139 (MQ=255)
                                                                                                                       cTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                                  <  1:1824556/47‑1 (MQ=255)
                                                                                                                       cTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                                  >  2:1824556/1‑47 (MQ=255)
                                                                                                                        tttCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACctg     <  1:847182/139‑1 (MQ=255)
                                                                                                                         ttCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGAt                                         <  1:928554/102‑1 (MQ=255)
                                                                                                                         ttCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGAt                                         >  2:928554/1‑102 (MQ=255)
                                                                                                                           cGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACctgctg  >  2:1687949/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACCTGCTG  >  NC_000913/1548231‑1548492

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