Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,909,137 T→A 5.2% intergenic (‑203/‑57) tal ← / → USA300HOU_RS09460 Transaldolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,909,1370TA 5.2% 911.5 / 27.2 267intergenic (‑203/‑57)tal/USA300HOU_RS09460Transaldolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (106/146);  new base A (2/12);  total (108/159)
Fisher's exact test for biased strand distribution p-value = 4.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 269 total aligned reads displayed.

TCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCT  >  USA300TCH1516_ALE/1909047‑1909226
                                                                                          |                                                                                         
 cACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTg                                                                                          <  2:2719643/91‑1 (MQ=255)
  aCGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGt                                                                                         <  2:2511032/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATGAAACAATAGttt                                                                                       <  1:3047716/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                       <  2:2186370/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                       <  2:1327972/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                       <  2:1328038/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                       <  1:3727929/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGttt                                                                                       <  1:3375505/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                       <  2:1268710/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                       <  1:1516672/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                       <  1:3444134/91‑1 (MQ=255)
    gCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGttt                                                                                       <  1:1732740/91‑1 (MQ=255)
     gttttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                      <  2:1533142/90‑1 (MQ=255)
     cTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGtttt                                                                                      <  1:3598722/91‑1 (MQ=255)
     cTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGgttt                                                                                      <  2:3504194/91‑1 (MQ=255)
      ttttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACAATAGTTTTa                                                                                     <  2:1026022/91‑1 (MQ=255)
        ttATAATATATAGCTAATTTTTTGCCTTTATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                   >  2:3103006/1‑91 (MQ=255)
        ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                   <  1:250883/91‑1 (MQ=255)
        ttATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATg                                                                                   <  1:1381103/91‑1 (MQ=255)
         tataATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGa                                                                                  >  2:3591584/1‑91 (MQ=255)
          ataataTATAGCTAATTTTTTGCCTATTTTTTTTAAAAATTACTTTAGCACACATTGTGTTTACAAGTTAATATTACTATTGTTTTATGAt                                                                                 >  2:58861/1‑91 (MQ=255)
          ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                                 >  2:2111444/1‑91 (MQ=255)
          ataataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAt                                                                                 >  1:2964996/1‑91 (MQ=255)
           taataTATAGCTAATTTTTTTTTTTTTTTTTTTAAAAATGACTTTAGCAAACATTTTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                                >  2:2716523/1‑91 (MQ=255)
           taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                                <  1:2207102/91‑1 (MQ=255)
           taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                                >  2:3215067/1‑91 (MQ=255)
           taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                                <  2:2317517/91‑1 (MQ=255)
           taataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGAtt                                                                                >  1:3473552/1‑91 (MQ=255)
            aataTATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTg                                                                               >  1:3427415/1‑91 (MQ=255)
               atataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGttt                                                                            <  1:460684/91‑1 (MQ=255)
                tataGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTa                                                                           <  2:2825931/91‑1 (MQ=255)
                   aGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGa                                                                        <  2:289885/91‑1 (MQ=255)
                    gCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGAt                                                                       >  2:645484/1‑91 (MQ=255)
                    gCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGAt                                                                       <  2:308240/91‑1 (MQ=255)
                     cTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATa                                                                      >  1:2063354/1‑91 (MQ=255)
                      tAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGTTTGTTTATGATAc                                                                     >  1:1662658/1‑91 (MQ=255)
                       aaTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACg                                                                    <  1:2144184/91‑1 (MQ=255)
                       aaTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACg                                                                    >  1:1432535/1‑91 (MQ=255)
                        aTTTTTTGCCTATATTTTTTAAATATTACTTTATCCAACATTGTGTCTACAAGTTATTAAAACTATTTTTTTATGATTGTTTATGATAcca                                                                   >  2:3373550/1‑89 (MQ=255)
                        aTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGa                                                                   <  2:3268553/91‑1 (MQ=255)
                        aTTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGa                                                                   <  1:3888495/91‑1 (MQ=255)
                          tttttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGata                                                                 <  2:1692029/91‑1 (MQ=255)
                           ttttGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGatat                                                                >  2:48567/1‑91 (MQ=255)
                           ttttGCCTATATTTTTGAAGAATGACTTTGGCAAACATGGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGatat                                                                <  2:2573234/91‑1 (MQ=255)
                              tGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTa                                                             <  1:3348266/91‑1 (MQ=255)
                               gCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaa                                                            >  1:573999/1‑91 (MQ=255)
                                ccTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaaa                                                           <  1:931476/91‑1 (MQ=255)
                                ccTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaaa                                                           >  1:821778/1‑91 (MQ=255)
                                ccTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTaaa                                                           <  1:741565/91‑1 (MQ=255)
                                  tataTTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATg                                                         <  2:3147924/91‑1 (MQ=255)
                                    tatTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGAt                                                       <  2:2399242/91‑1 (MQ=255)
                                      ttttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATaa                                                     >  1:129259/1‑91 (MQ=255)
                                        ttttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatc                                                   >  2:1939095/1‑91 (MQ=255)
                                        ttttAAAAATGACTTTAGCAAACATTGGGTCTACCAGGTAATAAAACTATTGTTTTATGATTGTTTATGATACGATTTGTAAATGATAatc                                                   <  2:1418684/91‑1 (MQ=255)
                                         tttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatca                                                  >  1:1606979/1‑91 (MQ=255)
                                          ttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcat                                                 >  1:1087188/1‑91 (MQ=255)
                                          ttAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAAATATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcat                                                 <  1:3134866/91‑1 (MQ=255)
                                           tAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAatcatc                                                <  2:1198753/91‑1 (MQ=255)
                                            aaaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATcc                                               <  2:1199441/91‑1 (MQ=255)
                                             aaaaTGAGTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCa                                              <  1:1968061/91‑1 (MQ=255)
                                             aaaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCa                                              >  2:3608915/1‑91 (MQ=255)
                                              aaaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCag                                             <  2:1810019/91‑1 (MQ=255)
                                               aaTGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCagg                                            >  1:90048/1‑91 (MQ=255)
                                                 tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                          <  1:1002621/91‑1 (MQ=255)
                                                 tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                          <  1:496772/91‑1 (MQ=255)
                                                 tGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCaggag                                          <  1:1846598/91‑1 (MQ=255)
                                                   aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                        <  2:1634955/91‑1 (MQ=255)
                                                   aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                        <  1:3644104/91‑1 (MQ=255)
                                                   aCTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGt                                        <  1:3643922/91‑1 (MQ=255)
                                                      ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTa                                     <  2:1572842/91‑1 (MQ=255)
                                                      ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTa                                     >  2:2688078/1‑91 (MQ=255)
                                                      ttAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTa                                     >  1:2890186/1‑91 (MQ=255)
                                                       tAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAg                                    >  2:2628619/1‑91 (MQ=255)
                                                       tAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAg                                    >  2:2040596/1‑91 (MQ=255)
                                                         gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                                  <  1:2018732/91‑1 (MQ=255)
                                                         gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                                  >  1:3787819/1‑91 (MQ=255)
                                                         gCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCa                                  >  2:2986997/1‑91 (MQ=255)
                                                             aCATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtc                              >  2:2426482/1‑91 (MQ=255)
                                                              cATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtcg                             >  2:3922631/1‑90 (MQ=255)
                                                               aTTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGtctc                            <  1:1944712/91‑1 (MQ=255)
                                                                ttGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCg                           <  1:1865696/91‑1 (MQ=255)
                                                                 tgtgTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGt                          <  1:2369933/91‑1 (MQ=255)
                                                                  gtgtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGtt                         >  2:377564/1‑91 (MQ=255)
                                                                    gtCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCa                       >  2:94591/1‑91 (MQ=255)
                                                                      cTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaa                     <  2:1969278/91‑1 (MQ=255)
                                                                        aCAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaa                   >  2:2358327/1‑91 (MQ=255)
                                                                        aCAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaa                   >  2:3403710/1‑91 (MQ=255)
                                                                        aCAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaa                   <  1:2306673/91‑1 (MQ=255)
                                                                         ctaGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaaa                  <  1:983023/89‑1 (MQ=255)
                                                                         cAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCaaaaaa                  <  1:2603408/91‑1 (MQ=255)
                                                                            gTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATa               >  1:329155/1‑91 (MQ=255)
                                                                              tcataaAACAATAGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACt             >  2:104659/3‑91 (MQ=255)
                                                                                 taaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt          <  2:1350177/91‑1 (MQ=255)
                                                                                 taaAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACtttt          >  2:695355/1‑91 (MQ=255)
                                                                                   aaaCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGCGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTAc        >  1:1468116/1‑91 (MQ=255)
                                                                                   aaaCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTAc        <  1:1714469/91‑1 (MQ=255)
                                                                                   aaaCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTATTAc        >  2:1477199/1‑91 (MQ=255)
                                                                                     aCTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACtt      >  2:1057258/1‑91 (MQ=255)
                                                                                       tATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTAt    <  1:2575598/91‑1 (MQ=255)
                                                                                         ttGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCt  >  2:864981/1‑91 (MQ=255)
                                                                                          |                                                                                         
TCACGCTTTTATAATATATAGCTAATTTTTTGCCTATATTTTTTAAAAATGACTTTAGCAAACATTGTGTCTACAAGTTAATAAAACTATTGTTTTATGATTGTTTATGATACGATATGTAAATGATAATCATCCAGGAGGTCTAGCATGTCTCGTTCAAAAAAATACTTTTACTTATCT  >  USA300TCH1516_ALE/1909047‑1909226

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