Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,221,660 T→C 8.7% intergenic (+75/‑141) rpoZ → / → coaBC DNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,221,6600TC 8.7% 617.5 / 36.4 213intergenic (+75/‑141)rpoZ/coaBCDNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC
Reads supporting (aligned to +/- strand):  ref base T (99/95);  new base C (7/12);  total (106/107)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 217 total aligned reads displayed.

TTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCAT  >  USA300TCH1516_ALE/1221570‑1221750
                                                                                          |                                                                                          
ttatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTAt                                                                                            <  1:1024094/91‑1 (MQ=255)
 tatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                           >  2:3326930/1‑91 (MQ=255)
 tatGGTAAAGAATAGTTTCGACCATTAAAAACATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                           >  1:3509017/1‑91 (MQ=255)
  atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                          >  2:3962909/1‑91 (MQ=255)
  atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                          >  2:2291441/1‑91 (MQ=255)
       aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                                     <  2:1196881/91‑1 (MQ=255)
       aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                                     <  1:1126152/91‑1 (MQ=255)
       aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtg                                                                                     <  1:1117713/91‑1 (MQ=255)
        aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                                    <  2:493716/91‑1 (MQ=255)
        aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                                    >  2:2986018/1‑91 (MQ=255)
        aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                    <  1:2095934/91‑2 (MQ=255)
        aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                    <  1:5735/91‑2 (MQ=255)
        aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                    <  1:991932/91‑2 (MQ=255)
         aGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                                   <  2:1029043/91‑1 (MQ=255)
         aGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                                   <  1:201329/91‑1 (MQ=255)
          gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                                  <  1:3598009/91‑1 (MQ=255)
            aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                                >  1:664053/1‑91 (MQ=255)
            aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                                >  2:3573719/1‑91 (MQ=255)
              aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCGAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttt                                                                              >  2:3245825/1‑91 (MQ=255)
               gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                             <  1:867888/91‑1 (MQ=255)
               gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                             >  2:3042927/1‑91 (MQ=255)
                 ttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTa                                                                           >  2:2814700/1‑91 (MQ=255)
                  tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                          <  2:2937727/91‑1 (MQ=255)
                  tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                          <  1:1861033/91‑1 (MQ=255)
                   cGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGt                                                                         >  2:1802695/1‑91 (MQ=255)
                   cGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGt                                                                         <  2:3136508/91‑1 (MQ=255)
                    gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                        >  2:1430759/1‑91 (MQ=255)
                    gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                        <  2:1083006/91‑1 (MQ=255)
                    gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                        >  2:2864973/1‑91 (MQ=255)
                       cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTAt                                                                     <  2:2141694/91‑1 (MQ=255)
                          tAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCa                                                                  <  1:1158271/91‑1 (MQ=255)
                          tAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCa                                                                  <  2:1493270/91‑1 (MQ=255)
                            aaaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                                >  1:595399/1‑91 (MQ=255)
                             aaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtg                                                               <  1:1213815/91‑1 (MQ=255)
                              aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                              <  2:2209528/91‑1 (MQ=255)
                              aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                              >  2:877278/1‑91 (MQ=255)
                              aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                              <  1:3275354/91‑1 (MQ=255)
                                tatGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGCTTTTTTAGTTTATCCAAtgatt                                                            >  2:3184380/1‑91 (MQ=255)
                                 atGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattg                                                           <  1:1072992/91‑1 (MQ=255)
                                  tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                          >  1:1755112/1‑91 (MQ=255)
                                  tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATACCCAAATCAAACGCCCCTATGTAGGTGTATTTTTTTAGTTTATCCAAtcattga                                                          >  1:3679492/1‑91 (MQ=255)
                                   gtaaaGCTAGACGTACATCAAATACGTCAAGGCATATCCATATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattgat                                                         <  2:1827057/88‑1 (MQ=255)
                                   gagAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattgat                                                         <  2:2635162/89‑1 (MQ=255)
                                      aaGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGtt                                                      <  1:1260865/91‑1 (MQ=255)
                                        gCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTtat                                                    <  1:2249165/91‑1 (MQ=255)
                                         cTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTtata                                                   <  2:3307094/91‑1 (MQ=255)
                                         cTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTtata                                                   <  2:2056174/91‑1 (MQ=255)
                                           aGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataat                                                 >  1:159307/1‑91 (MQ=255)
                                            gACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataata                                                <  1:881628/91‑1 (MQ=255)
                                              cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                              <  1:3837131/91‑1 (MQ=255)
                                              cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                              <  1:1754194/91‑1 (MQ=255)
                                              cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                              >  1:1885469/1‑91 (MQ=255)
                                              cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATGGATGTTataataat                                              >  1:3909882/1‑91 (MQ=255)
                                                tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                            <  2:504032/91‑1 (MQ=255)
                                                tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                            <  1:1787818/91‑1 (MQ=255)
                                                 aCATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACt                                           <  2:1010398/91‑1 (MQ=255)
                                                  cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                          >  2:722341/1‑91 (MQ=255)
                                                  cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                          <  1:3392889/91‑1 (MQ=255)
                                                  cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                          >  1:3814206/1‑91 (MQ=255)
                                                    tCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTaaa                                        >  2:732492/1‑91 (MQ=255)
                                                      aaaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                                      <  2:1532526/91‑1 (MQ=255)
                                                       aaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAttt                                     <  2:3822389/91‑1 (MQ=255)
                                                       aaTACGTCAAGGCATATCCAAATCAAACGCACATATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAttt                                     <  1:3355757/91‑1 (MQ=255)
                                                        aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                                    <  2:1714134/91‑1 (MQ=255)
                                                        aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                                    >  2:2202940/1‑91 (MQ=255)
                                                        aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                                    >  2:3767025/1‑91 (MQ=255)
                                                          aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                  >  2:756561/1‑91 (MQ=255)
                                                          aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                  >  1:2181748/1‑91 (MQ=255)
                                                          aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                  >  1:1144063/1‑91 (MQ=255)
                                                           cGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTAt                                 <  1:3090554/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                <  1:1510728/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                <  1:1511814/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                <  1:3022370/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                <  2:2587766/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                <  2:2554701/91‑1 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                >  2:2256365/1‑91 (MQ=255)
                                                            gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTAATAATAATACTAAATTTGTATc                                <  2:2996506/91‑1 (MQ=255)
                                                                 ggCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaa                           >  2:2169686/1‑91 (MQ=255)
                                                                  gCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaa                          <  1:3240845/91‑1 (MQ=255)
                                                                  gCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaa                          <  2:1363345/91‑1 (MQ=255)
                                                                   cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                         >  2:466571/1‑91 (MQ=255)
                                                                   cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                         >  1:308048/1‑91 (MQ=255)
                                                                    atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                        <  1:2823199/91‑1 (MQ=255)
                                                                     tatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaataaa                       >  1:3188624/1‑91 (MQ=255)
                                                                     tatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaaa                       >  2:1973611/1‑91 (MQ=255)
                                                                           aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg                 >  1:673709/1‑91 (MQ=255)
                                                                             tCAAACGCGCATATGTAGGTGTGTTGTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg               <  2:1656622/91‑1 (MQ=255)
                                                                             tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg               <  2:156876/91‑1 (MQ=255)
                                                                             tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg               >  2:780416/1‑91 (MQ=255)
                                                                              cAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGc              <  1:2541283/91‑1 (MQ=255)
                                                                               aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa             <  2:1596121/91‑1 (MQ=255)
                                                                                 acacACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt           >  1:2971576/4‑91 (MQ=255)
                                                                                  cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCCttt          >  1:3560403/1‑91 (MQ=255)
                                                                                  cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt          >  1:264958/1‑91 (MQ=255)
                                                                                   gCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtg         >  1:1903328/1‑91 (MQ=255)
                                                                                   acACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATGGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtg         >  1:2119897/2‑91 (MQ=255)
                                                                                     aCCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtgtg       >  2:238679/1‑91 (MQ=255)
                                                                                       cTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTgcg     >  2:304652/1‑91 (MQ=255)
                                                                                         aTGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCa   >  1:130415/1‑91 (MQ=255)
                                                                                         aTGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAATTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCa   >  2:3907161/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCAT  >  USA300TCH1516_ALE/1221570‑1221750

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