Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,478,133 T→A 6.7% intergenic (+20/‑555) 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,1330TA 6.7% 705.8 / 28.8 209intergenic (+20/‑555)USA300HOU_RS12760/USA300HOU_RS12765hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (110/85);  new base A (3/11);  total (113/96)
Fisher's exact test for biased strand distribution p-value = 1.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 232 total aligned reads displayed.

TGTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAACTAA  >  USA300TCH1516_ALE/2478043‑2478223
                                                                                          |                                                                                          
tGTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattatt                                                                                            <  1:2714660/91‑1 (MQ=255)
 gTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattatta                                                                                           >  2:1429789/1‑91 (MQ=255)
  ttattaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACtattattat                                                                                          <  1:2538832/91‑1 (MQ=255)
    attaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAg                                                                                        <  2:3046561/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)
     ttaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAgt                                                                                       <  1:261973/91‑2 (MQ=255)
     ttaACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAgg                                                                                       >  1:2418553/1‑91 (MQ=255)
       aaCGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgtt                                                                                     >  2:3212276/1‑91 (MQ=255)
        aCGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAACCTATAATAATAGtttg                                                                                    <  2:1974139/91‑1 (MQ=255)
         cGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgttgt                                                                                   >  1:1278930/1‑91 (MQ=255)
         cGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgttgt                                                                                   >  2:1320678/1‑91 (MQ=255)
          gTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTTCGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGgttgtt                                                                                  >  2:1048432/1‑91 (MQ=255)
            tataTTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTc                                                                                <  2:456939/91‑1 (MQ=255)
            tataTTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTc                                                                                >  1:3697310/1‑91 (MQ=255)
             atatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCt                                                                               <  2:3322353/91‑1 (MQ=255)
             atatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCt                                                                               <  2:2930017/91‑1 (MQ=255)
              tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                              <  1:1796350/91‑1 (MQ=255)
              tatTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTa                                                                              >  1:1960169/1‑91 (MQ=255)
               atTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTaa                                                                             <  2:41445/91‑1 (MQ=255)
               atTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTaa                                                                             <  1:3175304/91‑1 (MQ=255)
                 tGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAata                                                                           >  1:2808204/1‑91 (MQ=255)
                  gTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAatat                                                                          <  2:122918/91‑1 (MQ=255)
                   ttCAAGTAATTGGAATAAATGTCTTTACGCCGCTATTATTTATAATAGTTAATACACAAACAAACTATTAATATAGGAAGTTTataataat                                                                         >  2:3627930/1‑89 (MQ=255)
                       agtagtTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGt                                                                     >  2:709169/1‑91 (MQ=255)
                         tagtTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCa                                                                   >  1:1047624/1‑91 (MQ=255)
                           gtTGGTTTATATGTCTTTACGCCGCTTTTGTTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCaaa                                                                 <  2:741657/91‑1 (MQ=255)
                            ttGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATtata                                                                                         >  2:2670439/1‑66 (MQ=255)
                            ttGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATtata                                                                                         <  1:2670439/66‑1 (MQ=255)
                            ttGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAAc                                                                >  1:27516/1‑91 (MQ=255)
                             tGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACt                                                               >  1:820337/1‑91 (MQ=255)
                             tGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACt                                                               <  2:2255302/91‑1 (MQ=255)
                              ggTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTa                                                              >  1:2386986/1‑91 (MQ=255)
                               gTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTAt                                                             >  2:1849070/1‑91 (MQ=255)
                                 ttATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCt                                                           >  2:2062437/1‑91 (MQ=255)
                                   atatGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCttt                                                         <  1:1674612/91‑1 (MQ=255)
                                   atatGTCGTTACGCCGCTTTTATTTATAATAGTGAATCCACAAACAAACGAGGATTATAGGTGGTTTCTAAGATTGCTTCAAACTATCttt                                                         <  1:3084912/91‑1 (MQ=255)
                                    tatGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTa                                                        >  2:1459483/1‑91 (MQ=255)
                                       gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                     <  1:1598046/91‑1 (MQ=255)
                                       gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                     >  2:499275/1‑91 (MQ=255)
                                       gTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAat                                                     >  1:3665624/1‑91 (MQ=255)
                                        tCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAata                                                    >  1:1690986/1‑91 (MQ=255)
                                         cTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAatat                                                   <  1:3529997/91‑1 (MQ=255)
                                         cTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAatat                                                   <  2:3077726/91‑1 (MQ=255)
                                          tttACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATAtt                                                  >  1:2305163/1‑91 (MQ=255)
                                           ttACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTc                                                 >  1:958648/1‑91 (MQ=255)
                                           ttACGCCGCGTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTc                                                 <  2:476888/91‑1 (MQ=255)
                                            tACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCg                                                >  2:201972/1‑91 (MQ=255)
                                            tACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCg                                                <  1:2056433/91‑1 (MQ=255)
                                             aCGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGt                                               <  1:3023226/91‑1 (MQ=255)
                                              cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                              >  1:781084/1‑91 (MQ=255)
                                              cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                              <  2:859353/91‑1 (MQ=255)
                                              cgccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGtt                                              >  2:3360242/1‑91 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            <  2:2910689/91‑1 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            >  2:2755335/1‑91 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            <  2:112221/91‑1 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            <  1:2965808/91‑1 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            >  2:567341/1‑91 (MQ=255)
                                                ccgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTaa                                            >  1:1442233/1‑91 (MQ=255)
                                                 cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                           <  2:3462709/91‑1 (MQ=255)
                                                 cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                           <  2:152650/91‑1 (MQ=255)
                                                 cgcTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAAt                                           <  1:113625/91‑1 (MQ=255)
                                                    ttttATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCt                                        >  2:259437/1‑91 (MQ=255)
                                                      ttatttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTg                                      >  1:3152010/1‑91 (MQ=255)
                                                        atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGAt                                    <  1:1031578/91‑1 (MQ=255)
                                                        atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATAGTCGTTAATCCTTGAt                                    <  1:1033028/91‑1 (MQ=255)
                                                        atttatAATAGTTAATCCACAAACAAACTATTATTATAGGTAGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGAt                                    >  2:2943964/1‑91 (MQ=255)
                                                          ttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTa                                  <  1:3417589/91‑1 (MQ=255)
                                                          ttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTa                                  >  1:237024/1‑91 (MQ=255)
                                                          ttatAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTCAATCCTTGATTa                                  >  2:3234236/1‑91 (MQ=255)
                                                          ttatAATAGTTAATCCACAAACAAACTATTATTATAGCTTGTTTCTAATATCGCGTCAAACTATCGTTAAATATTCGTTAATCCTCGATTa                                  >  1:3612766/1‑91 (MQ=255)
                                                            ataataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGa                                <  2:887305/91‑1 (MQ=255)
                                                            ataataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGa                                <  1:2902389/91‑1 (MQ=255)
                                                              aataGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACa                              <  1:834549/91‑1 (MQ=255)
                                                                 aGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAAtt                           >  2:2872767/1‑91 (MQ=255)
                                                                 aGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAAtt                           <  1:2007755/91‑1 (MQ=255)
                                                                  gTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTg                          <  1:700773/91‑1 (MQ=255)
                                                                   ttAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGt                         >  1:1711483/1‑91 (MQ=255)
                                                                    tAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGtt                        >  1:3832831/1‑91 (MQ=255)
                                                                     aaTCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTta                       <  1:101664/91‑1 (MQ=255)
                                                                      aTCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTtat                      >  2:2586089/1‑91 (MQ=255)
                                                                       tCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTtata                     >  1:1474024/1‑91 (MQ=255)
                                                                        ccACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTatat                    >  1:2268257/1‑91 (MQ=255)
                                                                        ccACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTatat                    >  2:495021/1‑91 (MQ=255)
                                                                         cacaAACAAACTATTATTATAGGTTGTTTATAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATAtt                   <  1:3115513/91‑1 (MQ=255)
                                                                             aacaaccTATAATAATAGTTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGaaa               >  2:1558939/1‑91 (MQ=38)
                                                                             aacaaacTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGaaa               >  2:1684926/1‑91 (MQ=255)
                                                                               caaacTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaat             >  1:574909/1‑91 (MQ=255)
                                                                                 aacTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataa           <  1:3748480/91‑1 (MQ=255)
                                                                                  acTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaa          >  1:553164/1‑91 (MQ=255)
                                                                                  acTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATAGTTATATTGAaaataaa          >  2:1024457/1‑91 (MQ=255)
                                                                                  acTATTATTATAGGTTGTTTCTAATATTGCGTCAAACAATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaa          >  1:3268607/1‑91 (MQ=255)
                                                                                   cTATTATTATATGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaat         >  2:1748079/1‑91 (MQ=255)
                                                                                   cTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAaaataaat         >  2:788927/1‑91 (MQ=255)
                                                                                   atattattatAGGTTGTTTCTAATATTGCGTCAAACTATCTGTAAATATTCGTTAATCCTTGATTAGACAAGTGTTATATTGAaaataaat         <  2:3409804/90‑1 (MQ=255)
                                                                                    tattattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAaataaata        <  2:3160649/91‑1 (MQ=255)
                                                                                     attattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAataaataa       >  2:1778273/1‑91 (MQ=255)
                                                                                     attattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAataaataa       >  2:853220/1‑91 (MQ=255)
                                                                                       tattatAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAAc     >  1:1503965/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGTTATTAACGTTATATTGTTCAAGTAGTTGGTTTATATGTCTTTACGCCGCTTTTATTTATAATAGTTAATCCACAAACAAACTATTATTATAGGTTGTTTCTAATATTGCGTCAAACTATCTTTAAATATTCGTTAATCCTTGATTAGACAATTGTTATATTGAAAATAAATAAACTAA  >  USA300TCH1516_ALE/2478043‑2478223

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: