Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,095,871 C→T 8.1% V341V (GTC→GTT cydA → cytochrome d ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,095,8710CT 8.1% 141.1 / 2.7 49V341V (GTC→GTTcydAcytochrome d ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

AGTGTCAAAACTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGACAGCCTTGGAT  >  CP000730/1095740‑1096011
                                                                                                                                   |                                                                                                                                            
aGTGTCAAAACTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTa                                                                                                                                     >  2:129093/1‑141 (MQ=255)
 gtgtCAAAACTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTaa                                                                                                                                    >  1:72854/1‑141 (MQ=255)
     cAAAACTAAAGTGCAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATTGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTTTTAACTTTAATGtt                                                                                                                                <  1:181885/141‑1 (MQ=255)
      aaaaCTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGttt                                                                                                                               >  1:39177/1‑141 (MQ=255)
      aaaaCTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGttt                                                                                                                               >  2:449646/1‑141 (MQ=255)
      aaaaCTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGttt                                                                                                                               >  2:430551/1‑141 (MQ=255)
           tAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTaaaaa                                                                                                                          <  2:330962/141‑1 (MQ=255)
           tAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTaaaaa                                                                                                                          >  1:60152/1‑141 (MQ=255)
                 aCAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAg                                                                                                                    <  1:91305/141‑1 (MQ=255)
                  cAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGa                                                                                                                   <  1:390242/141‑1 (MQ=255)
                       gTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGta                                                                                                              <  2:277223/141‑2 (MQ=255)
                           aaTGGTTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCt                                                                                                          >  1:433692/1‑141 (MQ=255)
                                ttttCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGCATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCa                                                                                                     >  1:55501/1‑141 (MQ=255)
                                       aaaaGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATgg                                                                                              <  2:383682/141‑1 (MQ=255)
                                       aaaaGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATgg                                                                                              >  2:163557/1‑141 (MQ=255)
                                        aaaGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGa                                                                                             <  2:148576/141‑1 (MQ=255)
                                         aaGTTTACATCCGCCTATGATTGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTTTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGAt                                                                                            >  2:429236/1‑141 (MQ=255)
                                          aGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATg                                                                                           >  1:237290/1‑141 (MQ=255)
                                              tACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCttt                                                                                       >  2:448344/1‑141 (MQ=255)
                                               aCATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTa                                                                                      >  1:148055/1‑141 (MQ=255)
                                               aCATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTa                                                                                      >  1:223885/1‑141 (MQ=255)
                                                 aTCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATg                                                                                    <  2:5592/141‑1 (MQ=255)
                                                     gCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAat                                                                                <  1:396594/141‑1 (MQ=255)
                                                          tGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAAtattat                                                                           >  1:265459/1‑141 (MQ=255)
                                                              cGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAAc                                                                       >  2:286442/1‑141 (MQ=255)
                                                                   attattTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGcc                                                                  <  2:50922/141‑1 (MQ=255)
                                                                   attattTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGcc                                                                  <  2:377425/141‑1 (MQ=255)
                                                                    ttattTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGccc                                                                 >  1:249959/1‑141 (MQ=255)
                                                                    ttattTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGccc                                                                 <  2:442974/141‑1 (MQ=255)
                                                                            ttGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTTTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTACGGAATATTATTAACAGGACCAGCTTCaa                                                         >  2:89867/1‑141 (MQ=255)
                                                                            ttGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCaa                                                         >  2:389352/1‑141 (MQ=255)
                                                                              gATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATg                                                       <  2:74720/141‑1 (MQ=255)
                                                                                tttAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATgct                                                     >  1:345042/1‑141 (MQ=255)
                                                                                    aTGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTGCTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCAAGCTTCAATgctagct                                                 <  2:394995/141‑1 (MQ=255)
                                                                                    aTGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATgctagct                                                 >  1:118790/1‑141 (MQ=255)
                                                                                            gATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAAtt                                         <  1:419335/141‑1 (MQ=255)
                                                                                                ggAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGa                                     >  1:439050/1‑141 (MQ=255)
                                                                                                   aTTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATgg                                  <  2:193742/141‑1 (MQ=255)
                                                                                                     tttATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTTAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGtt                                >  1:376469/1‑141 (MQ=255)
                                                                                                          gTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTaa                           <  2:187837/141‑1 (MQ=255)
                                                                                                            tttGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACa                         <  1:286442/141‑1 (MQ=255)
                                                                                                             ttGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACag                        >  1:267038/1‑141 (MQ=255)
                                                                                                                    ttCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATggg                 >  1:93087/1‑141 (MQ=255)
                                                                                                                      cAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTGTTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTa               <  1:264162/141‑1 (MQ=255)
                                                                                                                         gtgtTTATGTTTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTACGGAATATTATTAACAGGACCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGAc            <  1:89867/141‑1 (MQ=255)
                                                                                                                         gtgtTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGAc            <  1:399287/141‑1 (MQ=255)
                                                                                                                         gtgtTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGAc            <  1:236487/141‑1 (MQ=255)
                                                                                                                         gtgtTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGAc            <  2:191185/141‑1 (MQ=255)
                                                                                                                             ttATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGACAGcc        >  2:307808/1‑141 (MQ=255)
                                                                                                                                   cTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGACAGCCTTGGAt  <  1:55306/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                            
AGTGTCAAAACTAAAGTACAAGGGTTAAATGATTTTCCAAAAAGTTTACATCCGCCTATGATCGTCCATTATTTCTTTGATTTAATGGTAACGATGGGAATTTTATGTTTTGTCATTTCAGGTGTTTATGTCTTAACTTTAATGTTTAAAAAGCTTAGAAAGTTTTCTACTCATAAATGGATGCTTTATGGAATATTATTAACAGGCCCAGCTTCAATGCTAGCTATAGAATTTGGATGGTTCTTAACAGAGATGGGTAGACAGCCTTGGAT  >  CP000730/1095740‑1096011

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