Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,156 A→T 11.4% intergenic (+34/+67) PP_1661 → / ← PP_1662 dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,1560AT11.4% 124.8 / 7.5 44intergenic (+34/+67)PP_1661/PP_1662dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base T (3/2);  total (25/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.49e-01

CAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCC  >  NC_002947/1859009‑1859303
                                                                                                                                                   |                                                                                                                                                   
cAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATg                                                                                                                                                   >  1:1230913/1‑150 (MQ=255)
      ggTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCggag                                                                                                                                             >  2:649872/1‑150 (MQ=255)
            ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTCTCATGACAACCCCACGATGCTCCGGCAATGCCGCAGCAtttt                                                                                                                                       >  1:441979/1‑150 (MQ=255)
               cAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTggg                                                                                                                                    <  2:810610/150‑1 (MQ=255)
                aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCAAAATGCTGCGGCATTGCCGGAGCAtcgtgggg                                                                                                                                   <  1:649872/150‑8 (MQ=38)
                  cTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGACAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGtt                                                                                                                                 >  1:546608/1‑150 (MQ=255)
                    gtgCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGc                                                                                                                               >  1:69364/1‑150 (MQ=255)
                      gCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGcca                                                                                                                             <  2:1036895/150‑3 (MQ=38)
                           aaTGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGa                                                                                                                        <  2:69364/150‑1 (MQ=255)
                           aaTGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGa                                                                                                                        <  1:612108/150‑1 (MQ=255)
                             tGACGACCTGGACGGGCGTTACGGGGCTGGACTGAGCCATAACAAGAGCAGAGTGTCGGGTGATATGTGGTTCTTTGGGTTGGAGGCCCGTTGATTGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGc                                                                                                                      >  2:162493/1‑150 (MQ=255)
                                        aCGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTc                                                                                                           >  1:59844/1‑150 (MQ=255)
                                            gCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAggg                                                                                                       >  2:390011/1‑150 (MQ=255)
                                                            cTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                                       >  1:610799/1‑150 (MQ=255)
                                                             tttgCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                                      >  1:410247/4‑150 (MQ=255)
                                                             tGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                                      >  1:215605/1‑150 (MQ=255)
                                                             tGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCa                                                                                      >  1:110140/1‑150 (MQ=255)
                                                                  gATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGACGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGGGTCAGGGGTTCGCCGCAATGGGCAgggac                                                                                 >  2:303760/1‑148 (MQ=255)
                                                                        cTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtg                                                                           >  1:181723/1‑150 (MQ=255)
                                                                          gggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                         <  1:564970/150‑1 (MQ=255)
                                                                                         agagGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTcc                                                          <  1:390011/150‑1 (MQ=255)
                                                                                          gagGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCa                                                         <  2:59844/150‑1 (MQ=255)
                                                                                             gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                         <  2:97864/99‑1 (MQ=255)
                                                                                             gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                         >  1:97864/1‑99 (MQ=255)
                                                                                                  tttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCa                                                 >  2:554170/1‑150 (MQ=255)
                                                                                                   ttCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCaa                                                <  1:29927/150‑1 (MQ=255)
                                                                                                    tCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAg                                               <  2:545152/150‑1 (MQ=255)
                                                                                                              gtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                    <  2:1219125/103‑1 (MQ=255)
                                                                                                              gtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                    >  1:1219125/1‑103 (MQ=255)
                                                                                                                   cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                    <  1:1068925/98‑1 (MQ=255)
                                                                                                                   cccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAgg                                                                                    >  2:1068925/1‑98 (MQ=255)
                                                                                                                    ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                               >  1:295061/1‑150 (MQ=255)
                                                                                                                    ccTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTgg                               >  1:1127244/1‑150 (MQ=255)
                                                                                                                       gTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGgtgt                            <  2:1165251/150‑1 (MQ=255)
                                                                                                                            agCAACCCCAAAATGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                       >  1:1341025/2‑150 (MQ=37)
                                                                                                                                 ccccaaaaTGCTGCGGCATTGCCGGAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                  >  2:814138/8‑150 (MQ=38)
                                                                                                                                 ccccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                  <  2:295061/150‑1 (MQ=255)
                                                                                                                                  cccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                 >  2:1054508/1‑150 (MQ=255)
                                                                                                                                  cccACGATGCTCCGGAAATGCCGTACCAGTTTGGGGTTGTTCACTGAGACCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTTAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                 <  2:582263/150‑1 (MQ=255)
                                                                                                                                   ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                       <  1:1401064/143‑1 (MQ=255)
                                                                                                                                   ccACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                       >  2:1401064/1‑143 (MQ=255)
                                                                                                                                    cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTgg                                                                        <  2:1375529/93‑1 (MQ=255)
                                                                                                                                    cACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTgg                                                                        >  1:1375529/1‑93 (MQ=255)
                                                                                                                                                 cAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCACTGGGCAGGGTAATTGAGCGGCGCGAGCGTTTTCCAGTTGTTCAAGCTCGACTCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGcc  <  2:1230913/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCGCAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCC  >  NC_002947/1859009‑1859303

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