Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,181,670 T→A 6.9% K19I (AAA→ATA)  U712_10920 ← Menaquinol‑cytochrome c reductase cytochrome b/c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,181,6700TA 6.9% 183.4 / 4.5 58K19I (AAA→ATA) U712_10920Menaquinol‑cytochrome c reductase cytochrome b/c subunit
Reads supporting (aligned to +/- strand):  ref base T (19/35);  new base A (2/2);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCTGATGCACAATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATG  >  CP006881/2181533‑2181804
                                                                                                                                         |                                                                                                                                      
gcTGATGCACAATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCttttt                                                                                                                                       >  2:209327/1‑139 (MQ=255)
gcTGATGCACAATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCttttt                                                                                                                                       >  1:89386/1‑139 (MQ=255)
  tGATGCACAATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCttttttt                                                                                                                                     >  7:98559/1‑139 (MQ=255)
       cacaATCGTCAAAACAAGATAACAGACCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCAAGAAAGCCTCTGTTTTACCGGGATAGTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGc                                                                                                                                 <  5:92909/138‑1 (MQ=255)
       cacaATCGTCAAAACAACAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGc                                                                                                                                 <  5:150555/138‑1 (MQ=255)
        acaATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCg                                                                                                                               <  2:269587/139‑1 (MQ=255)
           aTCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAACGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTa                                                                                                                            >  5:129891/1‑139 (MQ=255)
               tCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCt                                                                                                                        >  4:148703/1‑139 (MQ=255)
                    actaGAAAACCGATCCCAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGt                                                                                                                   <  3:366228/136‑1 (MQ=255)
                     cAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTc                                                                                                                  <  3:287778/139‑1 (MQ=255)
                            aCCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACa                                                                                                           <  4:94955/139‑1 (MQ=255)
                            aCCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACa                                                                                                           <  4:130266/139‑1 (MQ=255)
                               gATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAAt                                                                                                        >  5:132884/1‑139 (MQ=255)
                                aTCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTc                                                                                                                                    <  6:133402/110‑1 (MQ=255)
                                aTCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTc                                                                                                                                    >  5:133402/1‑110 (MQ=255)
                                aTCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAAtt                                                                                                       >  1:360435/1‑139 (MQ=255)
                                   acaacgACAGCTCCAACCATCCCTTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCa                                                                                                    <  1:159967/134‑1 (MQ=255)
                                   aGAAAGACTGGTCCAACCATCCATTCCTTCAACAAGTAATTCGGCGAGTAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTTTATTCGGCTTTTTTTCAGCCGGTATGCTTGAGTCGCCTACTAATTTCa                                                                                                    <  2:71485/139‑1 (MQ=255)
                                   aGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCa                                                                                                    <  5:348740/139‑1 (MQ=255)
                                        gACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCt                                                                                               <  8:282736/139‑1 (MQ=255)
                                         aCAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGTGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCtt                                                                                              <  8:98559/139‑1 (MQ=255)
                                            gCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGc                                                                                           >  2:76040/1‑139 (MQ=255)
                                            gCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGc                                                                                           >  6:166156/1‑139 (MQ=255)
                                             cTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATATCATCCCATTGcc                                                                                          >  7:266947/1‑139 (MQ=255)
                                                cAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         <  7:91799/137‑1 (MQ=255)
                                                cAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         <  5:7828/137‑1 (MQ=255)
                                                cAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         <  5:166156/137‑1 (MQ=255)
                                                cAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         <  1:76040/137‑1 (MQ=255)
                                                  aCCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGt                                                                                     <  6:27044/139‑1 (MQ=255)
                                                        catTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGt                                                                               <  8:30868/139‑1 (MQ=255)
                                                           tCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCAt                                                                            <  8:175122/139‑1 (MQ=255)
                                                               tCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGc                                                                                                                                 >  6:74490/1‑82 (MQ=255)
                                                               tCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGc                                                                                                                                 <  5:74490/82‑1 (MQ=255)
                                                               tCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATcccc                                                                        >  1:130312/1‑139 (MQ=255)
                                                                caacaaGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATccc                                                                         <  1:209327/137‑1 (MQ=255)
                                                                 aacaaGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCa                                                                                                    <  8:295019/109‑1 (MQ=255)
                                                                 aacaaGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCa                                                                                                    >  7:295019/1‑109 (MQ=255)
                                                                  acaaGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCccctcc                                                                     >  6:283616/1‑139 (MQ=255)
                                                                             cgCCAGAAAGCCTCTGGTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCg                                                           <  4:266505/137‑1 (MQ=255)
                                                                                    aaGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCaa                                                   <  4:160628/139‑1 (MQ=255)
                                                                                               ttACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         <  4:120222/138‑1 (MQ=255)
                                                                                               gacaggggATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         <  5:227537/133‑1 (MQ=255)
                                                                                                 aCCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCg                                       <  1:208142/138‑1 (MQ=255)
                                                                                                               aataatCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCt                        >  8:4089/1‑139 (MQ=255)
                                                                                                               aataatCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCAAAAAATTTAAACCCTTTGCCCCGGTGCATTTCATCCCCTCCCATTTAATAGCACATCAAAAGTTTTAAGCGGTCCTGaaaatgtct                        >  6:227537/1‑133 (MQ=255)
                                                                                                                ataatCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCaa                                                   >  8:143296/1‑111 (MQ=255)
                                                                                                                ataatCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCaa                                                   <  7:143296/111‑1 (MQ=255)
                                                                                                                   atCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCtt                    <  2:360435/139‑1 (MQ=255)
                                                                                                                    tCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGGATGTCATCCCCTCCCTTTTATTCGGAGCTGAAGAGTTTTAAGGGGGCCGGAAATCCCCTGCTTc                   >  3:23357/1‑139 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         <  6:189786/64‑1 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGccc                                                                                         >  5:189786/1‑64 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGa               <  2:96715/138‑1 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGa               >  1:157675/1‑138 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTATTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         >  5:196918/1‑112 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTATTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         <  6:201771/112‑1 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTATTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         <  6:196918/112‑1 (MQ=255)
                                                                                                                         cGGTATATTCGGCTTTATTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAg                                         >  5:201771/1‑112 (MQ=255)
                                                                                                                          ggTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGATc             <  2:130312/139‑1 (MQ=255)
                                                                                                                          ggTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGATc             <  3:148703/139‑1 (MQ=255)
                                                                                                                                    gCTTTTTTTCAGCCGGTATCCTTGAGGCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAAt   <  1:3801/139‑1 (MQ=255)
                                                                                                                                     cTTTTTTTCAGCCGGTATCCTTGAGTCGCGGACAAATTTGATCCCTTGGCCCCGGGGCATGTGCTGCCCGCCCAGTTTTTCGCAGAGAAAGAGAGTTAAGCGGTGCTGAAATCCCCGGCGTCCGGTTCATAATAAAAtg  >  2:361159/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                      
GCTGATGCACAATCGTCAAAACAAGAAAACCGATCAGAAAGACAGCTCCAACCATCCATTCCTTCAACAAGAAATTCGGCCAGAAAGCCTCTGTTTTACCGGGATACTCCGAATAATCCTTCGGTATATTCGGCTTTTTTTCAGCCGGTATCCTTGAGTCGCCTACAAATTTCATCCCTTTGCCCCGGTGCATGTCATCCCCTCCCTTTTATTCGCAGATCAAGAGTTTAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATG  >  CP006881/2181533‑2181804

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