Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,156,333 G→A 6.2% intergenic (+22/‑302) trxA_2 → / → uvrC Thioredoxin/UvrABC system protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,156,3330GA 6.2% 831.6 / 26.2 258intergenic (+22/‑302)trxA_2/uvrCThioredoxin/UvrABC system protein C
Reads supporting (aligned to +/- strand):  ref base G (134/108);  new base A (6/10);  total (140/118)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 271 total aligned reads displayed.

CCAGTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTGAGA  >  USA300TCH1516_ALE/1156243‑1156423
                                                                                          |                                                                                          
ccAGTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTggg                                                                                            <  2:2451433/91‑1 (MQ=255)
ccAGTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTggg                                                                                            <  1:2839527/91‑1 (MQ=255)
 cAGTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTgggg                                                                                           >  2:488381/1‑91 (MQ=255)
   gTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCa                                                                                         >  2:2967645/1‑91 (MQ=255)
       aTAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTc                                                                                     >  2:905448/1‑91 (MQ=255)
       aTAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTc                                                                                     <  1:291547/91‑1 (MQ=255)
        tAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCt                                                                                    >  2:3488252/1‑91 (MQ=255)
        tAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCt                                                                                    >  1:735090/1‑91 (MQ=255)
        tAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCt                                                                                    <  1:2967645/91‑1 (MQ=255)
          aaGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGACGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCttt                                                                                  >  1:1214569/1‑91 (MQ=255)
          aaGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCttt                                                                                  <  2:390487/91‑1 (MQ=255)
          aaGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCttt                                                                                  >  1:2736651/1‑91 (MQ=255)
           aGTTGTTAGTATCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTa                                                                                 <  2:1442305/91‑1 (MQ=255)
            gttgttGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTaa                                                                                >  1:105105/1‑91 (MQ=255)
            gttgttGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTACAAGTTACAACCAATGACGACTGGGGCATTTCTTTaa                                                                                >  2:2669254/1‑91 (MQ=255)
             ttgttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAAt                                                                               <  2:1743811/91‑1 (MQ=255)
              tgttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATg                                                                              <  2:407925/91‑1 (MQ=255)
              tgttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCAATTCATTaaag                                                                              <  1:1903319/91‑3 (MQ=255)
               gttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGa                                                                             >  2:3661514/1‑91 (MQ=255)
               gttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACAAATGACGACTGGGGCATTTCTTTAATGa                                                                             <  2:312083/91‑1 (MQ=255)
                ttgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGaa                                                                            >  1:3244768/1‑91 (MQ=255)
                 tgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTATAACCAATGACGACTGGGGCATTTCTTTAATGAAt                                                                           >  1:1000155/1‑91 (MQ=255)
                 tgGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAAt                                                                           <  1:2424331/91‑1 (MQ=255)
                  ggTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAAtt                                                                          >  1:2242553/1‑91 (MQ=255)
                     ttCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCt                                                                       >  2:2193587/1‑91 (MQ=255)
                       gaaaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTcc                                                                     >  1:1617407/3‑91 (MQ=255)
                       gaaaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTcc                                                                     >  1:620469/3‑91 (MQ=255)
                       gaaaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTcc                                                                     >  2:2909970/3‑91 (MQ=255)
                       gaaaccaaAAGAAAACGTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGATcc                                                                     >  2:2344283/3‑91 (MQ=255)
                        caaccaaAAGAAAACTTAGCTGAAGTTTTGGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCa                                                                    <  1:2170055/91‑1 (MQ=255)
                        caaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGAGCAATTCATTAAAGAAATGCCCCa                                                                    <  1:1071942/91‑1 (MQ=25)
                         aaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAg                                                                   <  1:739985/91‑1 (MQ=255)
                         aaccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGAGCAATTCATTAAAGAAATGCCCCAg                                                                   <  2:2231638/91‑1 (MQ=25)
                           ccaaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGtt                                                                 >  2:3237094/1‑91 (MQ=255)
                            caaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGttt                                                                >  2:1703834/1‑91 (MQ=255)
                            caaAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGttt                                                                >  2:1584182/1‑91 (MQ=255)
                             aaaaGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGtttt                                                               <  1:1434114/91‑1 (MQ=255)
                             aaaaGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGtttt                                                               >  2:1170630/1‑91 (MQ=255)
                             aaaaGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGtttt                                                               >  1:540442/1‑91 (MQ=255)
                                 gAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTtttgtt                                                           >  2:2488241/1‑91 (MQ=255)
                                 gAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTtttgtt                                                           <  2:1198997/91‑1 (MQ=255)
                                   aaaCTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTtttgtttg                                                         <  2:1611850/91‑1 (MQ=255)
                                    aaCTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTttgtttgt                                                        >  1:3197262/1‑91 (MQ=255)
                                     aCTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTtttgttt                                                          <  1:842766/88‑1 (MQ=255)
                                     aCTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTtttgttt                                                          >  2:842783/1‑88 (MQ=255)
                                      cTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAGTGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTgtgt                                                      <  2:3403346/91‑1 (MQ=255)
                                      cTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTgtgt                                                      <  1:2881149/91‑1 (MQ=255)
                                         aGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGtttt                                                   >  1:2529899/1‑91 (MQ=255)
                                          gCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGttttt                                                  >  1:2603375/1‑91 (MQ=255)
                                          gCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGttttt                                                  <  2:812895/91‑1 (MQ=255)
                                            tGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTaa                                                <  1:186032/91‑1 (MQ=255)
                                            tGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTaa                                                <  1:3376753/91‑1 (MQ=255)
                                             gAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAat                                               <  1:1584182/91‑1 (MQ=255)
                                             gAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAat                                               >  1:2944950/1‑91 (MQ=255)
                                              aaGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAata                                              <  2:2079495/91‑1 (MQ=255)
                                              aaGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAata                                              <  1:135487/91‑1 (MQ=255)
                                                gTTTTAGATAAACATTTATAAGTTACAACCAATGGCGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAAtata                                            <  2:3320976/91‑1 (MQ=255)
                                                 ttttAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATaa                                           <  1:1968866/91‑1 (MQ=255)
                                                   ttAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATaaaa                                         >  2:2797261/1‑91 (MQ=255)
                                                   ttAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATaaaa                                         <  1:488381/91‑1 (MQ=255)
                                                    tAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGCCTTGCTCCAGTTTTTGTTCGTGTTTTTAATCTacaaa                                        >  1:1671571/1‑91 (MQ=255)
                                                    tAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATaaaaa                                        >  1:1489826/1‑91 (MQ=255)
                                                     aGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAg                                       >  2:1324260/1‑91 (MQ=255)
                                                     aGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAg                                       <  2:2399008/91‑1 (MQ=255)
                                                        tAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTg                                    <  2:540442/91‑1 (MQ=255)
                                                               ttATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATa                             >  1:2457775/1‑91 (MQ=255)
                                                                tataAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATaa                            >  1:385862/1‑91 (MQ=255)
                                                                 ataAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAg                           >  1:3401384/1‑91 (MQ=255)
                                                                  taAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGt                          <  2:2860784/91‑1 (MQ=255)
                                                                   aaGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGtc                         >  2:3544091/1‑91 (MQ=255)
                                                                    aGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGGTAAGtca                        <  1:1613824/91‑1 (MQ=255)
                                                                    aGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGtca                        <  2:2050318/91‑1 (MQ=255)
                                                                    aGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGtca                        <  1:345404/91‑1 (MQ=255)
                                                                     gTTACAACCGATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGtcat                       <  1:3716073/91‑1 (MQ=255)
                                                                      ttACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGtcatc                      >  2:2033629/1‑91 (MQ=255)
                                                                          aaCCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCatat                  >  2:793173/1‑91 (MQ=255)
                                                                          aaCCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCatat                  >  2:3455054/1‑91 (MQ=255)
                                                                           aCCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATAtt                 >  2:3835966/1‑91 (MQ=255)
                                                                            ccAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTg                >  1:993179/1‑91 (MQ=255)
                                                                             cAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGt               <  2:2611620/91‑1 (MQ=255)
                                                                             cAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGt               <  1:3195291/91‑1 (MQ=255)
                                                                               aTGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGttt             <  1:2358229/91‑1 (MQ=255)
                                                                                 gacgacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTAc           >  1:262259/1‑91 (MQ=255)
                                                                                   cgacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGa         >  1:888115/1‑91 (MQ=255)
                                                                                    gacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGAc        <  2:270770/91‑1 (MQ=255)
                                                                                    gacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGAc        <  2:3738524/91‑1 (MQ=255)
                                                                                    gacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGAc        <  2:186858/91‑1 (MQ=255)
                                                                                    gacTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGAc        >  2:2841426/1‑91 (MQ=255)
                                                                                    aacTGGAGCAATTCATTAAAGAAATGCCCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGAc        >  1:1151995/2‑91 (MQ=25)
                                                                                     acTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACt       >  1:909835/1‑91 (MQ=255)
                                                                                     acTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACt       >  2:2282200/1‑91 (MQ=255)
                                                                                      cTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACtt      >  2:1753646/1‑91 (MQ=255)
                                                                                      cTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACtt      <  2:3870407/91‑1 (MQ=255)
                                                                                      cTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACtt      >  1:3591889/1‑91 (MQ=255)
                                                                                      cTGGAGCAATTCATTAAAGAAATGCCCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACtt      >  2:338521/1‑91 (MQ=25)
                                                                                        ggggCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTga    <  1:915060/91‑1 (MQ=255)
                                                                                        ggggCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTga    <  2:921092/91‑1 (MQ=255)
                                                                                          ggCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTgaga  >  1:3338243/1‑91 (MQ=255)
                                                                                          ggCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTgaga  <  1:3048238/91‑1 (MQ=255)
                                                                                          ggCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTgaga  <  1:2402750/91‑1 (MQ=255)
                                                                                          |                                                                                          
CCAGTTGATAAAGTTGTTGGTTTCCAACCAAAAGAAAACTTAGCTGAAGTTTTAGATAAACATTTATAAGTTACAACCAATGACGACTGGGGCATTTCTTTAATGAATTGCTCCAGTTTTTGTTTGTGTTTTTAATATAAAAAGTTGAATGATAAGTCATCATATTGTTTACGACTTGAGA  >  USA300TCH1516_ALE/1156243‑1156423

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: