Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 772,568 C→A Q371K (CAG→AAG)  cydA → cytochrome bd‑I ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913772,5680CA97.0% 111.3 / ‑2.8 33Q371K (CAG→AAG) cydAcytochrome bd‑I ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (15/17);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

TCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAT  >  NC_000913/772435‑772704
                                                                                                                                     |                                                                                                                                        
tctGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTc                                                                                                                                    >  1:33276/1‑140 (MQ=255)
  tGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCaa                                                                                                                                  >  1:242750/1‑140 (MQ=255)
        aacaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGGCTGAAGCGAAGATTCAACAGGCa                                                                                                                            >  1:615135/1‑140 (MQ=255)
         acaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCaa                                                                                                                           <  1:215734/140‑1 (MQ=255)
         acaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCaa                                                                                                                           <  2:587856/140‑1 (MQ=255)
              tGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCaaa                                                                                                                      <  2:242750/140‑1 (MQ=255)
                   tCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTc                                                                                                                 >  1:110512/1‑140 (MQ=255)
                              cGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtg                                                                                                      <  1:149512/140‑1 (MQ=255)
                               gACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtgt                                                                                                     >  2:511417/1‑140 (MQ=255)
                               gACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtgt                                                                                                     >  2:370813/1‑140 (MQ=255)
                                 ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg                                                                                                   <  2:552013/140‑1 (MQ=255)
                                 ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg                                                                                                   <  1:238843/140‑1 (MQ=255)
                                 ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg                                                                                                   <  1:511417/140‑1 (MQ=255)
                                  cAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAgc                                                                                                  >  2:145258/1‑140 (MQ=255)
                                   aGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAgcg                                                                                                 <  2:364518/140‑1 (MQ=255)
                                               aCCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACttt                                                                                     <  2:110512/140‑1 (MQ=255)
                                                   gTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGt                                                                                 >  1:21631/1‑140 (MQ=255)
                                                                 agaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTa                                                                           <  1:439155/132‑1 (MQ=255)
                                                                 agaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTa                                                                           >  2:439155/1‑132 (MQ=255)
                                                                       aCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgtgt                                                             >  1:655004/1‑140 (MQ=255)
                                                                       aCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgtgt                                                             <  1:219677/140‑1 (MQ=255)
                                                                             gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgt                                                               >  1:372336/1‑132 (MQ=255)
                                                                             gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgt                                                               <  2:372336/132‑1 (MQ=255)
                                                                             gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTc                                                       >  1:415134/1‑140 (MQ=255)
                                                                               tACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCt                                                     <  2:415134/140‑1 (MQ=255)
                                                                                 cGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGc                                                   >  1:236976/1‑140 (MQ=255)
                                                                                        ctgctgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGc                                            >  1:461753/1‑140 (MQ=255)
                                                                                            tgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAatc                                        <  2:461753/140‑1 (MQ=255)
                                                                                                  gCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATcgcg                                  <  2:414156/140‑1 (MQ=255)
                                                                                                   cTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCgcgc                                 <  2:33276/140‑1 (MQ=255)
                                                                                                   cTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCgcgc                                 <  2:21631/140‑1 (MQ=255)
                                                                                                                  ggCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAgtg                  >  2:68031/1‑140 (MQ=255)
                                                                                                                                 aGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCa   <  2:317525/140‑1 (MQ=255)
                                                                                                                                  gcgaAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAt  <  1:318676/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
TCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAT  >  NC_000913/772435‑772704

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