Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 477,248 C→G 14.5% intergenic (+182/‑665) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported477,2480CG14.5% 148.3 / 16.1 62intergenic (+182/‑665)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (22/31);  new base G (5/4);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATA  >  Exported/477180‑477312
                                                                    |                                                                
aaTATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCAc                                                                  <  1:2832282/69‑1 (MQ=255)
  tatatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCActct                                                               >  1:19652/1‑70 (MQ=255)
  tatatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                                                              <  1:2729503/71‑1 (MQ=255)
  tatatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                                                              <  1:5477406/71‑1 (MQ=255)
     atTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAg                                                           >  1:3436746/1‑71 (MQ=255)
        tgcatgggccatgtaggcctgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAAgccg                                                        <  1:4894872‑M1/17‑2 (MQ=255)
         gcatgggccatgtaggcctgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       <  1:2736229‑M1/18‑3 (MQ=255)
         gcatgggccatgtaggcctgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       <  1:3883591‑M1/18‑3 (MQ=255)
         gcatgggccatgtaggcctgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       <  1:4351900‑M1/18‑3 (MQ=255)
           gagaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCaa                                                            >  1:4870258/1‑64 (MQ=255)
           gagaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCaa                                                            >  1:2066866/1‑64 (MQ=255)
            agaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                                                    <  1:1778325/71‑1 (MQ=255)
            agaTGATGTAGACACTGGCAAGCGTAGCGCATAAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                                                    >  1:1263320/1‑71 (MQ=255)
             gatgatGTAGACACTGGCAAGCGTAGCTCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                                                    <  1:5219710/70‑1 (MQ=255)
                 atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                               <  1:1643478/71‑1 (MQ=255)
                 atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                               >  1:1575023/1‑71 (MQ=255)
                 atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                               <  1:570277/71‑1 (MQ=255)
                    tAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                 <  1:5336468/66‑1 (MQ=255)
                     taggcctgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCaagccgcg                                                      >  1:2086626‑M1/42‑57 (MQ=255)
                      gACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCtt                                          >  1:1953073/1‑71 (MQ=255)
                      gACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCtt                                          >  1:3663150/1‑71 (MQ=255)
                       acacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATACACGTGGCtt                                          <  1:2439347/70‑1 (MQ=25)
                         acTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGtata                                                  >  1:4348652/1‑60 (MQ=255)
                         acTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAt                                       >  1:3919259/1‑71 (MQ=255)
                          cTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                                                              <  1:2436253/47‑1 (MQ=255)
                          cTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                                                              <  1:5334242/47‑1 (MQ=255)
                          cTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                      >  1:1714366/1‑71 (MQ=255)
                           tgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       >  1:2247114‑M1/36‑51 (MQ=255)
                           tgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       >  1:5023272‑M1/36‑51 (MQ=255)
                           tgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       >  1:4223045‑M1/36‑51 (MQ=255)
                           tgataagcgtagcgcatcaggcaattttacctttgTCATCAGTCTCAagccgc                                                       >  1:3714582‑M1/36‑51 (MQ=255)
                            ggCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt                                    >  1:854474/1‑71 (MQ=255)
                            ggCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttt                                     >  1:3603881/1‑70 (MQ=255)
                             gCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCActct                                                               <  1:5147042/43‑1 (MQ=255)
                             gCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCActct                                                               <  1:4176340/43‑1 (MQ=255)
                              cAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt                                    <  1:3464564/69‑1 (MQ=255)
                                aGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                                                    <  1:2671208/51‑1 (MQ=255)
                                aGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGta                                                    <  1:5022141/51‑1 (MQ=255)
                                  cGTAGCGCATCATGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGtat                                                   >  1:4457864/1‑50 (MQ=255)
                                  cGTAGCGCATCAGGTAATTTTGCTTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaa                              <  1:5094143/71‑1 (MQ=255)
                                  cGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACaa                              <  1:948178/71‑1 (MQ=255)
                                      gcgcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCa                                                        >  1:2508926/1‑41 (MQ=255)
                                      gagcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCa                                                        >  1:2831984/3‑41 (MQ=255)
                                        gcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTgg                                             >  1:1031412/1‑50 (MQ=38)
                                             aGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAAccc                   <  1:3482073/71‑1 (MQ=255)
                                              ggTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                  <  1:5551177/71‑1 (MQ=255)
                                              ggTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                  <  1:1673090/71‑1 (MQ=255)
                                               gTAATTTTGCGTTTATCTTCACTCTCAAGCCACGt                                                     <  1:1164678/35‑1 (MQ=255)
                                               gTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                             <  1:3381039/43‑1 (MQ=255)
                                               gTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                             <  1:919856/43‑1 (MQ=255)
                                               gTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                             <  1:4398554/43‑1 (MQ=255)
                                                tAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTg                 <  1:1161026/70‑1 (MQ=255)
                                                 aaTTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                 >  1:2689727/1‑37 (MQ=255)
                                                 aaTTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                 >  1:4845500/1‑37 (MQ=255)
                                                   ttttGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGc                                            <  1:1206574/40‑1 (MQ=255)
                                                   ttttGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttttt                                   <  1:4490166/49‑1 (MQ=255)
                                                        cGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTa        <  1:1595869/71‑1 (MQ=255)
                                                         gTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTa                                        <  1:3191139/38‑1 (MQ=255)
                                                         gTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                      <  1:5643857/40‑1 (MQ=255)
                                                          tttATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaataa                         <  1:2106858/52‑1 (MQ=255)
                                                          tttATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaat                           <  1:3310780/50‑1 (MQ=255)
                                                           ttATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt                                    >  1:1753974/1‑40 (MQ=255)
                                                           ttATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAt      >  1:3210392/1‑70 (MQ=255)
                                                              tCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtat   >  1:2150869/1‑70 (MQ=255)
                                                                  cACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtata  <  1:1507974/67‑1 (MQ=255)
                                                                   aCTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                >  1:5387776/1‑52 (MQ=255)
                                                                    |                                                                
AATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATA  >  Exported/477180‑477312

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