Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,522,286 A→G P782P (CCT→CCC)  nirB ← nitrite reductase (NAD(P)H) large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,522,2860AG96.3% 84.3 / NA 27P782P (CCT→CCC) nirBnitrite reductase (NAD(P)H) large subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/15);  minor base T (0/1);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCTTTTTTCCAATTAGGGGTGTTAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAAGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCTGTTCTTTCAG  >  CP000730/2522154‑2522416
                                                                                                                                    |                                                                                                                                  
aaCCTTTTTTCCAATTAGGGGTGTTAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGc                                                                                                                            >  1:309909/1‑141 (MQ=255)
      ttttCCAATTAGGGGTGTTAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGAc                                                                                                                      >  1:268614/1‑141 (MQ=255)
                     tgtTAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATc                                                                                                       >  2:397383/1‑141 (MQ=255)
                        tAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCAt                                                                                                    <  1:399191/141‑1 (MQ=255)
                        tAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCAt                                                                                                    <  2:100260/141‑1 (MQ=255)
                            tCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGAGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATg                                                                                                >  1:288168/1‑141 (MQ=255)
                                   taGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCaa                                                                                         >  1:364741/1‑141 (MQ=255)
                                          cACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTc                                                                                  <  2:111854/141‑1 (MQ=255)
                                              ttAATTTTCTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCAttt                                                                              <  2:319171/141‑1 (MQ=255)
                                                aaTTTTTTTTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCATGGTTCAGCTTCGACTGCTTTTTTAGCACCCAGAATGCGTTCAAATAGTTCATTTTg                                                                            <  2:357731/141‑1 (MQ=255)
                                                aaTTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTg                                                                            <  2:271050/141‑1 (MQ=255)
                                                      ttCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTc                                                                      >  1:369747/1‑141 (MQ=255)
                                                         ttttGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTgg                                                                   <  2:282413/141‑1 (MQ=255)
                                                                   aTGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAg                                                         >  1:50626/1‑141 (MQ=255)
                                                                    tGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGa                                                        <  2:248254/141‑1 (MQ=255)
                                                                     gTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGAc                                                       <  1:317030/141‑1 (MQ=255)
                                                                           cacCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCttt                                                 <  2:268614/141‑1 (MQ=255)
                                                                                tctcCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACAt                                            <  1:397383/141‑1 (MQ=255)
                                                                                     aCTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCa                                       >  1:196778/1‑141 (MQ=255)
                                                                                     aCTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCa                                       <  2:296110/141‑1 (MQ=255)
                                                                                         cAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATc                                   <  2:267836/141‑1 (MQ=255)
                                                                                               ttCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGtc                             <  2:151610/141‑1 (MQ=255)
                                                                                                     gTGCTTGCGCATTGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTa                       >  1:245892/1‑141 (MQ=255)
                                                                                                                ttGCTAGTTATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCt            <  2:349594/141‑1 (MQ=255)
                                                                                                                        tATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCTGTTCTTTc    >  1:180541/1‑141 (MQ=255)
                                                                                                                        tATTGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCTGTTCTTTc    >  1:180536/1‑141 (MQ=255)
                                                                                                                          ttGCTTCCCAGGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCTGTTCTTTCAg  <  1:187019/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
AACCTTTTTTCCAATTAGGGGTGTTAATTCATCTATAGTTGTCACTTTAATTTTTTCTTTTGTTTCCATGTATTACACCTTCTCCACTTCAAAAATTCTACGTGCTTGCGCATTGCTAGTTATTGCTTCCCAAGGTTCAGCTTCGACTGCTTTTTTAGCATCCATAATGCGTTCAAATAGTTCATTTTGTCTTTCTGGGTCAAGTAAGACTTCTTTTACATTTTCAAATCCAAGTCTTCTTAACCATGGCGCTGTTCTTTCAG  >  CP000730/2522154‑2522416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑