Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 195,286 T→A 100% intergenic (‑97/‑82) USA300HOU_RS00900 ← / → USA300HOU_RS00905 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE195,2860TA89.7% 185.2 / 14.2 59intergenic (‑97/‑82)USA300HOU_RS00900/USA300HOU_RS00905hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base A (32/20);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAA  >  USA300TCH1516_ALE/195146‑195421
                                                                                                                                            |                                                                                                                                       
tttttGTTTCTAGCTAAACTTGCTACGACAGCTTCTTTATCCATAATCATCGCACCATATATATATGTTTATTTATTTATACCCTATCATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATAttt                                                                                                                                         <  2:114046/141‑1 (MQ=255)
  tttGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATATCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTaaa                                                                                                                                       >  1:95373/1‑141 (MQ=255)
    tGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAaaaa                                                                                                                                     >  1:192362/1‑141 (MQ=255)
         cTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAata                                                                                                                                >  2:97716/1‑141 (MQ=255)
             cTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCtt                                                                                                                            <  1:593042/141‑1 (MQ=255)
                  cTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAatat                                                                                                                       >  1:687694/1‑141 (MQ=255)
                  cTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAatat                                                                                                                       >  1:233879/1‑141 (MQ=255)
                    tGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATcc                                                                                                                     >  1:774430/1‑141 (MQ=255)
                    tGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATcc                                                                                                                     >  1:515240/1‑141 (MQ=255)
                           aCAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAAt                                                                                                              <  1:108120/141‑1 (MQ=255)
                           aCAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAAt                                                                                                              >  2:320066/1‑141 (MQ=255)
                           aCAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAAt                                                                                                              >  2:109137/1‑141 (MQ=255)
                                  cTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAACTATAATTAAAAGAATTTTAAATAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATaa                                                                                                       >  1:115578/1‑141 (MQ=255)
                                   tttATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAg                                                                                                      >  1:391286/1‑141 (MQ=255)
                                    ttATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAATATATTTTACAGAATTTTATTTAAATATTAAAAAAAAAATTTTAATATTCTTTTAATTCGGTAAGaa                                                                                                    >  2:235027/1‑141 (MQ=255)
                                     tATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTGTGGAAAGATATTTTACAGAATTTTATCTGAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGaa                                                                                                    >  1:470853/1‑141 (MQ=255)
                                        ccATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATta                                                                                                 >  2:103763/1‑141 (MQ=255)
                                            aaTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGt                                                                                             <  2:69219/141‑1 (MQ=255)
                                            aaTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGt                                                                                             >  1:476936/1‑141 (MQ=255)
                                            aaTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGt                                                                                             >  1:175845/1‑141 (MQ=255)
                                             aTGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTa                                                                                            >  1:183560/1‑141 (MQ=255)
                                                aTAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAata                                                                                         <  1:39050/141‑1 (MQ=255)
                                                  aGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATAtt                                                                                       >  1:120479/1‑141 (MQ=255)
                                                   gCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATAttt                                                                                      >  1:769919/1‑141 (MQ=255)
                                                            atatatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCa                                                                             >  1:355655/1‑141 (MQ=255)
                                                                   tttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTtat                                                                      <  2:154534/141‑1 (MQ=255)
                                                                     tatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATag                                                                    >  1:409761/1‑141 (MQ=255)
                                                                     tatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATag                                                                    <  1:61648/141‑1 (MQ=255)
                                                                       tttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATagga                                                                  >  1:257499/1‑141 (MQ=255)
                                                                         taCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATaggagg                                                                <  2:211854/141‑1 (MQ=255)
                                                                          aCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAGTCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGt                                                               <  1:722092/141‑1 (MQ=255)
                                                                          aCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGt                                                               >  1:207650/1‑141 (MQ=255)
                                                                           cTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc                                                              >  2:176181/1‑141 (MQ=255)
                                                                           cTTATACCCGAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc                                                              >  1:691040/1‑141 (MQ=255)
                                                                             tataCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCtt                                                            >  1:778560/1‑141 (MQ=255)
                                                                              ataCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCGACCATTGTTATAGGAGGTCtta                                                           >  2:278828/1‑141 (MQ=255)
                                                                               taCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCttat                                                          >  2:221096/1‑141 (MQ=255)
                                                                                        aTGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACat                                                 <  2:32555/141‑1 (MQ=255)
                                                                                          gATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACatta                                               <  1:703491/141‑1 (MQ=255)
                                                                                             ttttATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTAttt                                            >  1:433124/1‑141 (MQ=255)
                                                                                               ttATACTCTTCGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTAttttt                                          <  1:721319/141‑1 (MQ=255)
                                                                                                tataCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTtta                                         <  2:271557/141‑1 (MQ=255)
                                                                                                 ataCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttat                                        >  1:24294/1‑141 (MQ=255)
                                                                                                   aCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttatta                                      <  2:232347/141‑1 (MQ=255)
                                                                                                      ctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGaa                                   >  2:335125/1‑141 (MQ=255)
                                                                                                          gAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTa                               <  2:333021/141‑1 (MQ=255)
                                                                                                            aaaTATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAAc                             >  1:601157/1‑141 (MQ=255)
                                                                                                              atatatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACaa                           >  1:322707/1‑141 (MQ=255)
                                                                                                              atatatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACaa                           >  1:412954/1‑141 (MQ=255)
                                                                                                                 tatTTTACAGAATTTTATCTAAATATTCAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCt                        <  2:222124/141‑1 (MQ=255)
                                                                                                                 tatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCt                        <  1:640386/141‑1 (MQ=255)
                                                                                                                    tttACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATACTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGa                     >  2:110057/1‑141 (MQ=255)
                                                                                                                          gAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGc               <  1:178516/141‑1 (MQ=255)
                                                                                                                            aTTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCAt             <  1:223986/141‑1 (MQ=255)
                                                                                                                               ttATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCaa          >  1:178410/1‑141 (MQ=255)
                                                                                                                                  tCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa       >  2:40066/1‑141 (MQ=255)
                                                                                                                                  tCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa       >  1:461637/1‑141 (MQ=255)
                                                                                                                                  tCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa       <  1:181285/141‑1 (MQ=255)
                                                                                                                                    tAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGaaa     >  2:245761/1‑141 (MQ=255)
                                                                                                                                     aaaTATTAAAAAAAATATCTTAATATCCTTGTACTCCGATAAGAATTATAGTAATATTATTTCAACCATTGTTATAGGAGGGCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAaag    <  2:264997/141‑1 (MQ=255)
                                                                                                                                       atatTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAaagaa  <  2:73286/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TTTTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAA  >  USA300TCH1516_ALE/195146‑195421

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