Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,245,600 G→A *147* (TAG→TAA ycgY → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,6000GA94.4% 107.2 / ‑0.0 36coding (441/441 nt)ycgYuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/16);  minor base . (1/1);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACT  >  NC_000913/1245468‑1245728
                                                                                                                                    |                                                                                                                                
ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaa                                                                                                                              <  2:942181/137‑1 (MQ=255)
ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaa                                                                                                                              >  1:942181/1‑137 (MQ=255)
ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaa                                                                                                                               >  1:1081803/1‑136 (MQ=255)
ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaa                                                                                                                               <  2:1081803/136‑1 (MQ=255)
        ggCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATCACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaa                                                                                                                               <  1:1209190/128‑1 (MQ=255)
        ggCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATCACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaa                                                                                                                               >  2:1209190/1‑128 (MQ=255)
                    cGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCt                                                                                                        >  2:655178/1‑139 (MQ=255)
                    cGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCt                                                                                                        >  2:426209/1‑139 (MQ=255)
                              aaGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                              <  2:1123560/139‑1 (MQ=255)
                                    aaCACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaa                                                                                                                              <  2:135750/101‑1 (MQ=255)
                                    aaCACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaa                                                                                                                              >  1:135750/1‑101 (MQ=255)
                                                     aTGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACg                                                                       <  2:718909/139‑1 (MQ=255)
                                                     aTGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACg                                                                       >  1:746908/1‑139 (MQ=255)
                                                                        ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTgg                                                    >  1:1108977/1‑139 (MQ=255)
                                                                        ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTgg                                                    >  2:70857/1‑139 (MQ=255)
                                                                         cACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGt                                                   <  1:426209/139‑1 (MQ=255)
                                                                          aCGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACTCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                  <  2:951656/139‑1 (MQ=255)
                                                                          aCGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                  <  1:907016/139‑1 (MQ=255)
                                                                                  cAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGtgtg                                          >  2:206265/1‑139 (MQ=255)
                                                                                    aaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGt                                       >  2:409456/1‑139 (MQ=255)
                                                                                     aCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGt                                       >  1:998829/1‑139 (MQ=255)
                                                                                      cATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGgag                                                                                    >  1:1217054/1‑93 (MQ=255)
                                                                                      cATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGgag                                                                                    <  2:1217054/93‑1 (MQ=255)
                                                                                        tttGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                    <  1:70857/139‑1 (MQ=255)
                                                                                          tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATg                                                                                                                 >  1:711225/1‑60 (MQ=255)
                                                                                          tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATg                                                                                                                 <  2:711225/60‑1 (MQ=255)
                                                                                                aTACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTTCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCttt                            >  2:129942/1‑139 (MQ=255)
                                                                                                 tACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTg                           >  1:948722/1‑139 (MQ=255)
                                                                                                      cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTa                                                                                                                   <  1:239629/46‑1 (MQ=255)
                                                                                                      cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTa                                                                                                                   >  2:239629/1‑46 (MQ=255)
                                                                                                      cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                    >  2:1183715/1‑125 (MQ=255)
                                                                                                      cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                    <  1:1183715/125‑1 (MQ=255)
                                                                                                         gCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                         <  1:629026/37‑1 (MQ=38)
                                                                                                         gCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                         >  2:629026/1‑37 (MQ=255)
                                                                                                         gCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAggg                   >  1:1229060/1‑139 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGc                                  >  1:285508/1‑122 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGc                                  <  2:285508/122‑1 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCctc                               >  2:113193/1‑125 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCctc                               <  1:113193/125‑1 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCtt                             <  1:1151340/127‑1 (MQ=255)
                                                                                                           cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCtt                             >  2:1151340/1‑127 (MQ=255)
                                                                                                                 ttGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTCTGGTTGAGGGTTGcgtcg           <  2:746908/139‑1 (MQ=255)
                                                                                                                     aaaaCCGTCTGTGTAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCt      <  1:409456/139‑1 (MQ=255)
                                                                                                                          cGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACt  <  1:129942/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
CCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACT  >  NC_000913/1245468‑1245728

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