Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,259 Δ1 bp intergenic (‑20/+28) yqeB ← / ← yqeC XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold/putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,2580C.100.0% 303.1 / NA 74intergenic (‑19/+29)yqeB/yqeCXdhC‑CoxI family protein with NAD(P)‑binding Rossman fold/putative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (37/37);  total (37/37)

CACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCT  >  NC_000913/3014126‑3014385
                                                                                                                                    |                                                                                                                               
cACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAg                                                                                                                          <  2:1030351/139‑1 (MQ=255)
        tttGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCgg                                                                                                                  <  1:480937/139‑1 (MQ=255)
           gtgCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCg                                                                                                                   >  2:455761/1‑135 (MQ=255)
           gtgCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCg                                                                                                                   <  1:455761/135‑1 (MQ=255)
           gtgCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAAt                                                                                                               >  1:607506/1‑139 (MQ=255)
           gtgCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGCGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACACGGAAt                                                                                                               >  1:1386273/1‑139 (MQ=255)
             gCAGAATGGCGGGGAGTCGAGCCTAGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATAc                                                                                                             >  2:1121373/1‑139 (MQ=255)
              cAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       >  1:1137708/1‑128 (MQ=255)
              cAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       <  2:1137708/128‑1 (MQ=255)
                   tGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       >  2:676361/1‑123 (MQ=255)
                   tGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       <  1:676361/123‑1 (MQ=255)
                               gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                           >  1:1202519/1‑139 (MQ=255)
                                     cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                     >  1:1069391/1‑139 (MQ=255)
                                          aTCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       <  2:8924/100‑1 (MQ=255)
                                          aTCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTaa                                                                                                                       >  1:8924/1‑100 (MQ=255)
                                           tCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGcgcgcg                                                                               <  2:1069391/139‑1 (MQ=255)
                                                aaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                  <  2:329084/115‑4 (MQ=255)
                                                aaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                  >  1:329084/1‑115 (MQ=255)
                                                 aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                         <  1:747087/139‑1 (MQ=255)
                                                 aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                         <  1:343675/139‑1 (MQ=255)
                                                  tttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCa                                                                                                                           <  2:1108078/88‑1 (MQ=255)
                                                  tttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGATTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCa                                                                                                                           >  1:1108078/1‑88 (MQ=255)
                                                    tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                      <  1:100534/139‑1 (MQ=255)
                                                    tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                      <  2:1397341/139‑1 (MQ=255)
                                                    tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGCGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                      <  2:1573777/139‑1 (MQ=255)
                                                     gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAg                                                                                                                          <  1:1181445/86‑1 (MQ=255)
                                                     gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAg                                                                                                                          >  2:1181445/1‑86 (MQ=255)
                                                     gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                     >  2:656665/1‑139 (MQ=255)
                                                     gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                     >  2:1202888/1‑139 (MQ=255)
                                                     gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                     >  1:665385/1‑139 (MQ=255)
                                                      cgacATAGCAAACAGACAATTTTGCTCTTCGAGGTTTGCAGCCTCTGTGAAAATATTCATACCAACGTCAATGGATATCTTTCAGTAACCCGGAAAACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                    <  1:189966/136‑1 (MQ=255)
                                                       gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTa                                                                                                                        >  1:456987/1‑86 (MQ=255)
                                                       gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTa                                                                                                                        <  2:456987/86‑1 (MQ=255)
                                                           tAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                  >  1:671242/1‑101 (MQ=255)
                                                           tAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                  <  2:671242/104‑4 (MQ=255)
                                                           tAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                               >  1:364365/1‑139 (MQ=255)
                                                                   ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAta                                                                    <  2:1387744/126‑2 (MQ=255)
                                                                   ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAta                                                                    >  1:1387744/1‑125 (MQ=255)
                                                                     aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                           <  1:208805/117‑1 (MQ=255)
                                                                     aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                           >  2:208805/1‑117 (MQ=255)
                                                                            tGCTCTTCGAGTTTTGCAGCCTCTGTTAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCctagccag                                              <  2:106308/139‑1 (MQ=255)
                                                                             gCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGGATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                                             >  1:860788/1‑139 (MQ=255)
                                                                                   cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                       <  2:1028363/139‑1 (MQ=255)
                                                                                         ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgt                                 <  1:1483620/139‑1 (MQ=255)
                                                                                           gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                               <  2:607506/139‑1 (MQ=255)
                                                                                                   tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                       >  1:1408285/1‑91 (MQ=255)
                                                                                                   tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                       <  2:1408285/91‑1 (MQ=255)
                                                                                                   tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                       >  2:259801/1‑139 (MQ=255)
                                                                                                     tgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                     <  2:460260/139‑1 (MQ=255)
                                                                                                     cgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTTTTGTTGCa                     <  2:1386273/138‑1 (MQ=255)
                                                                                                             ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                 >  1:256233/1‑71 (MQ=255)
                                                                                                             ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                 <  2:256233/71‑1 (MQ=255)
                                                                                                             ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                 >  2:445829/1‑103 (MQ=255)
                                                                                                             ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                 <  1:445829/103‑1 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCgg                                                                                                                  >  1:1171798/1‑36 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCgg                                                                                                                  <  2:1171798/36‑1 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg              <  2:178290/136‑1 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTg              >  1:178290/1‑136 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg           >  1:1432976/1‑139 (MQ=255)
                                                                                                               cATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGAAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAg           >  1:478851/1‑139 (MQ=255)
                                                                                                                 tACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGa         >  1:1270668/1‑139 (MQ=255)
                                                                                                                  aCCAACCTCAATGGATATCTTTCAGTAACCCGGGATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg        <  2:860788/139‑1 (MQ=255)
                                                                                                                   ccAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                      >  2:1651744/1‑76 (MQ=255)
                                                                                                                   ccAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                      <  1:1651744/76‑1 (MQ=255)
                                                                                                                      aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc       <  2:1226953/136‑1 (MQ=255)
                                                                                                                      aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc       >  1:1226953/1‑136 (MQ=255)
                                                                                                                      aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg    >  2:1655605/1‑139 (MQ=255)
                                                                                                                       ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg        >  1:269576/1‑134 (MQ=255)
                                                                                                                       ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg        <  2:269576/134‑1 (MQ=255)
                                                                                                                        cTCGATGGATATCTTTCAGTGACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCTGTTGAAGGAGCTCGCt  >  1:1123431/1‑139 (MQ=255)
                                                                                                                        cTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                  <  1:1621688/43‑4 (MQ=255)
                                                                                                                        cTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                  >  2:1621688/1‑43 (MQ=255)
                                                                                                                            aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                   <  2:597827/54‑1 (MQ=255)
                                                                                                                            aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                   >  1:597827/1‑54 (MQ=255)
                                                                                                                                    |                                                                                                                               
CACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCT  >  NC_000913/3014126‑3014385

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