Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,275,552 C→T 5.4% T302I (ACA→ATA)  mtlR → Transcriptional regulator MtlR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,275,5520CT 5.4% 677.8 / 15.8 203T302I (ACA→ATA) mtlRTranscriptional regulator MtlR
Reads supporting (aligned to +/- strand):  ref base C (83/109);  new base T (3/8);  total (86/117)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 225 total aligned reads displayed.

GTCAAATATATGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTGGA  >  USA300TCH1516_ALE/2275462‑2275642
                                                                                          |                                                                                          
gTCAAATATATGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAc                                                                                            <  2:2526285/91‑1 (MQ=255)
 tCAAATATATGACGTCACGTTTAATCAGGCAGAGATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAca                                                                                           <  1:3184793/91‑1 (MQ=255)
 tCAAATATATGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAca                                                                                           <  1:2850638/91‑1 (MQ=255)
  cAAATATATGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAcat                                                                                          <  1:1020873/91‑1 (MQ=255)
     atatatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAcatcat                                                                                       >  1:2528213/1‑91 (MQ=255)
     atatatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAcatcat                                                                                       >  1:2526530/1‑91 (MQ=255)
     atatatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATAcatcat                                                                                       <  2:2017188/91‑1 (MQ=255)
      tatatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATGTATCAtt                                                                                      <  1:3501985/91‑1 (MQ=255)
        tatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTaa                                                                                    <  1:2735394/91‑1 (MQ=255)
        tatGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTaa                                                                                    <  1:573576/91‑1 (MQ=255)
          tGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATGTATCATTaaga                                                                                  <  1:1353203/91‑3 (MQ=255)
          tGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaat                                                                                  >  2:3621432/1‑91 (MQ=255)
          tGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaat                                                                                  <  1:2736260/91‑1 (MQ=255)
          tGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaat                                                                                  <  1:2818682/91‑1 (MQ=255)
           gACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTGCTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaata                                                                                 <  2:3133242/91‑1 (MQ=255)
            aCGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaataa                                                                                >  2:3435921/1‑91 (MQ=255)
             cGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAaataat                                                                               <  2:3790913/91‑1 (MQ=255)
             cGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGGGGAAATAAAAGAAAAAATCTTAATGATACATCATTAaataat                                                                               <  2:3790502/91‑1 (MQ=255)
              gTCCCGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATc                                                                              >  2:3039923/1‑91 (MQ=255)
              gTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATc                                                                              >  2:3685410/1‑91 (MQ=255)
               tCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCg                                                                             >  2:3093860/1‑91 (MQ=255)
                  cGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTtg                                                                          <  1:3339152/91‑1 (MQ=255)
                   gTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTtgt                                                                         >  2:2198455/1‑91 (MQ=255)
                     ttAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTga                                                                       >  1:2321353/1‑91 (MQ=255)
                     ttAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTga                                                                       >  2:2539732/1‑91 (MQ=255)
                     ttAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTga                                                                       <  1:1813052/91‑1 (MQ=255)
                      tAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTgaa                                                                      <  1:2128576/91‑1 (MQ=255)
                        aTCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTgaaga                                                                    <  2:553610/91‑1 (MQ=255)
                        aTCAGGCAGAAATTGCTTTCATTACTAACCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTgaaga                                                                    <  2:553476/91‑1 (MQ=255)
                         tCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTgaagaa                                                                   <  2:1732589/91‑1 (MQ=255)
                         tCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTgaagaa                                                                   <  1:1116726/91‑1 (MQ=255)
                           aGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGaaaa                                                                 >  1:2923383/1‑91 (MQ=255)
                            ggCAGAAATTGCTTTCATTACTATCCATTCACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAAc                                                                >  2:3377742/1‑91 (MQ=255)
                              cAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACaa                                                              <  2:1259549/91‑1 (MQ=255)
                                 aaaTTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAAt                                                           <  2:476907/91‑1 (MQ=255)
                                  aaTTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAAtt                                                          >  2:823025/1‑91 (MQ=255)
                                  aaTTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAAtt                                                          >  2:2941720/1‑91 (MQ=255)
                                  aaTTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAAtt                                                          >  1:3652078/1‑91 (MQ=255)
                                   aTTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTa                                                         >  2:1393442/1‑91 (MQ=255)
                                    ttGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTaa                                                        >  2:2642501/1‑91 (MQ=255)
                                     tGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTaaa                                                       <  2:2742270/91‑1 (MQ=255)
                                     tGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTaaa                                                       >  1:183110/1‑91 (MQ=255)
                                     tGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTaaa                                                       >  1:184569/1‑91 (MQ=255)
                                     tGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAAACTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTaaa                                                       <  1:1758903/91‑1 (MQ=255)
                                      gCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAg                                                      <  1:1530090/91‑1 (MQ=255)
                                      gCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAg                                                      >  1:2777504/1‑91 (MQ=255)
                                      gCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAg                                                      <  2:1108032/91‑1 (MQ=255)
                                      gCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAg                                                      >  2:1248358/1‑91 (MQ=255)
                                         ttCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGt                                                   >  2:2366237/1‑91 (MQ=255)
                                          tCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGtt                                                  >  1:804051/1‑91 (MQ=255)
                                             ttACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGt                                               >  2:3673159/1‑91 (MQ=255)
                                                cTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTaa                                            >  1:1258299/1‑91 (MQ=255)
                                                cTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTaa                                            <  1:2481631/91‑1 (MQ=255)
                                                     cATTTACGTGGAGATAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAg                                       <  2:2099156/91‑1 (MQ=255)
                                                         tACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGa                                   >  1:699953/1‑91 (MQ=255)
                                                          aCGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGaa                                  <  2:2732284/91‑1 (MQ=255)
                                                           cGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGaaa                                 <  2:1382526/91‑1 (MQ=255)
                                                           cGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGaaa                                 <  2:188283/91‑1 (MQ=255)
                                                            gTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAAt                                <  2:184931/91‑1 (MQ=255)
                                                             tgaaGCTATACGAAAGAATCTTAATGATACATCATTAAATAAACGTTGTGGAGAAAACAAAATGAAAGCGTTTGTTATCAAAGGAGAAATg                               <  2:3515292/88‑1 (MQ=255)
                                                             tGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATg                               <  1:3039799/91‑1 (MQ=255)
                                                              ggAGCTAAACGAATAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGa                              <  1:3882092/91‑1 (MQ=255)
                                                               gAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGAt                             <  2:3218114/91‑1 (MQ=255)
                                                               gAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGAt                             <  2:3408049/91‑1 (MQ=255)
                                                               gAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGAt                             <  1:3719399/91‑1 (MQ=255)
                                                                 gCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGAttt                           <  1:1273151/91‑1 (MQ=255)
                                                                   tAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTcc                         <  1:1248298/91‑1 (MQ=255)
                                                                     aaCGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCgg                       <  1:9370/91‑1 (MQ=255)
                                                                      aCGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGt                      <  1:359327/91‑1 (MQ=255)
                                                                       cGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTa                     >  1:3746847/1‑91 (MQ=255)
                                                                         taaaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATg                   <  2:3839953/90‑1 (MQ=255)
                                                                         aaaaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATg                   >  2:725013/1‑91 (MQ=255)
                                                                         aaaaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATg                   <  1:1183242/91‑1 (MQ=255)
                                                                         aaaaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATg                   <  2:1452082/91‑1 (MQ=255)
                                                                           aaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGAc                 >  2:756861/1‑91 (MQ=255)
                                                                           aaaaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGAc                 <  2:1506065/91‑1 (MQ=255)
                                                                             aaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAt               >  2:3265931/1‑91 (MQ=255)
                                                                             aaTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAt               >  2:3740016/1‑91 (MQ=255)
                                                                              aTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAtt              <  2:2499903/91‑1 (MQ=255)
                                                                              aTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAtt              >  2:2845556/1‑91 (MQ=255)
                                                                              aTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAtt              >  1:3097998/1‑91 (MQ=255)
                                                                              aTCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAtt              <  1:1015533/91‑1 (MQ=255)
                                                                               tCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACAttt             <  2:3468288/91‑1 (MQ=255)
                                                                                tttAATGATGTATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTg            >  1:1164601/2‑91 (MQ=255)
                                                                                cTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTg            <  2:1614068/91‑1 (MQ=255)
                                                                                 ttAATGATGTATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGc           >  2:2285883/1‑91 (MQ=255)
                                                                                 ttAATGATGTATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGc           >  2:3859326/1‑91 (MQ=255)
                                                                                 ttAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGc           >  1:2647494/1‑91 (MQ=255)
                                                                                  tAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCa          <  2:332334/91‑1 (MQ=255)
                                                                                   aaTGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAg         >  2:404177/1‑91 (MQ=255)
                                                                                   aaTGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTAGCAg         >  1:3234875/1‑91 (MQ=255)
                                                                                    aTGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGa        <  2:1000191/91‑1 (MQ=255)
                                                                                      gATACATCATTAAATAATCGTTGTGAATAAAACAAAATTAAACCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGAtt      <  1:1063605/91‑1 (MQ=255)
                                                                                      gATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGAtt      <  1:1922255/91‑1 (MQ=255)
                                                                                      gATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGAtt      <  1:2174829/91‑1 (MQ=255)
                                                                                      gATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGAtt      <  2:1239025/91‑1 (MQ=255)
                                                                                        tACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTg    <  2:589183/91‑1 (MQ=255)
                                                                                        tACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTg    <  1:3579990/91‑1 (MQ=255)
                                                                                          catcatTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTGGa  >  1:211502/1‑91 (MQ=255)
                                                                                          catcatTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTGGa  <  1:556130/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTCAAATATATGACGTCACGTTTAATCAGGCAGAAATTGCTTTCATTACTATCCATTTACGTGGAGCTAAACGAAAAAATCTTAATGATACATCATTAAATAATCGTTGTGAAGAAAACAAAATTAAAGCGTTTGTTAACAAAGTAGAAATGATTTCCGGTATGACATTTGCAGATTTGGA  >  USA300TCH1516_ALE/2275462‑2275642

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: