Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,221,661 G→A 8.8% intergenic (+76/‑140) 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,6610GA 8.8% 655.6 / 33.6 212intergenic (+76/‑140)rpoZ/coaBCDNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC
Reads supporting (aligned to +/- strand):  ref base G (98/95);  new base A (7/12);  total (105/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 = 1.00e+00

Only 100 of 217 total aligned reads displayed.

TATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATA  >  USA300TCH1516_ALE/1221571‑1221751
                                                                                          |                                                                                          
tatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                            >  2:3326930/1‑91 (MQ=255)
tatGGTAAAGAATAGTTTCGACCATTAAAAACATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                            >  1:3509017/1‑91 (MQ=255)
 atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                           >  2:2291441/1‑91 (MQ=255)
 atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                           >  2:3962909/1‑91 (MQ=255)
      aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtg                                                                                      <  1:1126152/91‑1 (MQ=255)
      aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtg                                                                                      <  1:1117713/91‑1 (MQ=255)
      aaaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtg                                                                                      <  2:1213693/91‑1 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                                     >  2:2986018/1‑91 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                                     <  2:493716/91‑1 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                     <  1:1176286/91‑2 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                     <  1:1707471/91‑2 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                     <  1:991931/91‑2 (MQ=255)
       aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                                     <  1:2095934/91‑2 (MQ=255)
        aGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                                    <  2:1029043/91‑1 (MQ=255)
         gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                                   <  1:3598009/91‑1 (MQ=255)
           aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                                 >  2:3573719/1‑91 (MQ=255)
           aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                                 >  1:664053/1‑91 (MQ=255)
             aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCGAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttttt                                                                               >  2:3245825/1‑91 (MQ=255)
              gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                              >  2:3042927/1‑91 (MQ=255)
              gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                              <  1:867888/91‑1 (MQ=255)
              gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtctttt                                                                              <  1:3125433/91‑1 (MQ=255)
                ttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTa                                                                            >  2:2814700/1‑91 (MQ=255)
                 tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                           <  1:1861033/91‑1 (MQ=255)
                 tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                           <  2:2937727/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:2864973/1‑91 (MQ=255)
                   gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                         <  2:1083006/91‑1 (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)
                            aaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATTCAAtg                                                                >  2:110784/1‑91 (MQ=255)
                             aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                               >  2:877278/1‑91 (MQ=255)
                             aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                               <  1:3275354/91‑1 (MQ=255)
                             aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                               <  2:2209528/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:2002768/91‑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:1885469/1‑91 (MQ=255)
                                             cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATGGATGTTataataat                                               >  1:3909882/1‑91 (MQ=255)
                                               tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                             <  1:1787818/91‑1 (MQ=255)
                                               tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                             <  2:504032/91‑1 (MQ=255)
                                                aCATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACt                                            <  2:1010398/91‑1 (MQ=255)
                                                 cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                           <  1:3392889/91‑1 (MQ=255)
                                                 cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                           >  2:722341/1‑91 (MQ=255)
                                                 cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                           >  1:3814206/1‑91 (MQ=255)
                                                   tCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTaaa                                         >  2:732492/1‑91 (MQ=255)
                                                    gatgtACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAt                                        >  1:462826/5‑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:3767025/1‑91 (MQ=255)
                                                       aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                                     >  2:2202940/1‑91 (MQ=255)
                                                         aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                   >  1:1144063/1‑91 (MQ=255)
                                                         aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                   >  1:2181748/1‑91 (MQ=255)
                                                         aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                                   >  2:756561/1‑91 (MQ=255)
                                                          cGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTAt                                  <  1:3090554/91‑1 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 <  1:1511814/91‑1 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 <  2:2554701/91‑1 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 <  2:2587766/91‑1 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 >  2:2256365/1‑91 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 <  1:1510728/91‑1 (MQ=255)
                                                           gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                                 <  1:3022370/91‑1 (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                          >  1:308048/1‑91 (MQ=255)
                                                                  cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                          >  2:466571/1‑91 (MQ=255)
                                                                   atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                         <  1:2823199/91‑1 (MQ=255)
                                                                    tatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaaa                        >  2:1973611/1‑91 (MQ=255)
                                                                          aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg                  <  1:3512303/91‑1 (MQ=255)
                                                                          aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg                  >  1:673709/1‑91 (MQ=255)
                                                                            tCAAACGCGCATATGTAGGTGTGTTGTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg                <  2:1656622/91‑1 (MQ=255)
                                                                            tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg                >  2:780416/1‑91 (MQ=255)
                                                                            tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg                <  2:156876/91‑1 (MQ=255)
                                                                             cAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGc               <  1:2541283/91‑1 (MQ=255)
                                                                              aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa              <  2:1596121/91‑1 (MQ=255)
                                                                                acacACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt            >  1:2971576/4‑91 (MQ=255)
                                                                                 cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt           >  1:264745/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)
                                                                                          |                                                                                          
TATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATA  >  USA300TCH1516_ALE/1221571‑1221751

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: