Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,172,573:1 +C 100% coding (931/1083 nt) PP_1925 ← monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,172,5731.C92.7% 130.2 / 4.0 41coding (931/1083 nt)PP_1925monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (20/18);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TCCGGCTGCTCCAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑C‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGAGCCAAGCGACCCCGTCC  >  NC_002947/2172428‑2172712
                                                                                                                                                  |                                                                                                                                            
tCCGGCTGCTCCAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑c‑cc                                                                                                                                          <  7:300563/149‑1 (MQ=255)
      tgctCCAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑c‑ccccccat                                                                                                                                    >  1:186563/1‑146 (MQ=255)
         tCCAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCaa                                                                                                                                  >  5:54155/1‑149 (MQ=255)
           cAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAAc                                                                                                                                 <  5:171785/148‑1 (MQ=255)
                 gTTCTGGTGTTCTACTACTTCGTTGACCGTGCTGCGCGCGTTACGGGTGACGCCGATCAAGGTCGCCGAGGGCAGGCCGGTCCAGTCGCCATATCCGACGAGCCTGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGAc                                                                                                                           <  4:251525/148‑1 (MQ=255)
                                 actTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑C‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAg                                                                                                         >  8:199911/1‑149 (MQ=255)
                                  ctTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAg                                                                                                         >  8:79742/1‑149 (MQ=255)
                                       ttGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGt                                                                                                    <  6:10262/149‑1 (MQ=255)
                                          aCCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTc                                                                                                   >  6:89469/1‑147 (MQ=255)
                                          aCCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCa                                                                                                 <  5:89469/149‑1 (MQ=255)
                                           ccGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑C‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAgt                                                                                               <  7:199911/149‑1 (MQ=255)
                                                  gcgcgcGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAgtgt                                                                                             <  3:281919/145‑1 (MQ=255)
                                                  gcgcgcGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAgtgt                                                                                             >  4:281919/1‑145 (MQ=255)
                                                   cgcgcgGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCa                                                                                        <  2:186563/149‑1 (MQ=255)
                                                            cGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCa                                                                               <  6:54155/149‑1 (MQ=255)
                                                                    gCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTc                                                                                            >  1:121571/1‑128 (MQ=255)
                                                                    gCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTc                                                                                            <  2:121571/128‑1 (MQ=255)
                                                                    gCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGc                                                                       >  3:51507/1‑149 (MQ=255)
                                                                     ccGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCa                                                                      <  7:14366/149‑1 (MQ=255)
                                                                          cAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATTGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGc                                                                       >  4:271062/1‑143 (MQ=255)
                                                                          cAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATTGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGc                                                                       <  3:271062/143‑1 (MQ=255)
                                                                           aaCGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAAcc                                                                <  7:79742/149‑1 (MQ=255)
                                                                            aCGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAgt                                                                                               <  6:89784/117‑1 (MQ=255)
                                                                            aCGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAgt                                                                                               >  5:89784/1‑117 (MQ=255)
                                                                                       gccaggccGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGa                                                    >  5:30878/1‑149 (MQ=255)
                                                                                       gccaggccGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGa                                                    >  5:327499/1‑149 (MQ=255)
                                                                                            gccGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACcaca                                               <  4:51507/149‑1 (MQ=255)
                                                                                                   cAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCa                                                      <  8:266315/135‑1 (MQ=255)
                                                                                                   cAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCa                                                      >  7:266315/1‑135 (MQ=255)
                                                                                                    aGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGAc                                       >  2:290989/1‑149 (MQ=255)
                                                                                                    aGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGAc                                       >  2:26099/1‑149 (MQ=255)
                                                                                                      tCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCt                                     <  6:30878/149‑1 (MQ=255)
                                                                                                        gCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTgg                                   <  1:290989/149‑1 (MQ=255)
                                                                                                             atCCGACGAGCCAGAGGCGCGGCTAGCGAATAGAGTG‑CACCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGAGCAGGCAGGAACCCAGTACACCAGACCAAATCATCTACCTgtgtat                               >  6:267894/1‑146 (MQ=255)
                                                                                                                cgacgaGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTGGGTTTtctc                           <  1:26099/149‑1 (MQ=255)
                                                                                                                     aGCCAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTCCCGGTGTCATCCAGCACGGCCAGGGTCGTCAGCTGATCCAGCGCCGGCAGGAACCCAGTACCCCAGACCCCAGCATCGACCTGGGGTTTCTCCCCGt                      >  1:189325/1‑149 (MQ=255)
                                                                                                                        cAGAGGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGa                   <  5:267894/149‑1 (MQ=255)
                                                                                                                          gagGCGCGGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACCCCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGAGc                 >  7:96697/1‑149 (MQ=255)
                                                                                                                               gcgGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGTATCGACCTGGGTTTTCTCACCGTTGAGCCAAGc            >  5:58454/1‑149 (MQ=255)
                                                                                                                                cgGCTCGCGAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACAAGCCCACAGCATCGACCTGGGTTTTATCTCCGTTGCGCCAAGCg           >  7:153550/1‑149 (MQ=255)
                                                                                                                                      gcgAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGAGCCAAGCGACCCCg     >  7:133581/1‑149 (MQ=255)
                                                                                                                                      gcgAATAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTATTCAGATGATCCAGCGCCGGCATGAACCAAGTACACCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGAGCCAAGCGACCCCg     >  3:60820/1‑149 (MQ=255)
                                                                                                                                         aaTAGAGTGCC‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGCCCACAGCATCGACCTGGGTTTTCTCACCGTTGAGCCAAGCGACCCCGTcc  >  3:98/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TCCGGCTGCTCCAAGCTGTTCAGGTGTTCAACTACTTCGTTGACCGTGCTGCGCGCGGTACGGGTGACGCCGATCAACGTCGCCGAGGCCAGGCCGGTCCAGTCGCCATATCCGACGAGCCAGAGGCGCGGCTCGCGAATAGAGTG‑C‑CCCCCATCAACCTCGACCCTACCGTTGTCATCCAGCACGTCCAGTGTCTTCAGATGATCCAGCGCCGGCAGGAACCCAGTACACCAGACCACAGCATCGACCTGGGTTTTCTCACCGTTGAGCCAAGCGACCCCGTCC  >  NC_002947/2172428‑2172712

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