Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,495,153 T→C G381G (GGT→GGC)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,495,1530TC100.0% 176.9 / NA 47G381G (GGT→GGC) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/24);  total (23/24)

ACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCT  >  NC_000913/3495053‑3495286
                                                                                                    |                                                                                                                                     
aCCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCtgctgc                                                                                                                             <  1:274122/111‑1 (MQ=255)
aCCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCtgctgc                                                                                                                             >  2:274122/1‑111 (MQ=255)
aCCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGcc                                                                                                 >  1:527691/1‑139 (MQ=255)
  cTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCgg                                                                                               <  2:527691/139‑1 (MQ=255)
           aGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                      <  1:324089/139‑1 (MQ=255)
           aGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                      <  2:266751/139‑1 (MQ=255)
           aGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCgg                                                                                      >  2:646875/1‑139 (MQ=255)
                    aTCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGa                                                                             >  1:494543/1‑139 (MQ=255)
                     tCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGac                                                                            >  2:153601/1‑138 (MQ=255)
                     tCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGAt                                                                            <  1:339130/139‑1 (MQ=255)
                             cGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAg                                                                    <  1:646875/139‑1 (MQ=255)
                              gCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGc                                                                   >  1:49276/1‑139 (MQ=255)
                               ccTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGca                                                                  <  1:89626/139‑1 (MQ=255)
                                 tGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAg                                                                    >  2:740441/1‑135 (MQ=255)
                                 tGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAg                                                                    <  1:740441/135‑1 (MQ=255)
                                            gAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCg                                                     >  1:278693/1‑139 (MQ=255)
                                             aaGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGa                                                                                              <  1:332676/97‑1 (MQ=255)
                                             aaGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGa                                                                                              >  2:332676/1‑97 (MQ=255)
                                              aGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGc                                                   <  2:278693/139‑1 (MQ=255)
                                                 caacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                <  1:525892/139‑1 (MQ=255)
                                                 caacaaAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                <  1:7207/139‑1 (MQ=255)
                                                     aaaaCCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGca                                                                  <  1:56306/117‑1 (MQ=255)
                                                     aaaaCCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGca                                                                  >  2:56306/1‑117 (MQ=255)
                                                      aaaCCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTc                                           <  2:49276/139‑1 (MQ=255)
                                                            tGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCgg                                     <  1:665642/139‑1 (MQ=255)
                                                                   tGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGaa                                                                                                                <  2:519437/57‑1 (MQ=255)
                                                                   tGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGaa                                                                                                                >  1:519437/1‑57 (MQ=255)
                                                                   tGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCg                                      >  1:563817/1‑131 (MQ=255)
                                                                   tGCGGTACTGGTGGGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCg                                      <  2:563817/131‑1 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAAcc                                                                                         >  1:648487/1‑70 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAAcc                                                                                         <  2:648487/70‑1 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATa                             >  2:117962/1‑130 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATa                             <  1:117962/130‑1 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCgg                    >  1:202831/1‑139 (MQ=255)
                                                                             gtggGCGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCgg                    >  1:668509/1‑139 (MQ=255)
                                                                                  cGATACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAgc               >  1:342270/1‑139 (MQ=255)
                                                                                    aTACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGcca                                             <  1:781624/107‑2 (MQ=255)
                                                                                    aTACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGcca                                             >  2:781624/1‑106 (MQ=255)
                                                                                     tACCAGCGACTACGGCCACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                <  1:517215/103‑1 (MQ=255)
                                                                                     tACCAGCGACTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGTAACCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                >  2:517215/1‑103 (MQ=255)
                                                                                            gACTACGGCAACCCGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCa                                                                     <  2:284956/75‑1 (MQ=255)
                                                                                            gACTACGGCAACCCGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCa                                                                     >  1:284956/1‑75 (MQ=255)
                                                                                             aCTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCg                                                     <  2:448212/90‑1 (MQ=255)
                                                                                             aCTACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCg                                                     >  1:448212/1‑90 (MQ=255)
                                                                                               tACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                <  2:369025/93‑1 (MQ=255)
                                                                                               tACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAg                                                >  1:369025/1‑93 (MQ=255)
                                                                                               tACGGCAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCt  >  2:339197/1‑139 (MQ=255)
                                                                                                    |                                                                                                                                     
ACCTCGACGAAAGTAAAGAGATCTACAAACGCCTGATTGTCAGCGAAGACAACAAAACCCTGCTCGGTGCGGTACTGGTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCT  >  NC_000913/3495053‑3495286

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