Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,479,047 A→G 9.6% intergenic (+3/+137) USA300HOU_RS12765 → / ← hrtA_2 hypothetical protein/Putative hemin import ATP‑binding protein HrtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,479,0470AG 9.6% 560.7 / 38.4 187intergenic (+3/+137)USA300HOU_RS12765/hrtA_2hypothetical protein/Putative hemin import ATP‑binding protein HrtA
Reads supporting (aligned to +/- strand):  ref base A (81/88);  new base G (10/8);  total (91/96)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 216 total aligned reads displayed.

TGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACTTTTA  >  USA300TCH1516_ALE/2478957‑2479137
                                                                                          |                                                                                          
tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                            >  2:468294/1‑91 (MQ=255)
tGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaatta                                                                                            >  2:1197126/1‑91 (MQ=255)
 ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                           <  2:2768511/91‑1 (MQ=255)
 ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                           <  1:1803343/91‑1 (MQ=255)
 ggTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaa                                                                                           <  1:17918/91‑1 (MQ=255)
  gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          >  1:260695/1‑91 (MQ=255)
  gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          <  1:1068952/91‑1 (MQ=255)
  gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          >  1:3057986/1‑91 (MQ=255)
  gTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          <  2:1418935/91‑1 (MQ=255)
      tAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGaa                                                                                      >  1:2935991/1‑91 (MQ=255)
        aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          <  2:842488/85‑1 (MQ=255)
        aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAAtaattaat                                                                                          >  1:842488/1‑85 (MQ=255)
        aaCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAAt                                                                                    >  1:1060834/1‑91 (MQ=255)
         aCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                                   <  2:122926/91‑1 (MQ=255)
         aCGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATc                                                                                   <  2:1060834/91‑1 (MQ=255)
          cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCa                                                                                  <  2:1705964/91‑1 (MQ=18)
          cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                                  >  2:1609794/1‑91 (MQ=255)
          cGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCa                                                                                  >  2:424681/1‑91 (MQ=255)
           gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                                 >  1:1853848/1‑91 (MQ=255)
           gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                                 >  2:1774074/1‑91 (MQ=255)
           gTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAg                                                                                 <  1:2526831/91‑1 (MQ=255)
            tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                                >  1:1483230/1‑91 (MQ=255)
            tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                                >  1:746971/1‑91 (MQ=255)
            tGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGc                                                                                >  2:767315/1‑91 (MQ=255)
              aCTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtt                                                                              <  1:3862655/91‑1 (MQ=255)
               cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCttt                                                                             <  1:2879575/91‑1 (MQ=17)
               cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATTAATTATTTGATAGAAATCAGCttt                                                                             >  2:1024368/1‑91 (MQ=17)
               cTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttt                                                                             <  1:3760924/91‑1 (MQ=255)
                 tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                           <  1:2711356/91‑1 (MQ=255)
                 tCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCttttt                                                                           <  1:136110/91‑1 (MQ=255)
                  cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                          <  2:3720869/91‑1 (MQ=255)
                  cAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCtttttt                                                                          >  2:2336848/1‑91 (MQ=255)
                    aGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTAc                                                                        <  2:1320562/91‑1 (MQ=255)
                     gCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACa                                                                       >  2:2410939/1‑91 (MQ=255)
                        atatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTg                                                                    <  2:1527770/91‑1 (MQ=40)
                        atatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTg                                                                    <  1:3880213/91‑1 (MQ=40)
                         tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                                   >  1:2462682/1‑91 (MQ=39)
                         tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                                   <  1:562312/91‑1 (MQ=39)
                         tatGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGc                                                                   >  2:1934909/1‑91 (MQ=39)
                           tGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCt                                                                 <  1:2578519/91‑1 (MQ=37)
                            gAATGCGACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                                >  2:1844795/1‑91 (MQ=18)
                            gAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTa                                                                <  1:367659/91‑1 (MQ=35)
                             aaTGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTaa                                                               <  2:2829875/91‑1 (MQ=34)
                               tGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGa                                                             <  1:3222661/91‑1 (MQ=32)
                                              atgataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTGTACATTGCCTAAGAACTTAATGTCCCAGc                                              <  2:292131/88‑1 (MQ=12)
                                              aTTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                              <  2:2111535/91‑1 (MQ=32)
                                              aTTATAAAATACATAACAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGc                                              <  2:1136334/91‑1 (MQ=17)
                                               ttATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGcc                                             <  1:116819/91‑1 (MQ=33)
                                                tataAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCt                                            <  1:1789356/91‑1 (MQ=34)
                                                  taAAATATATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGt                                          <  1:1938756/91‑1 (MQ=37)
                                                  taAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGt                                          <  1:1914788/91‑1 (MQ=37)
                                                   aaaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtt                                         <  1:18715/91‑1 (MQ=38)
                                                    aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                                        >  1:1419992/1‑91 (MQ=39)
                                                    aaaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttt                                        <  1:1885215/91‑1 (MQ=39)
                                                     aaTACATATCAAAAAAGCTGATTTGTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                       <  2:3606546/91‑1 (MQ=255)
                                                     aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                       <  2:3198461/91‑1 (MQ=40)
                                                     aaTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGtttt                                       <  2:3864051/91‑1 (MQ=39)
                                                      aTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttttt                                      <  1:2330018/91‑1 (MQ=255)
                                                      aTACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGttttt                                      <  2:3344124/91‑1 (MQ=255)
                                                        aCATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTg                                    >  2:621491/1‑91 (MQ=255)
                                                          atatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTa                                  >  2:3009654/1‑91 (MQ=255)
                                                           tatCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaa                                 <  1:1272158/91‑1 (MQ=255)
                                                            atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                                <  2:1458381/91‑1 (MQ=255)
                                                            atCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTaaa                                >  1:845779/1‑91 (MQ=255)
                                                             tCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAc                               >  2:2726906/1‑91 (MQ=255)
                                                               aaaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAAccc                             >  2:1871591/1‑91 (MQ=255)
                                                                aaaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCa                            <  1:3517478/91‑1 (MQ=255)
                                                                 aaaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaa                           >  2:2497034/1‑91 (MQ=255)
                                                                  aaaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCaaa                          >  1:3260182/1‑91 (MQ=255)
                                                                   aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAAt                         >  1:3794521/1‑91 (MQ=255)
                                                                   aaGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAAt                         <  2:1101366/91‑1 (MQ=255)
                                                                    aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                        >  1:2117872/1‑91 (MQ=255)
                                                                    aGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATc                        <  2:3057986/91‑1 (MQ=255)
                                                                     gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                       <  1:1329139/91‑1 (MQ=255)
                                                                     gCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                       <  1:1013434/91‑1 (MQ=255)
                                                                     gCTGATTTCTATCAAATAATTAATAGAAATCAGATTTTTTACATTGCCTAAGAACTTAAGGTCCCAGCCTGTTTTTTGTAAACCCAAATCg                       <  1:3444983/91‑1 (MQ=255)
                                                                        gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                    >  1:2195508/1‑91 (MQ=255)
                                                                        gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                    <  2:785046/91‑1 (MQ=255)
                                                                        gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                    >  2:351461/1‑91 (MQ=255)
                                                                        gATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAt                    >  1:2079682/1‑91 (MQ=255)
                                                                         aTTTCTATTAATTATTTGATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt                   >  2:3247919/1‑91 (MQ=21)
                                                                         aTTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTAtt                   <  2:1697955/91‑1 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCTAGCCTGTTTTTTGTAAACTCAAATCGTATTc                  >  2:3653000/1‑91 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc                  >  2:2559874/1‑91 (MQ=255)
                                                                          tttCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTCAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc                  >  2:918040/1‑91 (MQ=255)
                                                                            tCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTg                <  2:675889/91‑1 (MQ=255)
                                                                             cTATCAAATAATTAATAGAAATCAGCTTTTTTGCATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGt               >  2:3816206/1‑91 (MQ=255)
                                                                              tATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGtt              <  1:2824296/91‑1 (MQ=255)
                                                                               aTCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTa             >  2:702726/1‑91 (MQ=255)
                                                                                tCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTAt            <  2:1499454/91‑1 (MQ=255)
                                                                                 cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc                  >  2:2680331/1‑84 (MQ=255)
                                                                                 cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTc                  <  1:2680331/84‑1 (MQ=255)
                                                                                 cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATc           <  2:116846/91‑1 (MQ=255)
                                                                                 cAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATc           >  1:2059429/1‑91 (MQ=255)
                                                                                  aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt          <  1:1806004/91‑1 (MQ=255)
                                                                                  aaaTAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCt          >  1:1520569/1‑91 (MQ=255)
                                                                                       attaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACtt     >  2:3138442/1‑91 (MQ=255)
                                                                                        ttaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACttt    <  1:660603/91‑1 (MQ=255)
                                                                                        ttaatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACttt    >  1:1749343/1‑91 (MQ=255)
                                                                                          aatAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACTTTta  >  1:2524641/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGGTTTTAAACGTGACTTCAAGCTATATGAATGCAACGAGCTGTACATTATAAAATACATATCAAAAAAGCTGATTTCTATCAAATAATTAATAGAAATCAGCTTTTTTACATTGCCTAAGAACTTAATGTCCCAGCCTGTTTTTTGTAAACCCAAATCGTATTCTGTTATCTAACTTTTA  >  USA300TCH1516_ALE/2478957‑2479137

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: