Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,221,653 G→A 8.8% intergenic (+68/‑148) 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,6530GA 8.8% 312.6 / 16.2 102intergenic (+68/‑148)rpoZ/coaBCDNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC
Reads supporting (aligned to +/- strand):  ref base G (42/51);  new base A (2/7);  total (44/58)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 103 total aligned reads displayed.

TGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221564‑1221743
                                                                                         |                                                                                          
tGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGc                                                                                           <  1:2770545/91‑1 (MQ=255)
tGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATATGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGc                                                                                           >  1:3594078/1‑91 (MQ=255)
tGAGTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGc                                                                                           <  1:3262929/91‑1 (MQ=255)
 gATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATTCAAATCAAACGCa                                                                                          <  1:2930139/91‑1 (MQ=255)
 gATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCa                                                                                          <  1:3697460/91‑1 (MQ=255)
   ttattaTGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAcc                                                                                        >  1:3977870/1‑91 (MQ=255)
   ttattaTGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAcc                                                                                        <  1:1253693/91‑1 (MQ=255)
   ttattaTGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAcc                                                                                        <  1:3919151/91‑1 (MQ=255)
      ttatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTAt                                                                                     >  1:2187968/1‑91 (MQ=255)
       tatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                    <  1:1626077/91‑1 (MQ=255)
        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                   >  1:203578/1‑91 (MQ=255)
        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  1:2039841/91‑4 (MQ=255)
        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  1:1623094/91‑4 (MQ=255)
        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  1:1623080/91‑4 (MQ=255)
        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  1:2367638/91‑4 (MQ=255)
         tGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTacata                                                                                  <  1:796549/91‑5 (MQ=255)
         tGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTacata                                                                                  <  1:796764/91‑5 (MQ=255)
         tGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacata                                                                                  <  1:217092/91‑5 (MQ=255)
          ggTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAg                                                                                 >  1:2379117/1‑91 (MQ=255)
             aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                              <  1:1000163/91‑1 (MQ=255)
             aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                              >  1:3118735/1‑91 (MQ=255)
             aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                              <  1:3052916/91‑1 (MQ=255)
              aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                             <  1:3975252/91‑1 (MQ=255)
              aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                             <  1:1534934/91‑1 (MQ=255)
              aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                             <  1:990285/91‑2 (MQ=255)
              aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                             <  1:990338/91‑2 (MQ=255)
               aggaTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                            <  1:1543973/88‑1 (MQ=255)
               aGAATAGTTACGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                            >  1:2089606/1‑91 (MQ=255)
                gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                           >  1:743965/1‑91 (MQ=255)
                gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                           >  1:1449870/1‑91 (MQ=255)
                   tAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttt                                                                        >  1:2863089/1‑91 (MQ=255)
                     gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  1:3745926/91‑1 (MQ=255)
                     gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  1:1808364/91‑1 (MQ=255)
                     gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  1:2509473/91‑1 (MQ=255)
                      tttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttttt                                                                     <  1:1398993/91‑1 (MQ=255)
                        tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   <  1:2867988/91‑1 (MQ=255)
                        tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   >  1:1941274/1‑91 (MQ=255)
                        tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   >  1:2417460/1‑91 (MQ=255)
                        tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   >  1:2417443/1‑91 (MQ=255)
                          gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                 >  1:1638035/1‑91 (MQ=255)
                          gACCATTAAAAATATATGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTTGGTGTGTTTTTTTAGtt                                                                 >  1:3545560/1‑91 (MQ=255)
                           aCCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGttt                                                                <  1:3027009/91‑1 (MQ=255)
                            ccATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTa                                                               >  1:1210667/1‑91 (MQ=255)
                             cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTAt                                                              <  1:1964758/91‑1 (MQ=255)
                             cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCCATGTAGGTGTGTTTTTTTAGTTTAt                                                              >  1:3865515/1‑91 (MQ=255)
                             cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCAATGTAGGTGTGTTTTTTTAGTTTAt                                                              >  1:3865450/1‑91 (MQ=255)
                               ttAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATcc                                                            <  1:487913/91‑1 (MQ=255)
                               ttAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGGGTTTTTTTAGTTTATcc                                                            >  1:907339/1‑91 (MQ=255)
                                  aaaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                         >  1:84848/1‑91 (MQ=255)
                                   aaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtg                                                        >  1:3079044/1‑91 (MQ=255)
                                     atatGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgat                                                      <  1:1015044/91‑1 (MQ=255)
                                       atGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattg                                                    <  1:828913/91‑1 (MQ=255)
                                        tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  1:1956360/1‑91 (MQ=255)
                                        tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  1:2563906/1‑91 (MQ=255)
                                            aaGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGtt                                               >  1:388020/1‑91 (MQ=255)
                                                 aGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataat                                          <  1:1819679/91‑1 (MQ=255)
                                                 aGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataat                                          <  1:1964911/91‑1 (MQ=255)
                                                  gACGTACATCAAATACGTCAAGGCATATCCAAATCAAATGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataata                                         >  1:1181898/1‑91 (MQ=255)
                                                   aCGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataa                                        <  1:2592337/91‑1 (MQ=255)
                                                    cGTACATCAAATACGTCAAGGCATATCCCAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       <  1:1693473/91‑1 (MQ=255)
                                                      tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                     <  1:2771612/91‑1 (MQ=255)
                                                        cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                   >  1:3237322/1‑91 (MQ=255)
                                                        cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                   >  1:2405881/1‑91 (MQ=255)
                                                           cAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAt                                >  1:2706175/1‑91 (MQ=255)
                                                            aaaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                               <  1:2889067/91‑1 (MQ=255)
                                                              aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                             <  1:3349611/91‑1 (MQ=255)
                                                              aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                             <  1:3350479/91‑1 (MQ=255)
                                                              aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                             <  1:2177534/91‑1 (MQ=255)
                                                                aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                           <  1:1636647/91‑1 (MQ=255)
                                                                 cGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGTTGTTATAATAATACTAAATTTGTAt                          >  1:3873144/1‑91 (MQ=255)
                                                                  gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:3156644/91‑1 (MQ=255)
                                                                  gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:1945834/91‑1 (MQ=255)
                                                                  gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:1219788/91‑1 (MQ=255)
                                                                  gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         >  1:429553/1‑91 (MQ=255)
                                                                  gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:3256679/91‑1 (MQ=255)
                                                                    cAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCta                       >  1:2327594/1‑91 (MQ=255)
                                                                    cAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCta                       >  1:548410/1‑91 (MQ=255)
                                                                     aaGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtat                      >  1:416243/1‑91 (MQ=255)
                                                                       ggCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaa                    <  1:3667338/91‑1 (MQ=255)
                                                                          atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:3631273/91‑1 (MQ=255)
                                                                          atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:361408/91‑1 (MQ=255)
                                                                          atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:260980/91‑1 (MQ=255)
                                                                          atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:261124/91‑1 (MQ=255)
                                                                           tatCCAAATCAAACGCACCTATGTAGGTGTGCTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaaa                >  1:489454/1‑91 (MQ=255)
                                                                            atCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAg               <  1:42254/91‑1 (MQ=255)
                                                                                aaaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAAt           >  1:3803768/1‑91 (MQ=255)
                                                                                aaaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAAt           >  1:3803842/1‑91 (MQ=255)
                                                                                 aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          >  1:2267318/1‑91 (MQ=255)
                                                                                  aaaaaaCACACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGa         >  1:1133369/4‑91 (MQ=37)
                                                                                  aTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGa         <  1:1686231/91‑1 (MQ=255)
                                                                                  aTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGa         >  1:60784/1‑91 (MQ=255)
                                                                                   tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        <  1:2479603/91‑1 (MQ=255)
                                                                                   tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        <  1:2494598/91‑1 (MQ=255)
                                                                                   aaaaaCACACCTACATAGGCGCGTTTTTTTAGTTTATACAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        >  1:3495513/3‑91 (MQ=255)
                                                                                     aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa      <  1:314065/91‑1 (MQ=255)
                                                                                     aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa      <  1:308177/91‑1 (MQ=255)
                                                                                       aCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt    <  1:158964/91‑1 (MQ=255)
                                                                                        cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt   >  1:1541431/1‑91 (MQ=255)
                                                                                        cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt   >  1:1541396/1‑91 (MQ=255)
                                                                                         gCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtg  >  1:639840/1‑91 (MQ=255)
                                                                                         |                                                                                          
TGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221564‑1221743

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: