Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 4,563,812 C→T S693S (AGC→AGT CJ019_04322 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,563,8120CT96.2% 166.8 / ‑6.6 53S693S (AGC→AGTCJ019_04322hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (21/29);  minor base A (1/1);  total (22/31)
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

CAGACCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAAT  >  ACB143/4563681‑4563944
                                                                                                                                   |                                                                                                                                    
cAGACCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAgtgcgc                                                                                                                                  >  5:307399/1‑136 (MQ=255)
cAGACCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAgtgcgc                                                                                                                                  >  6:53590/1‑136 (MQ=255)
   tcccTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGa                                                                                                                               <  1:178214/135‑1 (MQ=255)
   aCCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGTGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGa                                                                                                                               <  2:104715/136‑1 (MQ=255)
   aCCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGa                                                                                                                               <  2:153506/136‑1 (MQ=255)
          tCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATTCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGc                                                                                                                        >  3:70569/1‑136 (MQ=255)
          tCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGCACGAGGc                                                                                                                        >  2:285715/1‑136 (MQ=255)
          tCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGc                                                                                                                        >  3:199965/1‑136 (MQ=255)
                      ccTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGgcgc                                                                                                            <  1:122858/136‑1 (MQ=255)
                      ccTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGgcgc                                                                                                            <  5:345416/136‑1 (MQ=255)
                          ggtccTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCt                                                                                                          <  7:155417/131‑1 (MQ=255)
                                   ccTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTAc                                                                                               >  2:148163/1‑136 (MQ=255)
                                   ccTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTAc                                                                                               >  5:141452/1‑136 (MQ=255)
                                    cTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTAc                                                                                               >  4:23529/1‑135 (MQ=255)
                                     tGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTAc                                                                                                                    <  7:188070/113‑1 (MQ=255)
                                     tGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTAc                                                                                                                    >  8:188070/1‑113 (MQ=255)
                                       aaCGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCt                                                                                             <  5:173602/134‑1 (MQ=255)
                                          gactacAATTCCCTTCTTCAATAATCCGAGCGTGTTCGGCAGTGGAAGCAGAAGGCCCTGCGCGAGGCAAATATTCGCAGTAGATGTAGTGCGCCGAAAGAGGCCTAGGAAGGCGACTGATCTCATTACCTCGAct                                                                                        <  4:235289/136‑2 (MQ=255)
                                          gacgacAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        <  7:133010/136‑1 (MQ=255)
                                          gacgacAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        <  6:153630/136‑1 (MQ=255)
                                          gacgacAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        <  1:359303/136‑1 (MQ=255)
                                          gacgacAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGAGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        <  2:227499/136‑1 (MQ=255)
                                            cgacAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGc                                                                                      <  6:67752/136‑1 (MQ=255)
                                                    cccTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCctgctgct                                                                              <  4:199965/136‑1 (MQ=255)
                                                    cccTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCctgctgct                                                                              <  2:15114/136‑1 (MQ=255)
                                                    cccTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCctgctgct                                                                              >  3:22488/1‑136 (MQ=255)
                                                    cccTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCctgctgct                                                                              <  3:236055/136‑1 (MQ=255)
                                                          gTCAATACGCCGAGCGTGTTCGCCAGTTTCAGAAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAg                                                                        >  2:200002/1‑136 (MQ=255)
                                                          gTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGTGTGCGCCGAACGAGGCCTACGACGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAg                                                                        <  4:14859/136‑1 (MQ=255)
                                                            cAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCt                                                                                                     >  8:342087/1‑105 (MQ=255)
                                                            cAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCt                                                                                                     <  7:342087/105‑1 (MQ=255)
                                                            cAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGAGCGCAGAACGAGGCCTACGAAGGCGCCTGCCCTCATTACCTCGACGGCCTGCTGCTGGACAGc                                                                       >  6:60471/1‑135 (MQ=255)
                                                                      aGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGc                                                             >  2:245189/1‑135 (MQ=255)
                                                                         gtgtTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACt                                                         >  7:173401/1‑136 (MQ=255)
                                                                         gtgtTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACt                                                         >  2:207053/1‑136 (MQ=255)
                                                                           gtTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTgc                                                       <  3:350400/136‑1 (MQ=255)
                                                                                ccAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAg                                                  <  6:51477/136‑1 (MQ=255)
                                                                                       cagcagAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTgg                                           >  5:177376/1‑136 (MQ=255)
                                                                                          cagAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGTCAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCg                                        <  2:219458/136‑1 (MQ=255)
                                                                                          cagAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCg                                        <  1:114649/136‑1 (MQ=255)
                                                                                                  ccTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGCCAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACgcggcg                                >  3:325864/1‑136 (MQ=255)
                                                                                                       gcgAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgct                           >  3:172522/1‑136 (MQ=255)
                                                                                                       gcgAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgct                           >  8:207167/1‑136 (MQ=255)
                                                                                                       gcgAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgct                           >  4:46568/1‑136 (MQ=255)
                                                                                                       gcgAGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgct                           >  3:344611/1‑136 (MQ=255)
                                                                                                          aGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgctgct                        <  3:326404/136‑1 (MQ=255)
                                                                                                          aGGCCAAGCTGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgctgct                        <  6:315081/136‑1 (MQ=255)
                                                                                                          aGGCCAAGCGGCGCAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAgctgct                        <  7:327146/136‑1 (MQ=255)
                                                                                                                     cgcAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCg             <  1:200002/136‑1 (MQ=255)
                                                                                                                      gcAGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCgg            <  8:189471/136‑1 (MQ=255)
                                                                                                                        aGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        <  4:237889/58‑1 (MQ=255)
                                                                                                                        aGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACg                                                                                        >  3:237889/1‑58 (MQ=255)
                                                                                                                        aGTAGCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGgcg          >  3:68364/1‑136 (MQ=255)
                                                                                                                            gCTGGAGTGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGTTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGctct      <  8:265168/136‑1 (MQ=255)
                                                                                                                                gAGTGCGCCGAACGAGGCCTACGAAGGCGGCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAAt  <  7:262658/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CAGACCCTGCTCGGCAGTTGGCCCTTGGACCTCGACCTGAACGACGACAATGCCCTGCGTCAATACGCCGAGCGTGTTCGCCAGTGGCAGCAGAAGGCCCTGCGCGAGGCCAAGCTGCGCAGTAGCTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAAT  >  ACB143/4563681‑4563944

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