Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 94,785 A→G 6.9% T212T (ACA→ACG)  gloB → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073094,7850AG 6.9% 155.8 / 2.3 58T212T (ACA→ACG) gloBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (2/2);  total (31/27)
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.13e-01

TTCTACAGAAATTAGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTGGTATGAATT  >  CP000730/94647‑94924
                                                                                                                                          |                                                                                                                                           
ttCTACAGAAATTAGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAAt                                                                                                                                           <  1:107823/141‑1 (MQ=255)
         aaTTAGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAGGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGtagcga                                                                                                                                  >  2:683937/1‑136 (MQ=255)
         aaTTAGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATga                                                                                                                                  >  2:695354/1‑141 (MQ=255)
             aGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGACATACCTATATCTACAATAGGTCATGAGaaa                                                                                                                              >  2:566918/1‑141 (MQ=255)
               cgcgCAACAAATGTATGAGTCCGTTCAAAATATTAAAAAATTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAATCCCTTCTGGTAAAGCATGCGGTGCCATCCCTAAATCTAAAATAGGTTATGACAAAAt                                                                                                                            >  2:247826/1‑141 (MQ=255)
                     aCAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCGGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATtaataa                                                                                                                      <  1:60117/141‑1 (MQ=255)
                        aaTGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTg                                                                                                                   <  2:621459/141‑1 (MQ=255)
                        aaTGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTg                                                                                                                   <  2:358429/141‑1 (MQ=255)
                        aaTGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTCATGAGAAAATTAATAACTg                                                                                                                   <  1:485885/141‑1 (MQ=255)
                         aTGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTgg                                                                                                                  <  1:774436/141‑1 (MQ=255)
                         aTGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGCGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTgg                                                                                                                  <  1:418365/141‑1 (MQ=255)
                          tGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACGATAGGTTATGAGAAAATTAATAACTggg                                                                                                                 >  1:216066/1‑141 (MQ=255)
                          tGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACGATAGGTTATGAGAAAATTAATAACTggg                                                                                                                 >  1:791143/1‑141 (MQ=255)
                             aTGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGATATGAGAAAATTAATAACTGGGCAt                                                                                                              >  2:198017/1‑141 (MQ=255)
                                 gTCCGTTTAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTaa                                                                                                          <  2:591126/141‑1 (MQ=255)
                                 gTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCGATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTaa                                                                                                          <  2:597033/141‑1 (MQ=255)
                                   tcttgcgaactattatATAGTTACCAGACTATGTTCGAATCTAGCCGTGTCATGGTGCTGGAAGCCCTTGTGGTACAGCATTAGGGGCCATACCTATATCTACAATAGGTTATGAGAAAAGTAATAACGGGGCATTTAATg                                                                                                        <  2:25934/131‑1 (MQ=255)
                                    cGTTCAAAATATTAAAAAGTTACCCGACTAGGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGa                                                                                                       <  1:318909/141‑1 (MQ=255)
                                        cAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAAtt                                                                                                   >  1:662414/1‑141 (MQ=255)
                                           aaTATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGATAATTGATATCTGGGCATTTAATGTAATtgat                                                                                                >  1:43338/1‑141 (MQ=255)
                                           aaTATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCAATTAATGGAATtgat                                                                                                >  2:49169/1‑141 (MQ=255)
                                             tatTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATtgatga                                                                                              >  2:268371/1‑141 (MQ=255)
                                                tAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGAc                                                                                           <  2:745435/141‑1 (MQ=255)
                                                  aaaaTTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTa                                                                                         >  2:606557/1‑141 (MQ=255)
                                                   aaaTTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAACAGGTGATGAGAAAATTAATACCTGGGCATTTAATGAAATTGATGAGACTaa                                                                                        >  1:665240/1‑141 (MQ=255)
                                                         aCCAGACTTGGTTCAAATCTGGCCGGTTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTAt                                                                                  <  2:719042/141‑1 (MQ=255)
                                                           ccgACTAGGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCCTTAGGTGCCATACCTCTATCTACAAAAGGTTATAAGAAATTGAATAACTGGGCATTTACTGAAATTGATGAGACTAAATTTATTg                                                                                <  2:334145/139‑1 (MQ=255)
                                                           cAGACTATGTTCAAATCTGGCCGGGGCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGCGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTg                                                                                <  2:358172/141‑1 (MQ=255)
                                                               cTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAGTAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTAATTGAAta                                                                            >  2:771565/1‑140 (MQ=255)
                                                               cTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATc                                                                            >  1:560817/1‑141 (MQ=255)
                                                                  tgtgCAAACCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCAtt                                                                         <  1:360789/141‑1 (MQ=255)
                                                                  tGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCAtt                                                                         <  1:6288/141‑1 (MQ=255)
                                                                              gCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACatcaaatca                                                             >  1:670085/1‑141 (MQ=255)
                                                                                cgggcaaTGGTGCTGGAAGCCATTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAAGTAATAACTGGGCATTTAATGAAATGGATGAGACAAAATTTATTGACTAATTAACATCAAATGAAc                                                           >  2:138170/7‑141 (MQ=255)
                                                                                 gggTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAAGCAAcc                                                          <  2:731560/141‑1 (MQ=255)
                                                                                  ggTCCTGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGGGCCATACCTATATCTACAATAGGTTATGAGAAAGTTAATAACTGGGCATTTAATGAAATTGAAGAGACTAAATTTACTGAATAATTAACATCAAATCAAACa                                                         >  2:615776/1‑141 (MQ=255)
                                                                                  ggTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCa                                                         >  2:313340/1‑141 (MQ=255)
                                                                                            gtgtatgtctttGTGGCAAAGTGTTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACTAGCACCACCGc                                               <  2:114488/131‑1 (MQ=255)
                                                                                            cTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACGATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGc                                               <  2:791143/141‑1 (MQ=255)
                                                                                            cTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACGATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGc                                               <  2:216066/141‑1 (MQ=255)
                                                                                            cTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGc                                               >  1:540742/1‑141 (MQ=255)
                                                                                             tGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCCGCACCACCGca                                              >  2:709721/1‑141 (MQ=255)
                                                                                                         tgGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGcaca                                  <  1:606557/141‑1 (MQ=255)
                                                                                                              aaGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGa                             >  1:97244/1‑141 (MQ=255)
                                                                                                               aGCATTAGCTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGCCATTTAATCAAATTGATCAGACTAAATTTATTGAATCATTAACATCAAATCCACCAGCACCACCGCATCATTTTGCACAAATGaa                            <  2:390378/141‑1 (MQ=255)
                                                                                                                  aTTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGCAACa                         >  1:614088/1‑141 (MQ=255)
                                                                                                                    tAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAACACCACCGCATCATTTGGCACAAATGAAACAAg                       >  1:315728/1‑141 (MQ=255)
                                                                                                                         gCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAAt                  >  2:286023/1‑141 (MQ=255)
                                                                                                                             tACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGt              <  2:97244/141‑1 (MQ=255)
                                                                                                                               ccTATATCGACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAACGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGttt            <  1:783152/141‑1 (MQ=255)
                                                                                                                                cTATATCTACAATAGGTTATGAGCAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAGTTTATTGAATCATTAACATCAAATCAACCAGCACACGCGCATCATTTTGCAAAAATGAAACAAGTTACTCAGtgtg           >  2:22841/1‑141 (MQ=255)
                                                                                                                                cTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTg           >  1:174877/1‑141 (MQ=255)
                                                                                                                                cTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTg           >  1:417935/1‑141 (MQ=255)
                                                                                                                                 tataTCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTgg          >  2:584433/1‑141 (MQ=255)
                                                                                                                                   tatCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATAAGTTTGGTa        >  1:32293/1‑141 (MQ=255)
                                                                                                                                        aCAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTGGTATGAAt   <  2:659994/141‑1 (MQ=255)
                                                                                                                                        aCAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTCATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTGGTATGAAt   <  2:659696/141‑1 (MQ=255)
                                                                                                                                         cAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAACCAAGTTAATCAGTTTGGTATGAAtt  >  2:43682/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TTCTACAGAAATTAGCGCGAAACAAATGTATGAGTCCGTTCAAAATATTAAAAATTTACCAGACTATGTTCAAATCTGGCCGGGTCATGGTGCTGGAAGCCCTTGTGGTAAAGCATTAGGTGCCATACCTATATCTACAATAGGTTATGAGAAAATTAATAACTGGGCATTTAATGAAATTGATGAGACTAAATTTATTGAATCATTAACATCAAATCAACCAGCACCACCGCATCATTTTGCACAAATGAAACAAGTTAATCAGTTTGGTATGAATT  >  CP000730/94647‑94924

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