Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,617,688 T→C V225A (GTA→GCA)  nikD → Ni(2(+)) ABC transporter ATP binding subunit NikD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,617,6880TC100.0% 171.0 / NA 50V225A (GTA→GCA) nikDNi(2(+)) ABC transporter ATP binding subunit NikD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/36);  total (14/36)

TCGATCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGTAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATGCGCACGGT  >  NC_000913/3617553‑3617823
                                                                                                                                       |                                                                                                                                        
tCGATCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGaa                                                                                                                                      >  2:317648/1‑140 (MQ=255)
 cGATCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGaaa                                                                                                                                     <  1:700565/140‑1 (MQ=255)
   aTCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACg                                                                                                                                   <  1:1006172/140‑1 (MQ=255)
   aTCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACg                                                                                                                                   <  1:85930/140‑1 (MQ=255)
        ctgGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGtt                                                                                                                             >  1:812864/1‑140 (MQ=255)
        ctgGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGGGATGTCTTACGGGAAGATTGTCGAACAGGGGGATGCAGAAAAGCTGtt                                                                                                                             >  2:1180596/1‑140 (MQ=255)
         tgGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGGAAGATTGTCGAACAGGGGGATGCAGAAAAGCTGttt                                                                                                                            >  2:691629/1‑140 (MQ=255)
             aaGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACg                                                                                                                        <  2:1178067/140‑1 (MQ=255)
              aGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGGTGTCTGACGGTAAGATTGTCGAACAGGGGGGAGGAGAAAAGCTGttttaagc                                                                                                                       >  2:85930/1‑135 (MQ=255)
              aGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGc                                                                                                                       >  2:728703/1‑140 (MQ=255)
                 aTTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGcccc                                                                                                                    <  2:812864/140‑1 (MQ=255)
                 aTTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGcccc                                                                                                                    <  1:13867/140‑1 (MQ=255)
                  ttATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGccccc                                                                                                                   >  1:1035335/1‑140 (MQ=255)
                  ttATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGccccc                                                                                                                   <  2:995803/140‑1 (MQ=255)
                  ttATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGccccc                                                                                                                   >  1:1050577/1‑140 (MQ=255)
                    aTGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCaa                                                                                                                 <  1:691629/140‑1 (MQ=255)
                     tGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCaaa                                                                                                                >  1:510568/1‑140 (MQ=255)
                      gCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAAc                                                                                                               <  1:96103/140‑1 (MQ=255)
                          aaaCAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGGAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAAAATAc                                                                                                           >  1:1038264/1‑140 (MQ=255)
                                 cgcCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGAcg                                                                                                    <  2:123682/140‑1 (MQ=255)
                                   ccGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGAcgcg                                                                                                  <  2:817605/140‑1 (MQ=255)
                                          tgctgcTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTg                                                                                           <  1:1180596/140‑1 (MQ=255)
                                          tgctgcTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGGTATTTTTGACGGGAAGAATGTCGAAAAAGGCGGTGCAGGAAAACTGTTTAAAGCCCCCAAAAATAAAGAGACGGGCCGCCTg                                                                                           >  2:277909/1‑140 (MQ=255)
                                            ctgctgGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTTGCGGGAAGATTGTCGAACACGGCGATGCAGAAACGCTGTTTTAAGCCCCCCAAAATAAAGAGACGCGCCGCCTggg                                                                                         >  2:572855/1‑139 (MQ=255)
                                            ctgctgGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGt                                                                                         <  1:317648/140‑1 (MQ=255)
                                                  gTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGGGATGTTTGACGGTAAGATTGTCGAACAGGGCGAGGCAGAAAAGCTGGTTTACGCCCCCAAACAAACAGTGACGCGCAGCCCGGTTTCCGc                                                                                   >  2:249993/1‑140 (MQ=255)
                                                    agcccATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTc                                                                                 <  1:556737/138‑1 (MQ=255)
                                                      cccATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCCAAAATACAGGGACGCGCAGCCTGGTTTTCGCCCAt                                                                               >  2:109729/1‑140 (MQ=255)
                                                       ccATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCAtc                                                                              <  2:510568/140‑1 (MQ=255)
                                                        cATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCAtct                                                                             <  1:728703/140‑1 (MQ=255)
                                                         aTGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCAtctc                                                                            <  1:572855/140‑1 (MQ=255)
                                                          tGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCg                                                                           <  2:1092871/140‑1 (MQ=255)
                                                          tGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCAAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCg                                                                           <  2:633520/140‑1 (MQ=255)
                                                          tGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGGAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCCAACATACAGGGACGCGCAGCCTGGGTTCCGCTCAtctcc                                                                           >  2:144416/1‑139 (MQ=255)
                                                                 ggCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTa                                                                    <  1:912670/140‑1 (MQ=255)
                                                                 ggCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTa                                                                    <  1:306179/140‑1 (MQ=255)
                                                                 ggCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTa                                                                    <  1:464191/140‑1 (MQ=255)
                                                                 ggCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGGAAGAATGTCGAACAAGGCGGTGCAGAAACGCTGTTTTAACCCCCCCAAAATACAATGACGCGCAGCCCGGGTTCCGCTCATCTCGCCctccc                                                                    >  2:556737/1‑138 (MQ=255)
                                                                    gTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGGAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAAAATAAAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACgg                                                                 >  1:10620/1‑140 (MQ=255)
                                                                     ttGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGGCGGGAAGATTGTCGGACAGGGCGATGCAGAAAAGCTGTTTAACGCCCCCAAAAATACAGTGACGCGCAGCCTGGTTTTCGCTCATCTCGCCCCCTACggg                                                                >  2:99316/1‑139 (MQ=255)
                                                                     ttGTGGCGCGTCTGGCGGATGACGTGGCGGGGATGTTTTACCGGAAGAATGTCGAACAGGGCGATGCAGAAAAGCTGTTTTACGCCCCCAAAAATACAAAGACGCGCAGCCTGGGTTTCCCTCATCTCCCCctcctcggg                                                                >  2:912670/1‑134 (MQ=255)
                                                                      tgtgGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTa                                                               <  1:144416/140‑1 (MQ=255)
                                                                       gtgGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGccc                                                                        >  1:136082/1‑130 (MQ=255)
                                                                       gtgGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGccc                                                                        <  2:136082/130‑1 (MQ=255)
                                                                       gtgGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTa                                                                    <  2:612006/134‑1 (MQ=255)
                                                                       gtgGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTa                                                                    >  1:612006/1‑134 (MQ=255)
                                                                              gTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTg                                                       <  1:705135/140‑1 (MQ=255)
                                                                                  ggCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGcat                                                   <  2:103918/140‑1 (MQ=255)
                                                                                                 ggTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTaa                                    <  2:1035335/140‑1 (MQ=255)
                                                                                                 ggTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTaa                                    <  2:1050577/140‑1 (MQ=255)
                                                                                                           gACGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGc                          >  1:1127712/1‑140 (MQ=255)
                                                                                                            aCGGTAAGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGcc                         <  1:187118/140‑1 (MQ=255)
                                                                                                                 aaGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTcc                    <  1:249993/140‑1 (MQ=255)
                                                                                                                  aGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGcc                         >  1:1055915/1‑134 (MQ=255)
                                                                                                                  aGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGcc                         <  2:1055915/134‑1 (MQ=255)
                                                                                                                  aGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTccc                   >  1:1069591/1‑140 (MQ=255)
                                                                                                                  aGATTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTccc                   <  2:959915/140‑1 (MQ=255)
                                                                                                                    aTTGTCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCAt                 <  1:99316/140‑1 (MQ=255)
                                                                                                                        tCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACt             >  2:1193931/1‑140 (MQ=255)
                                                                                                                        tCGAACAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACt             >  2:1008830/1‑140 (MQ=255)
                                                                                                                            aCAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATgc         >  1:246166/1‑140 (MQ=255)
                                                                                                                             cAGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATgcg        >  2:490655/1‑140 (MQ=255)
                                                                                                                              aGGGCGATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATgcgc       >  2:42789/1‑140 (MQ=255)
                                                                                                                                   gATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATGCGCACGGt  <  1:125364/140‑1 (MQ=255)
                                                                                                                                   gATGCAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATGCGCACGGt  <  2:1061504/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TCGATCTGCTGGAAAGCATTATGCAAAAACAAGCGCCGGGAATGCTGCTGGTGACCCATGATATGGGCGTTGTGGCGCGTCTGGCGGATGACGTGGCGGTGATGTCTGACGGTAAGATTGTCGAACAGGGCGATGTAGAAACGCTGTTTAACGCCCCCAAACATACAGTGACGCGCAGCCTGGTTTCCGCTCATCTCGCCCTCTACGGTATGGAGCTGGCATCATGACTTTACTTAACATCTCCGGCCTTTCCCATCACTATGCGCACGGT  >  NC_000913/3617553‑3617823

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