Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,478,130 T→A 6.6% intergenic (+17/‑558) USA300HOU_RS12760 → / → USA300HOU_RS12765 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,478,1300TA 6.6% 714.7 / 28.2 214intergenic (+17/‑558)USA300HOU_RS12760/USA300HOU_RS12765hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (111/88);  new base A (3/11);  total (114/100)
Fisher's exact test for biased strand distribution p-value = 2.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 232 total aligned reads displayed.

GATTGTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAAC  >  USA300TCH1516_ALE/2478040‑2478220
                                                                                          |                                                                                          
    gTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattatta                                                                                        >  2:1431071/1‑91 (MQ=255)
     ttattaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattattat                                                                                       <  1:3415710/91‑1 (MQ=255)
     ttattaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattattat                                                                                       <  1:2538832/91‑1 (MQ=255)
      tattaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAAGAAACTATTATtata                                                                                      <  2:3296557/91‑1 (MQ=255)
      tattaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATtata                                                                                      <  1:3250215/91‑1 (MQ=255)
       attaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAg                                                                                     <  1:1534613/91‑1 (MQ=255)
       attaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAg                                                                                     <  2:3796584/91‑1 (MQ=255)
        ttaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAgt                                                                                    <  1:1231327/91‑2 (MQ=255)
        ttaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAgg                                                                                    >  1:2418553/1‑91 (MQ=255)
         taACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgt                                                                                   >  2:3372772/1‑91 (MQ=255)
          aaCGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgtt                                                                                  >  1:3160649/1‑91 (MQ=255)
            cGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAGtttgt                                                                                <  2:1272279/91‑1 (MQ=38)
            cGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgttgt                                                                                >  1:1278930/1‑91 (MQ=255)
               tataTTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTc                                                                             >  1:3697310/1‑91 (MQ=255)
               tataTTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTc                                                                             <  2:456939/91‑1 (MQ=255)
                atatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCt                                                                            <  2:2930445/91‑1 (MQ=255)
                 tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                           >  1:1960169/1‑91 (MQ=255)
                 tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                           <  2:3047784/91‑1 (MQ=255)
                 tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                           <  1:2574434/91‑1 (MQ=255)
                 tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                           <  1:1796350/91‑1 (MQ=255)
                  atTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTaa                                                                          <  2:41445/91‑1 (MQ=255)
                   ttGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAat                                                                         >  2:1047073/1‑91 (MQ=255)
                    tGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTACTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAata                                                                        >  2:3193504/1‑91 (MQ=255)
                     gTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAatat                                                                       <  2:658023/91‑1 (MQ=255)
                      ttCAAGTAATTGGAATAAATGTCTTTACGCCGCTATTATTTATAATAGTTAATACACAAACAAACTATTAATATAGGAAGTTTataataat                                                                      >  2:3627930/1‑89 (MQ=255)
                            tagtTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCa                                                                >  1:1047624/1‑91 (MQ=255)
                              gtTGGTTTATATGTCTTTACGCCGCTTTTGTTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCaaa                                                              <  2:741657/91‑1 (MQ=255)
                               ttGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAAc                                                             >  1:2790581/1‑91 (MQ=255)
                                tGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACt                                                            >  1:820337/1‑91 (MQ=255)
                                 ggTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTa                                                           >  1:2386986/1‑91 (MQ=255)
                                    ttATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCt                                                        >  2:2062428/1‑91 (MQ=255)
                                    ttATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCt                                                        >  1:2711195/1‑91 (MQ=255)
                                    ttATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTCGTTTCTAATATTACGTCAAACTATCt                                                        <  2:3218571/91‑1 (MQ=255)
                                      atatGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCttt                                                      <  1:1674612/91‑1 (MQ=255)
                                          gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                  >  1:3665624/1‑91 (MQ=255)
                                          gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                  >  2:499275/1‑91 (MQ=255)
                                          gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                  <  1:1598046/91‑1 (MQ=255)
                                          gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                  >  1:2996651/1‑91 (MQ=255)
                                           tCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAata                                                 <  2:1201238/91‑1 (MQ=255)
                                           tCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAata                                                 >  1:1690986/1‑91 (MQ=255)
                                            cTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTACatat                                                >  2:3343003/1‑91 (MQ=255)
                                             tttACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATCTTGCGTCAAACTATCTTTAAATAtt                                               >  1:2655248/1‑91 (MQ=255)
                                             tttACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATAtt                                               >  1:2305163/1‑91 (MQ=255)
                                              ttACGCCGCGTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTc                                              <  2:476888/91‑1 (MQ=255)
                                               tACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCg                                             <  1:2056433/91‑1 (MQ=255)
                                                aCGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGt                                            <  1:3024516/91‑1 (MQ=255)
                                                 cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                           >  2:191473/1‑91 (MQ=255)
                                                 cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                           <  2:859353/91‑1 (MQ=255)
                                                 cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                           >  1:781084/1‑91 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAGATATTCGTTaa                                         >  1:3516330/1‑91 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         <  2:2911354/91‑1 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         >  2:567341/1‑91 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         >  2:2755114/1‑91 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         <  1:962045/91‑1 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         >  1:1442233/1‑91 (MQ=255)
                                                   ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                         >  2:2007755/1‑91 (MQ=255)
                                                    cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                        <  1:113625/91‑1 (MQ=255)
                                                    cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                        <  2:3462709/91‑1 (MQ=255)
                                                    cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                        <  2:153652/91‑1 (MQ=255)
                                                      cTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATcc                                      <  1:911061/91‑1 (MQ=255)
                                                          tatttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGa                                  <  2:2594246/91‑1 (MQ=255)
                                                           atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGAt                                 <  1:1031578/91‑1 (MQ=255)
                                                           atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGAt                                 >  2:1426632/1‑91 (MQ=255)
                                                           atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATAGTCGTTAATCCTTGAt                                 <  1:1033028/91‑1 (MQ=255)
                                                           atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTGGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGAt                                 >  2:2943979/1‑91 (MQ=255)
                                                             ttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTa                               >  1:237024/1‑91 (MQ=255)
                                                             ttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTa                               <  2:2861591/91‑1 (MQ=255)
                                                              tataATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAg                              <  2:288255/91‑1 (MQ=255)
                                                              tataATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAg                              >  1:3148703/1‑91 (MQ=255)
                                                              tataATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAg                              <  1:2829384/91‑1 (MQ=255)
                                                               ataataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGa                             <  2:887305/91‑1 (MQ=255)
                                                                 aataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACa                           <  1:834549/91‑1 (MQ=255)
                                                                  ataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACaa                          <  1:2740552/91‑1 (MQ=255)
                                                                  ataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACaa                          <  1:356233/91‑1 (MQ=255)
                                                                   taGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAAt                         <  2:2064650/91‑1 (MQ=255)
                                                                    aGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAAtt                        <  1:2007755/91‑1 (MQ=255)
                                                                    aGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAAtt                        >  2:1572672/1‑91 (MQ=255)
                                                                     gTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTg                       <  2:2655271/91‑1 (MQ=255)
                                                                     gTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTg                       <  2:1724607/91‑1 (MQ=255)
                                                                     gTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTg                       <  1:700773/91‑1 (MQ=255)
                                                                      ttAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGt                      >  1:1711483/1‑91 (MQ=255)
                                                                      ttAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGt                      >  2:3104689/1‑91 (MQ=255)
                                                                       tAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGtt                     >  1:3832831/1‑91 (MQ=255)
                                                                        aaTCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTta                    <  1:101664/91‑1 (MQ=255)
                                                                         aTCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTtat                   <  2:1762321/91‑1 (MQ=255)
                                                                         aTCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTtat                   >  2:2585758/1‑91 (MQ=255)
                                                                          tCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTtata                  >  1:1474024/1‑91 (MQ=255)
                                                                           ccACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTatat                 >  2:495021/1‑91 (MQ=255)
                                                                           ccACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTatat                 >  1:2268257/1‑91 (MQ=255)
                                                                            cacacacaAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATAtt                <  1:3114849/91‑1 (MQ=255)
                                                                            cacaAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTCTCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATAtt                >  2:1107889/1‑91 (MQ=255)
                                                                             acaaacaaACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTg               <  2:1791429/91‑1 (MQ=255)
                                                                              gaaacaaccTATAATAATAGTTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGa              >  1:293950/2‑91 (MQ=38)
                                                                                  caaacTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaat          >  1:574909/1‑91 (MQ=255)
                                                                                    aacTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataa        <  1:3748480/91‑1 (MQ=255)
                                                                                     acTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaa       >  1:553164/1‑91 (MQ=255)
                                                                                     acTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAACATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaa       >  2:1495447/1‑91 (MQ=255)
                                                                                      cTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaat      >  2:788927/1‑91 (MQ=255)
                                                                                        attattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAataaataa    >  2:853220/1‑91 (MQ=255)
                                                                                          tattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAAc  >  1:1503965/1‑91 (MQ=255)
                                                                                          |                                                                                          
GATTGTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAAC  >  USA300TCH1516_ALE/2478040‑2478220

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