Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,581 C→T 15.4% intergenic (‑428/+18) PP_0022 ← / ← PP_0023 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294727,5810CT15.4% 84.6 / 9.2 39intergenic (‑428/+18)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (4/2);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG  >  NC_002947/27440‑27728
                                                                                                                                               |                                                                                                                                                   
tcgtcgTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAATCAGCTCCCGAGAgct                                                                                                                                             >  1:262915/1‑150 (MQ=255)
           acggAAGGGATAACCCATCCGGAGGTAGGTTCAGAATCCACCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGGGAGACGGTCGTTAGGGTTTGCAAACAGCCCCAAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTc                                                                                                                                    <  1:343003/147‑1 (MQ=255)
               aaGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGTCTCACGACGCGGTCCGATTGATCAGCGAGAGGGTCGTGAGGGTTTGCACAAAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc                                                                                                                              <  1:436755/150‑1 (MQ=255)
                        ccGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCtttt                                                                                                                      >  1:40361/1‑150 (MQ=255)
                         cGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTtg                                                                                                                    >  1:171723/1‑150 (MQ=255)
                                gAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtgt                                                                                                             <  1:532571/150‑1 (MQ=255)
                                        ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCttttt                                                                                                                     >  1:482376/1‑134 (MQ=255)
                                        ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCttttt                                                                                                                     <  2:482376/134‑1 (MQ=255)
                                         tCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                     >  1:126836/1‑150 (MQ=255)
                                             aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa                                                                                                                                 <  2:403589/117‑1 (MQ=255)
                                             aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa                                                                                                                                 >  1:403589/1‑117 (MQ=255)
                                                   tccagccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGCCGGTCGTCAGGGATTGCACACAGCCCC‑‑AAAAAAAAAGCATCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg                                                                                          <  1:169808/147‑1 (MQ=255)
                                                          aagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                   >  2:422592/1‑150 (MQ=255)
                                                               gAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCtt                                                                              <  1:15267/150‑1 (MQ=255)
                                                                    cTCACGACGCGGTCCTATTGAGCAGCGAGACGGTCGTGAGGGTTTGCACACAGCGCC‑‑GAAAAGAAAGTAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAg                                                                         <  2:552714/150‑1 (MQ=255)
                                                                       acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                       <  2:126836/150‑1 (MQ=255)
                                                                          acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg                                                                        >  2:220790/1‑145 (MQ=255)
                                                                          acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgg                                                                        <  1:220790/145‑1 (MQ=255)
                                                                               gTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg                                                                                          >  1:594531/1‑122 (MQ=255)
                                                                               gTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg                                                                                          <  2:594531/122‑1 (MQ=255)
                                                                                           aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                   <  1:519927/101‑1 (MQ=255)
                                                                                           aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                   >  2:519927/1‑101 (MQ=255)
                                                                                                      cGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATc                                       <  1:422592/150‑1 (MQ=255)
                                                                                                              tttGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                   >  1:561490/1‑98 (MQ=255)
                                                                                                              tttGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                   <  2:561490/98‑1 (MQ=255)
                                                                                                                   acacaGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTgcgc                          >  1:360449/1‑150 (MQ=255)
                                                                                                                     acaGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCCTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCtt                        >  2:351036/1‑150 (MQ=255)
                                                                                                                         ccca‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTggg                    >  1:452759/1‑150 (MQ=255)
                                                                                                                              taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                     <  2:708961/142‑1 (MQ=39)
                                                                                                                              taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                     >  1:708961/5‑146 (MQ=39)
                                                                                                                              taagaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTTTTTTTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                 >  2:505986/5‑150 (MQ=255)
                                                                                                                                agaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg               >  1:112647/3‑150 (MQ=255)
                                                                                                                                aaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg               >  1:158843/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGTCTGGc              >  1:555197/1‑150 (MQ=255)
                                                                                                                                     aaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTTCGCTTTGGGCCTGGCTTGa          >  2:282587/1‑150 (MQ=255)
                                                                                                                                      aaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGAt         <  2:171723/150‑1 (MQ=255)
                                                                                                                                          agcTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTc     <  2:112647/150‑1 (MQ=255)
                                                                                                                                            cTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCAt   <  2:360449/150‑1 (MQ=255)
                                                                                                                                             tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGGGGGGTTTCAGGTCCGGATGGGCCTTCTGCGTCTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCGTGATGCTCGTg  >  1:100969/1‑150 (MQ=255)
                                                                                                                                               ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGTTGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                       >  1:548968/1‑79 (MQ=255)
                                                                                                                                               ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                       <  2:548968/79‑1 (MQ=255)
                                                                                                                                               ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc                                                           >  1:674831/1‑91 (MQ=255)
                                                                                                                                               ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc                                                           <  2:674831/91‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG  >  NC_002947/27440‑27728

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: