Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 76,969 G→A intergenic (+509/‑107) USA300HOU_RS00140 → / → USA300HOU_RS00350 hypothetical protein/Monoacylglycerol lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE76,9690GA98.4% 194.2 / NA 62intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (44/17);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAA  >  USA300TCH1516_ALE/76854‑77105
                                                                                                                   |                                                                                                                                        
cATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCa                                                                                                                 >  1:636872/1‑141 (MQ=255)
cATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCa                                                                                                                 >  1:627593/1‑141 (MQ=255)
 aTGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCaa                                                                                                                >  1:434825/1‑141 (MQ=255)
      aCGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAtt                                                                                                           <  1:129466/141‑1 (MQ=255)
      aCGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAtt                                                                                                           <  1:705403/141‑1 (MQ=255)
       cGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGTGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTg                                                                                                          <  1:405023/141‑1 (MQ=255)
                        tAGCAAGTTAAGCGAACACTGACAGGATAAATTAGTGGTTAGCTATATTTGTTTACTTTGCAACAGAACCTAGAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTAATGAAAACATTGATTTGGGTCAAAAAATa                                                                                         >  1:289823/1‑141 (MQ=25)
                            aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTCTTACTGAAAACATTGATTTGGCTCAAAAATTAgagg                                                                                     >  2:51691/1‑141 (MQ=255)
                            aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgggg                                                                                     >  2:895124/1‑141 (MQ=255)
                                  aGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAACCATTAAATTGGTACAAAAATTAGGGCTATTaa                                                                               >  2:157056/1‑141 (MQ=25)
                                  aGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaa                                                                               >  2:108601/1‑141 (MQ=255)
                                    cGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaa                                                                             >  2:739820/1‑141 (MQ=255)
                                        cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTATGGGTGTGAAAAAaat                                                                         >  1:150329/1‑141 (MQ=255)
                                        cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaat                                                                         >  1:821491/1‑141 (MQ=255)
                                        cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaat                                                                         >  1:743640/1‑141 (MQ=255)
                                         acTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatg                                                                        >  1:162120/1‑141 (MQ=255)
                                           tGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaa                                                                      >  1:921576/1‑141 (MQ=255)
                                           tGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaa                                                                      >  1:448895/1‑141 (MQ=255)
                                            gACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTAAAAAATTAGGGGTATTAAAAAaatgaat                                                                     >  1:15560/1‑141 (MQ=255)
                                            gACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGGCAAAAATTAGGGGTATTAAAAAaatgaat                                                                     >  1:22196/1‑141 (MQ=255)
                                              cATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAatgaatga                                                                   >  1:748344/1‑141 (MQ=255)
                                                tGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAACGTTAATTTATATAGTATTACTCCAAACATTGATTTGGGTCAAGAATTAGGGGTAATAAAAAAAtgaatgaaa                                                                 >  2:137521/1‑140 (MQ=38)
                                                tGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACATATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTCGGGTCACAGACTAGGGGTATTAAAAAAatggatgagt                                                                 >  1:347506/1‑139 (MQ=37)
                                                   tAAATTAGTGGTTAGCTATATTTTTTTTCTTTTTTTCAGAACCTATAAAGAACAGGGGCTTAAAATGTTAGTTTATATAGTAGTACTGCGAACATTGATTTGGGTCGAAAATTAGGGGTATTAAAAAAATGAGTGAATCAt                                                              >  2:700749/1‑141 (MQ=255)
                                                      aTTAGTGGTTAGCTATATTTTTTTACTTTGCTACAGAACCTATAAAAAACATAGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGACAAAAATTAGGGGTATTAAAAAAATGAACGAATCCTAAt                                                           >  2:482049/1‑141 (MQ=38)
                                                      aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATATTATTACTGCAAAAATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGGATGAATCTTAAt                                                           >  2:474902/1‑141 (MQ=38)
                                                      aTTAGTGGTTAGCTATATTTTTTTACTTTGCAAAAGAACCTATAAAAAACATGGGCTTAAAAAGTTAATTTATATAGTATTACTGCAAAAATTGATATGGGTCAAAAATTAAGGGTATTAAAAAAATGAATGAATCATAAt                                                           >  2:628049/1‑141 (MQ=38)
                                                            ggTTATCTATATTATTTTACTTTGCAACAGAACCTATAAATAACATGAGCTAAAAAGGTTAATTTATATAGTATTACTGCAAACATTCATGTGTGACAAAAAATAGGAGCATTAAAAAAATGAATTAATCATAATTTTATc                                                     <  2:561045/141‑1 (MQ=255)
                                                            ggTTAGCTATATTTTTGTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATc                                                     <  2:294999/141‑1 (MQ=255)
                                                                 gttataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCt                                                <  1:371827/139‑1 (MQ=255)
                                                                 gCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCt                                                <  1:615771/141‑1 (MQ=255)
                                                                  cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGACTCATAATTTTATCAAGCTg                                               >  1:106542/1‑141 (MQ=255)
                                                                  cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTg                                               >  1:793100/1‑141 (MQ=255)
                                                                    atatTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTTGGTCAAAAATTAGGGGTATTTAAAAAATGAATGAATCATAATTTTATCAAACttgt                                             >  2:155841/1‑138 (MQ=255)
                                                                    atatTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGAt                                             >  1:786986/1‑141 (MQ=255)
                                                                       tttttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAAGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTGTCATGGTGAtgga                                          >  2:765570/1‑138 (MQ=255)
                                                                         tttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgag                                        >  1:328781/1‑141 (MQ=255)
                                                                         tttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgag                                        >  2:25481/1‑141 (MQ=255)
                                                                          ttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGgaga                                       <  2:162120/141‑1 (MQ=255)
                                                                           tttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGagag                                      <  1:18039/141‑1 (MQ=255)
                                                                            ttACTGTGCAACAGAACCGATAAAGAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAAAAAGAAAAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAgg                                     <  2:888310/141‑1 (MQ=255)
                                                                             tACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGt                                    <  1:704741/141‑1 (MQ=255)
                                                                               cTTTGCAACAGAACCTATAGAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                                  >  1:689243/1‑141 (MQ=255)
                                                                               cTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTa                                  >  2:742173/1‑141 (MQ=255)
                                                                                 ttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTaaa                                <  2:737768/141‑1 (MQ=255)
                                                                                            ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAACCATTGATTTGGGTCAACCATTTGGGGTATTAATAAAATGAATGAATTATAATTTTATTAAGCTGATTGGTTATGTTAAAATGCATtatat                     >  2:633935/1‑141 (MQ=255)
                                                                                            ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatat                     >  2:620107/1‑141 (MQ=255)
                                                                                            ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAGTTTATCAAGCTGATTGGAGAGGTTAAAATGAATtatat                     >  2:689464/1‑141 (MQ=255)
                                                                                             cTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatata                    >  2:291219/1‑141 (MQ=255)
                                                                                             cTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATtatata                    >  2:480420/1‑141 (MQ=255)
                                                                                              tataAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaa                   >  1:518621/1‑141 (MQ=255)
                                                                                                taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa                 <  1:739820/141‑1 (MQ=255)
                                                                                                aaaaaaaCATGGGCTCAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGAGTATAAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATaaaa                 <  2:777943/140‑1 (MQ=255)
                                                                                                   aaaaCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAAttt              >  2:71924/1‑141 (MQ=255)
                                                                                                       cATGGGCTTAAAAGGTTAATTTATTTTCAATAATTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTg          >  1:445846/1‑141 (MQ=255)
                                                                                                       cATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTg          >  1:663958/1‑141 (MQ=255)
                                                                                                        aTGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGa         >  2:620103/1‑141 (MQ=255)
                                                                                                          gggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGt       <  1:474742/141‑1 (MQ=255)
                                                                                                           ggCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTc      <  1:96193/141‑1 (MQ=255)
                                                                                                              ttAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCaaa   >  1:183369/1‑141 (MQ=255)
                                                                                                               tAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCaaaa  >  1:602578/1‑141 (MQ=255)
                                                                                                               tAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCaaaa  <  1:155841/141‑1 (MQ=255)
                                                                                                                   |                                                                                                                                        
CATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCTGATTGGAGAGGTTAAAATGCATTATATAAAATTTATTGAGTCAAAA  >  USA300TCH1516_ALE/76854‑77105

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