Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,309,807 Δ1 bp 10.4% coding (62/1032 nt) yfiZ_2 ← putative siderophore transport system permease protein YfiZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,309,8070G.10.4% 771.6 / 84.4 231coding (62/1032 nt)yfiZ_2putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand):  ref base G (107/100);  new base . (9/15);  total (116/115)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 242 total aligned reads displayed.

TATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTT  >  USA300TCH1516_ALE/2309718‑2309897
                                                                                           |                                                                                          
   cGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGtgttgt                                                                                        <  2:2789696/91‑1 (MQ=255)
    gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTtgtg                                                                                       <  1:520259/91‑1 (MQ=25)
    gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTtgtg                                                                                       >  2:772813/1‑91 (MQ=255)
     cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGc                                                                                      <  2:2664332/91‑1 (MQ=255)
      tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCg                                                                                     >  2:3837805/1‑91 (MQ=255)
      tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:252818/91‑3 (MQ=17)
      tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:2468624/91‑3 (MQ=17)
      tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  2:2117023/91‑3 (MQ=17)
       aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  2:1603282/91‑3 (MQ=21)
       aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTTGTGctt                                                                                    <  1:2333842/91‑3 (MQ=21)
       aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGt                                                                                    <  1:3736597/91‑1 (MQ=255)
       aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGcttt                                                                                   <  1:3636025/91‑4 (MQ=12)
        aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtc                                                                                   <  1:275243/91‑1 (MQ=255)
        aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctttc                                                                                  <  2:557892/91‑5 (MQ=9)
        aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGACGCACAACACTCACGTATATAGTGAGTTTGTGctttc                                                                                  <  1:3487516/91‑5 (MQ=9)
         aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATAGTGAGTTTGTGctttc                                                                                  <  2:2423895/91‑5 (MQ=21)
         aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtct                                                                                  <  1:3720102/91‑1 (MQ=255)
              attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                             <  1:3852239/91‑1 (MQ=255)
              attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                             <  2:379097/91‑1 (MQ=255)
              attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttc                                                                             >  1:768196/1‑91 (MQ=255)
               ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATCCGTGAGTGTTGTGCGtctttttct                                                                            >  2:3682411/1‑91 (MQ=255)
               ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                            >  2:20820/1‑91 (MQ=255)
               ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttct                                                                            >  1:1132302/1‑91 (MQ=255)
                tttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctt                                                                           >  1:3550661/1‑91 (MQ=255)
                  tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                         <  1:874771/91‑1 (MQ=255)
                  tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGtctttttctttt                                                                         >  1:242024/1‑91 (MQ=255)
                    gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTAt                                                                       <  2:2831337/91‑1 (MQ=255)
                     aagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATg                                                                      <  1:56633/91‑1 (MQ=255)
                      agaaCCAATCGCCATATTTAAATATATACGAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGa                                                                     <  2:511500/91‑1 (MQ=255)
                       gaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGAt                                                                    >  2:2951750/1‑91 (MQ=255)
                       gaaCCAATCGCCATATTTAAATATATACAAAAAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGAt                                                                    >  1:3512689/1‑91 (MQ=255)
                         aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCt                                                                  <  1:1612811/91‑1 (MQ=255)
                          ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                                 >  1:2639021/1‑91 (MQ=255)
                          ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                                 >  1:308611/1‑91 (MQ=255)
                             aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATgg                                                              <  2:2229970/91‑1 (MQ=255)
                              tCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGc                                                             >  2:3842766/1‑91 (MQ=255)
                              tCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGc                                                             <  1:433185/91‑1 (MQ=255)
                                   atatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGAtt                                                        >  1:197711/1‑91 (MQ=255)
                                   atatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGAtt                                                        <  1:688681/91‑1 (MQ=255)
                                   atatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGCGTTGTGCGTCTTTTTCTATTATGATCTATGGCAGAtt                                                        >  1:1043593/1‑91 (MQ=255)
                                    tatTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTg                                                       <  2:2720261/91‑1 (MQ=255)
                                      tttAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGc                                                     <  1:158162/91‑1 (MQ=255)
                                       ttAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                                    <  2:185554/91‑1 (MQ=255)
                                       ttAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                                    >  2:1486773/1‑91 (MQ=255)
                                        tAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTa                                                   <  2:1549342/91‑1 (MQ=255)
                                         aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  >  2:470797/1‑91 (MQ=255)
                                         aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  <  1:1269482/91‑1 (MQ=255)
                                         aaaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                                  >  2:1307787/1‑91 (MQ=255)
                                          aaTATATACAAATAAGTAGAAAGCACAAAGTCACTATATTCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                                 <  2:3866868/91‑1 (MQ=255)
                                          aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                                 <  2:1126498/91‑1 (MQ=255)
                                          aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                                 >  1:205486/1‑91 (MQ=255)
                                          aaTATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                                 >  1:3600779/1‑91 (MQ=255)
                                           atatatACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACtt                                                <  1:431294/91‑1 (MQ=255)
                                            tatataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                               >  1:2092059/1‑91 (MQ=255)
                                            tatataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                               >  1:2939283/1‑91 (MQ=255)
                                             atataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAt                                              <  1:277234/91‑1 (MQ=255)
                                               ataCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                                            >  2:2780124/1‑91 (MQ=255)
                                                 aCAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCt                                          <  1:1389189/91‑1 (MQ=255)
                                                   aaaTAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCttt                                        >  1:588176/1‑91 (MQ=255)
                                                       aaaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                                    <  1:2672384/91‑1 (MQ=255)
                                                        aaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAg                                   >  1:2427784/1‑91 (MQ=255)
                                                        aaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAGTTCTTTTATAg                                   >  1:3584696/1‑91 (MQ=255)
                                                         aaagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGc                                  <  1:2210745/91‑1 (MQ=255)
                                                         aaagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATCTCTTTCATAGc                                  >  2:1832044/1‑91 (MQ=255)
                                                          aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                                 <  1:3335573/91‑1 (MQ=255)
                                                          aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                                 <  1:1521376/91‑1 (MQ=255)
                                                          aagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGca                                 <  2:1401614/91‑2 (MQ=255)
                                                           agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                                <  1:3307957/91‑1 (MQ=255)
                                                            gaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAt                               >  2:462479/1‑91 (MQ=255)
                                                             aaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                              <  1:1064512/91‑1 (MQ=255)
                                                             aaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                              <  1:3228388/91‑1 (MQ=255)
                                                              aagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt                             >  2:253248/1‑91 (MQ=255)
                                                               agCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTa                            <  2:63598/91‑1 (MQ=255)
                                                                gCACAACACTCACGTATATAGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                           >  1:2539788/1‑91 (MQ=25)
                                                                gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                           >  2:1286443/1‑91 (MQ=255)
                                                                gCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc                           >  2:3214787/1‑91 (MQ=255)
                                                                  acaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt                         >  2:3136810/1‑91 (MQ=255)
                                                                  acaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt                         >  1:341470/1‑91 (MQ=255)
                                                                  acaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt                         >  2:3067490/1‑91 (MQ=255)
                                                                   caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc                        >  1:2238552/1‑91 (MQ=255)
                                                                       aCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc                     >  1:1875147/1‑91 (MQ=255)
                                                                         tCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTCATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt                   <  2:3907542/91‑1 (MQ=255)
                                                                          cacgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa                 >  1:4036/5‑91 (MQ=255)
                                                                          cacgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa                 >  2:3250788/5‑91 (MQ=255)
                                                                           acgtatataGTGAGTTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc                >  1:2563624/4‑91 (MQ=255)
                                                                           aCTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGTAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa                 <  2:3492072/91‑1 (MQ=255)
                                                                           aCTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa                 >  2:1370446/1‑91 (MQ=255)
                                                                            cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATTGCAGATTGTCTACTTATTTCTTTTATAGCCATTTACTTTCACCtttc                >  2:820275/1‑89 (MQ=255)
                                                                            cTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc                >  2:736639/1‑91 (MQ=255)
                                                                                tataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACttt             <  2:2680949/91‑1 (MQ=255)
                                                                                  taCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg           >  2:305658/1‑91 (MQ=255)
                                                                                    cGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGtt         >  2:3584155/1‑91 (MQ=255)
                                                                                     gTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTc        <  2:1620726/91‑1 (MQ=255)
                                                                                     gTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTc        >  2:2021997/1‑91 (MQ=255)
                                                                                      tGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCt       <  1:1385298/91‑1 (MQ=255)
                                                                                      tGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCt       <  2:592269/91‑1 (MQ=255)
                                                                                       gAGTGTTGTGCGTCTTTTTTTTTTATTTTCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCtt      >  2:3088317/1‑91 (MQ=255)
                                                                                       gAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCtt      <  1:3822300/91‑1 (MQ=255)
                                                                                         gtgtTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCtttt    >  1:1085042/1‑91 (MQ=255)
                                                                                           gttgtGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCtttttt  >  1:3906016/1‑91 (MQ=255)
                                                                                           |                                                                                          
TATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTT  >  USA300TCH1516_ALE/2309718‑2309897

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: