Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,368,905 T→G 7.0% intergenic (+309/+71) USA300TCH1516_01277 → / ← lexA_1 hypothetical protein/LexA repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,368,9050TG 7.0% 443.7 / 19.9 129intergenic (+309/+71)USA300TCH1516_01277/lexA_1hypothetical protein/LexA repressor
Reads supporting (aligned to +/- strand):  ref base T (68/52);  new base G (2/7);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 7.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 133 total aligned reads displayed.

TATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGTACAAA  >  USA300TCH1516_ALE/1368817‑1368994
                                                                                        |                                                                                         
tatTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACttt                                                                                         >  1:3640370/1‑91 (MQ=255)
 atTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACtttt                                                                                        >  1:828929/1‑91 (MQ=255)
 atTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACtttt                                                                                        >  1:2030399/1‑91 (MQ=255)
  ttAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACttttt                                                                                       >  1:3167382/1‑91 (MQ=255)
   tAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTc                                                                                      <  1:2711839/91‑1 (MQ=255)
   tAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTc                                                                                      <  1:601476/91‑1 (MQ=39)
    aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                     <  1:3711016/91‑1 (MQ=39)
    aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                     <  1:3710260/91‑1 (MQ=39)
      aGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCAag                                                                                   <  1:53575/91‑1 (MQ=38)
         aTGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAgg                                                                                <  1:3819836/91‑1 (MQ=255)
         aTGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAgg                                                                                <  1:2151677/91‑1 (MQ=255)
          tGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCAAGAGGt                                                                               <  1:2189709/91‑1 (MQ=38)
          tGAGATTATCTAACAATTCGCAATTTCTGTGTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGt                                                                               <  1:1539200/91‑1 (MQ=255)
           gagaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTc                                                                              >  1:2341031/1‑91 (MQ=255)
           gagaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTc                                                                              >  1:4022615/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:266982/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:2548663/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:1593079/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:1674274/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:218483/1‑91 (MQ=255)
            agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                             >  1:516127/1‑91 (MQ=255)
              aTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCaa                                                                           >  1:3763943/1‑91 (MQ=255)
              aTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCaa                                                                           >  1:3763874/1‑91 (MQ=255)
                   cTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAAttt                                                                      <  1:218441/91‑1 (MQ=255)
                   cTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAAttt                                                                      >  1:2608438/1‑91 (MQ=255)
                     aaCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtat                                                                    >  1:2130909/1‑91 (MQ=255)
                      aCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtata                                                                   >  1:3663573/1‑91 (MQ=255)
                      aCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtata                                                                   >  1:2947780/1‑91 (MQ=255)
                       cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                                  <  1:25460/91‑1 (MQ=255)
                       cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                                  >  1:293896/1‑91 (MQ=255)
                       cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                                  <  1:25817/91‑1 (MQ=255)
                       cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                                  >  1:1488994/1‑91 (MQ=255)
                        aaTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaa                                                                 >  1:917700/1‑91 (MQ=255)
                         aTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaa                                                                <  1:3704842/91‑1 (MQ=255)
                         aTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaa                                                                <  1:3704805/91‑1 (MQ=255)
                          ttCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaaa                                                               <  1:2833164/91‑1 (MQ=255)
                            cGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATg                                                             <  1:3356870/91‑1 (MQ=255)
                            cGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATg                                                             <  1:2121568/91‑1 (MQ=255)
                             gCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAATTTATAAAAAagag                                                           >  1:3351076/1‑91 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           >  1:1739400/1‑91 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           >  1:700727/1‑91 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           >  1:2876103/1‑91 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           >  1:898763/1‑91 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           <  1:2361855/91‑1 (MQ=255)
                              cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                           <  1:2361748/91‑1 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:585622/1‑91 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:3627403/1‑91 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:1061998/1‑91 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:2573247/1‑91 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:1743216/1‑91 (MQ=255)
                               aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                          >  1:1742967/1‑91 (MQ=255)
                                aTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAggg                                                         >  1:1710859/1‑91 (MQ=255)
                                 tttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGc                                                        >  1:92256/1‑91 (MQ=255)
                                 tttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGc                                                        >  1:3723410/1‑91 (MQ=255)
                                 tttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGc                                                        <  1:1241386/91‑1 (MQ=255)
                                  ttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCt                                                       <  1:1111666/91‑1 (MQ=255)
                                  ttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCt                                                       <  1:4066479/91‑1 (MQ=255)
                                  ttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCCCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCt                                                       <  1:4066406/91‑1 (MQ=255)
                                   tCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCtt                                                      >  1:3478270/1‑91 (MQ=255)
                                      gTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAAc                                                   <  1:2885509/91‑1 (MQ=255)
                                       tttGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACa                                                  >  1:939530/1‑91 (MQ=255)
                                         tGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAAt                                                >  1:3234985/1‑91 (MQ=255)
                                          gAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAAtt                                               >  1:19093/1‑91 (MQ=255)
                                            aCTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTAc                                             <  1:883362/91‑1 (MQ=255)
                                             cTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACg                                            >  1:670795/1‑91 (MQ=255)
                                               gAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGat                                          <  1:2523927/91‑1 (MQ=255)
                                                aaTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGata                                         <  1:2966217/91‑1 (MQ=255)
                                                aaTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGata                                         >  1:491360/1‑91 (MQ=255)
                                                 aTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGataa                                        <  1:2899309/91‑1 (MQ=255)
                                                     agattagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAatata                                    <  1:1387788/91‑1 (MQ=255)
                                                      gattagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAatatat                                   >  1:4056482/1‑91 (MQ=255)
                                                       attagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAAtatata                                  <  1:1100596/91‑1 (MQ=255)
                                                          agaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACg                               >  1:686336/1‑91 (MQ=255)
                                                           gaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACgg                              <  1:2996652/91‑1 (MQ=255)
                                                              taaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGttt                           <  1:854126/91‑1 (MQ=255)
                                                              taaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGttt                           <  1:2766728/91‑1 (MQ=255)
                                                               aaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTa                          >  1:1692383/1‑91 (MQ=255)
                                                                aataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaa                         >  1:1502759/1‑91 (MQ=255)
                                                                aataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaa                         >  1:1279301/1‑91 (MQ=255)
                                                                 ataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaaa                        <  1:2928006/91‑1 (MQ=255)
                                                                 ataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaaa                        >  1:1866332/1‑91 (MQ=255)
                                                                   aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                      >  1:4032743/1‑91 (MQ=255)
                                                                   aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                      <  1:3524767/91‑1 (MQ=255)
                                                                   aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                      >  1:764029/1‑91 (MQ=255)
                                                                     gaCCTCTTGAAAAAGTTACGTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAatat                    >  1:1070025/2‑91 (MQ=37)
                                                                     gaCCTCTTGAAAAAGTTACGTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAatat                    >  1:1388330/2‑91 (MQ=37)
                                                                     aaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAatat                    <  1:1309822/91‑1 (MQ=255)
                                                                      aCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATAtt                   >  1:1478986/1‑91 (MQ=255)
                                                                       ccTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTa                  >  1:462050/1‑91 (MQ=255)
                                                                       ccTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTa                  <  1:942451/91‑1 (MQ=255)
                                                                          ctTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAt               <  1:90668/91‑1 (MQ=255)
                                                                           ttGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAtt              >  1:793525/1‑91 (MQ=255)
                                                                           ttGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAtt              >  1:2849178/1‑91 (MQ=255)
                                                                             gAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTc            >  1:3902178/1‑91 (MQ=255)
                                                                               aaaCGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTcgc          >  1:1770612/1‑91 (MQ=255)
                                                                               aaaCGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTcgc          <  1:1271853/91‑1 (MQ=255)
                                                                                  cGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGt       <  1:3468348/91‑1 (MQ=255)
                                                                                     aaCTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGTACa    <  1:3427191/91‑1 (MQ=255)
                                                                                      aCTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGTACaa   >  1:967444/1‑91 (MQ=255)
                                                                                       cTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGTACaaa  <  1:2727720/91‑1 (MQ=255)
                                                                                        |                                                                                         
TATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGTACAAA  >  USA300TCH1516_ALE/1368817‑1368994

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