Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 195,286 T→A 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,2860TA98.2% 199.0 / ‑6.2 57intergenic (‑97/‑82)USA300HOU_RS00900/USA300HOU_RS00905hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (24/31);  minor base C (1/0);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAGAACTAG  >  USA300TCH1516_ALE/195148‑195428
                                                                                                                                          |                                                                                                                                              
tttGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAata                                                                                                                                       >  1:125198/1‑148 (MQ=255)
     ttCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTACATATTAGAAAAAATATCTTa                                                                                                                                  >  1:780950/1‑148 (MQ=255)
     ttCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTa                                                                                                                                  >  1:657856/1‑148 (MQ=255)
      tCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTAACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTCAAAAAAATATCTTaa                                                                                                                                 >  1:761175/1‑148 (MQ=255)
                                   tATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGt                                                                                                    <  1:94223/148‑1 (MQ=255)
                                              aTAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTc                                                                                         <  1:1076762/148‑1 (MQ=255)
                                                aGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCaa                                                                                       <  1:444033/148‑1 (MQ=255)
                                                   ccccTCTATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCa                                                                                    <  1:53745/148‑1 (MQ=255)
                                                   ccccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATTTCCTTGTAATCCGCTAAGAATTATAGTAGTATTTTTTCAGCCa                                                                                    >  1:265971/1‑148 (MQ=255)
                                                   ccccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCa                                                                                    <  1:781484/148‑1 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGTTAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:939582/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:982449/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:682619/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:587188/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:519300/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                                   >  1:316168/1‑148 (MQ=255)
                                                    cccTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGAAAAGAATTATAGTAATATTTTTTCaccaat                                                                                   >  1:1106201/1‑143 (MQ=255)
                                                         tatatatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTta                                                                              <  1:440341/148‑1 (MQ=255)
                                                         tatatatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAACATCCTTGGAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTta                                                                              <  1:1243652/148‑1 (MQ=255)
                                                            atatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATag                                                                           <  1:1044728/148‑1 (MQ=255)
                                                              atGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATTCTCCAGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTTTTATagga                                                                         >  1:544000/1‑148 (MQ=255)
                                                              atGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATagga                                                                         <  1:202862/148‑1 (MQ=255)
                                                                 tttatttaCTTATCCCCTAACATTATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATATTAATATTTTTTCAACCATTGTTATAGGAGGt                                                                      <  1:803279/148‑1 (MQ=255)
                                                                  ttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc                                                                     <  1:936466/148‑1 (MQ=255)
                                                                  ttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc                                                                     <  1:178612/148‑1 (MQ=255)
                                                                  ttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTc                                                                     <  1:1257656/148‑1 (MQ=255)
                                                                   tatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTATTATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCt                                                                    <  1:835367/148‑1 (MQ=255)
                                                                    ttttaCTTATACCCTAACATTATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCtt                                                                   <  1:272895/147‑1 (MQ=255)
                                                                    ttttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCtt                                                                   <  1:696711/147‑1 (MQ=255)
                                                                    atttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCtt                                                                   <  1:481756/148‑1 (MQ=255)
                                                                      ttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATACGAATTATAGTAATATTTTTTGAACCATTGTTTTAGGAGGTCttat                                                                 >  1:165037/1‑148 (MQ=255)
                                                                       taCTTATACTCTTTCTTGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCttatt                                                                <  1:23814/148‑1 (MQ=255)
                                                                       taCTTATACACTAACTAGATTTTTATATTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTCGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCttatt                                                                <  1:129484/148‑1 (MQ=255)
                                                                        aCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCttatta                                                               <  1:619386/148‑1 (MQ=255)
                                                                          ttATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAAt                                                             >  1:230186/1‑148 (MQ=255)
                                                                              aCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACa                                                         <  1:1187773/148‑1 (MQ=255)
                                                                               cccTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACat                                                        >  1:647735/1‑148 (MQ=255)
                                                                                 cTAATATGATTTTTATAATCATTGAAAATATATTTTACAGAATTTTATCTAAATCTTAAAAAAAATATCCTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACatta                                                      <  1:572106/148‑1 (MQ=255)
                                                                                 cTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACatta                                                      <  1:942387/148‑1 (MQ=255)
                                                                                 cTAACACGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACatta                                                      <  1:576434/148‑1 (MQ=255)
                                                                                  tAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACattat                                                     >  1:759530/1‑148 (MQ=255)
                                                                                     cATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTAtttt                                                  >  1:704849/1‑148 (MQ=255)
                                                                                     cATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATACTTGTAATCAGATAAGAATTATAGTAATATTTTTTCTACCATTGTTTTAGGTGGTCTTATTAATGACATtttttt                                                  >  1:876446/1‑148 (MQ=255)
                                                                                        gATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATGATAGTAATATTTTTTCAACCAGTGTTATAGGAGGTCTTATTAATGACATTATTTtgat                                               >  1:941259/1‑148 (MQ=255)
                                                                                          tttttATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttatta                                             >  1:789487/1‑148 (MQ=255)
                                                                                          tttttATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttatta                                             >  1:910132/1‑148 (MQ=255)
                                                                                                 tatctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTa                                      <  1:1116589/146‑1 (MQ=255)
                                                                                                       tGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAAtt                                >  1:159769/1‑147 (MQ=255)
                                                                                                          aaaTATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTg                             <  1:141932/148‑1 (MQ=255)
                                                                                                            atatatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAt                           <  1:8842/148‑1 (MQ=255)
                                                                                                              atatTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGAttt                         <  1:891078/148‑1 (MQ=255)
                                                                                                                    tACAGAATTTTATCTAAATATTGAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTATTGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATc                   <  1:869779/148‑1 (MQ=255)
                                                                                                                    tACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATc                   <  1:466080/148‑1 (MQ=255)
                                                                                                                         aaTTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa              >  1:890969/1‑148 (MQ=255)
                                                                                                                            tttATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACATTCTTGATTTTGCATCAACGAaag           <  1:621685/148‑1 (MQ=255)
                                                                                                                             ttATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAaaga          <  1:514283/148‑1 (MQ=255)
                                                                                                                                     atatTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATTAACGAAAGAAGAACTAg  >  1:828650/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TTTGTTTCTAGCTAAACTTGCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAGAACTAG  >  USA300TCH1516_ALE/195148‑195428

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