Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,221,653 G→A 7.1% 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 7.1% 626.0 / 23.2 196intergenic (+68/‑148)rpoZ/coaBCDNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC
Reads supporting (aligned to +/- strand):  ref base G (85/97);  new base A (2/12);  total (87/109)
Fisher's exact test for biased strand distribution p-value = 2.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 205 total aligned reads displayed.

GTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221563‑1221743
                                                                                          |                                                                                          
gTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACg                                                                                            <  2:2773677/91‑1 (MQ=255)
gTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACg                                                                                            <  1:3231223/91‑1 (MQ=255)
  gATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCa                                                                                          >  2:3475208/1‑91 (MQ=255)
   attattATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAc                                                                                         >  1:1018288/1‑91 (MQ=255)
    ttattaTGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAcc                                                                                        <  1:3474385/91‑1 (MQ=255)
     tattatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCt                                                                                       >  2:3220240/1‑91 (MQ=255)
      attatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTa                                                                                      >  2:1770377/1‑91 (MQ=255)
      attatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTa                                                                                      >  2:1791766/1‑91 (MQ=255)
      attatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTa                                                                                      <  2:1103061/91‑1 (MQ=255)
      attatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTa                                                                                      >  1:1704747/1‑91 (MQ=255)
       ttatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTAt                                                                                     >  2:1321480/1‑91 (MQ=255)
       ttatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTAt                                                                                     >  2:2503387/1‑91 (MQ=255)
         atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                   >  1:1624980/1‑91 (MQ=255)
         atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  1:1751435/91‑4 (MQ=255)
         atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  2:3179638/91‑4 (MQ=255)
            gTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAgg                                                                                <  1:3592703/91‑1 (MQ=255)
            gTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAgg                                                                                >  2:2672861/1‑91 (MQ=255)
              aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtg                                                                              <  1:3660695/91‑1 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAGTCAAACGCCCCTATGTAGgtgt                                                                             <  1:2570038/91‑1 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                             >  2:3545007/1‑91 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                             <  1:2623768/91‑1 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                             <  1:980206/91‑2 (MQ=255)
                 gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                           <  1:2551894/91‑1 (MQ=255)
                 gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                           <  1:955359/91‑1 (MQ=255)
                 gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgcgt                                                                           <  1:653850/91‑4 (MQ=255)
                   aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                         <  2:1930003/91‑1 (MQ=255)
                     aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttt                                                                       <  2:1407126/91‑1 (MQ=255)
                     aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttt                                                                       <  1:311230/91‑1 (MQ=255)
                      gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  2:3580600/91‑1 (MQ=255)
                      gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  1:3299118/91‑1 (MQ=255)
                      gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      >  1:1347849/1‑91 (MQ=255)
                       tttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttttt                                                                     <  1:3850275/91‑1 (MQ=255)
                       tttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttttt                                                                     <  2:3376182/91‑1 (MQ=255)
                         tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   <  1:1475781/91‑1 (MQ=255)
                          cGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGt                                                                  >  1:3415677/1‑91 (MQ=255)
                          cGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGt                                                                  <  1:336129/91‑1 (MQ=255)
                           gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                 <  2:2326647/91‑1 (MQ=255)
                            aCCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGttt                                                                <  1:2016666/91‑1 (MQ=255)
                              cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTAt                                                              <  2:3089752/91‑1 (MQ=255)
                               aTTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATc                                                             >  2:2799529/1‑91 (MQ=255)
                                   aaaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                         <  2:384649/91‑1 (MQ=255)
                                        atGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattg                                                    <  2:2046546/91‑1 (MQ=255)
                                        atGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattg                                                    >  2:1030501/1‑91 (MQ=255)
                                         tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  1:416891/1‑91 (MQ=255)
                                         tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  2:2564640/1‑91 (MQ=255)
                                           tgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATg                                                 >  2:2029979/1‑91 (MQ=255)
                                            gAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGt                                                >  2:3353247/1‑91 (MQ=255)
                                            gAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGt                                                <  2:1786490/91‑1 (MQ=255)
                                                 tAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataa                                           >  2:630601/1‑91 (MQ=255)
                                                   gACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataata                                         <  1:719897/91‑1 (MQ=255)
                                                    aCGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGCTTGATGTTataataa                                        >  2:1624505/1‑91 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTATTTTATCCAATGATTGATGTTataataat                                       >  2:857702/1‑91 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       <  1:3182170/91‑1 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       >  1:2101824/1‑91 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       <  1:2100148/91‑1 (MQ=255)
                                                       tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                     <  2:3739622/91‑1 (MQ=255)
                                                       tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                     <  1:2770578/91‑1 (MQ=255)
                                                       tACATAAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGAGGTTATAATAATAc                                     <  2:822309/91‑1 (MQ=255)
                                                         cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                   <  2:2358546/91‑1 (MQ=255)
                                                           tCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTaaa                                 >  1:1144125/1‑91 (MQ=255)
                                                             atgtACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                               >  1:817017/4‑91 (MQ=255)
                                                             aaaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                               >  1:1109362/1‑91 (MQ=255)
                                                                 aCGTCAAGGCATGTCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTCGTa                           >  1:3003113/1‑91 (MQ=255)
                                                                 aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                           <  1:1491460/91‑1 (MQ=255)
                                                                 aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                           >  1:2939341/1‑91 (MQ=255)
                                                                  cGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTAt                          >  2:3322032/1‑91 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         >  2:18871/1‑91 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         >  1:882584/1‑91 (MQ=255)
                                                                    tCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCt                        <  2:431112/91‑1 (MQ=255)
                                                                     cAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCta                       <  2:3466190/91‑1 (MQ=255)
                                                                      aaGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtat                      <  2:1519291/91‑1 (MQ=255)
                                                                       aGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtata                     >  1:859736/1‑91 (MQ=255)
                                                                        ggCATATCCTAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATAGATGTTATAATAAGACTAAATTTGTATCTACaa                    <  2:1585262/91‑1 (MQ=255)
                                                                        ggCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaa                    >  1:1784159/1‑91 (MQ=255)
                                                                          cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                  <  1:2417958/91‑1 (MQ=255)
                                                                          cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                  <  2:128874/91‑1 (MQ=255)
                                                                            tatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaaa                >  2:1378552/1‑91 (MQ=255)
                                                                              tCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGt              <  2:724867/91‑1 (MQ=255)
                                                                               ccAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTa             >  2:2622214/1‑91 (MQ=255)
                                                                               ccAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTa             >  2:1018502/1‑91 (MQ=255)
                                                                                cAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTaa            >  2:897749/1‑91 (MQ=255)
                                                                                cAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTaa            <  2:1747474/91‑1 (MQ=255)
                                                                                cAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTaa            <  2:2245186/91‑1 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          >  2:3271881/1‑91 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          >  1:2710798/1‑91 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  1:3706418/91‑1 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  1:3439123/91‑1 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  1:1747884/91‑1 (MQ=255)
                                                                                   aTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGa         <  1:230876/91‑1 (MQ=255)
                                                                                    tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        >  1:3008053/1‑91 (MQ=255)
                                                                                    tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAGGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        <  2:291658/91‑1 (MQ=255)
                                                                                    aaaaaCACACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTCTATAAAAGTAATGAg        >  2:1168777/3‑91 (MQ=37)
                                                                                      aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa      >  1:2825051/1‑91 (MQ=255)
                                                                                       aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAt     >  1:227540/1‑91 (MQ=255)
                                                                                       aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAt     >  1:1454500/1‑91 (MQ=255)
                                                                                        aCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt    >  1:773026/1‑91 (MQ=255)
                                                                                         cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt   <  2:149814/91‑1 (MQ=255)
                                                                                          gCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTTTATCTATAAAAAAGTAATGAGCATTtg  >  2:556097/1‑91 (MQ=255)
                                                                                          gCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtg  >  1:2216650/1‑91 (MQ=255)
                                                                                          gCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTtg  >  2:303742/1‑91 (MQ=255)
                                                                                          |                                                                                          
GTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221563‑1221743

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: