Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,270,460 C→G 5.2% intergenic (+211/‑206) trmFO → / → xerC_1 Methylenetetrahydrofolate‑‑tRNA‑(uracil‑5‑)‑ methyltransferase TrmFO/Tyrosine recombinase XerC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,270,4600CG 5.2% 342.3 / 8.2 96intergenic (+211/‑206)trmFO/xerC_1Methylenetetrahydrofolate‑‑tRNA‑(uracil‑5‑)‑ methyltransferase TrmFO/Tyrosine recombinase XerC
Reads supporting (aligned to +/- strand):  ref base C (39/52);  new base G (3/2);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

AATTTAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTAAC  >  USA300TCH1516_ALE/1270372‑1270549
                                                                                        |                                                                                         
aaTTTAGAATTCGCATTCAAAAAATATCTAAATAAATAACACTCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACtt                                                                                         <  1:1333645/91‑1 (MQ=255)
aaTTTAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACtt                                                                                         <  1:3030550/91‑1 (MQ=255)
   ttAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACtttat                                                                                      <  1:578529/91‑1 (MQ=255)
   ttAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACtttat                                                                                      <  1:4119482/91‑1 (MQ=255)
    tAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAAGtttatt                                                                                     <  1:2335774/91‑1 (MQ=255)
     aGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACtttattt                                                                                    >  1:421240/1‑91 (MQ=255)
       aaTTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAAGTTTATTTaa                                                                                  <  1:49910/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGGGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:3162137/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:2997567/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:2797683/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:433944/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:3525033/91‑1 (MQ=255)
          tCGCATTCAAAAAATATCTAAATAAATAACACGCAAAAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAg                                                                               <  1:2796560/91‑1 (MQ=255)
           cGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAgg                                                                              >  1:2161354/1‑91 (MQ=255)
           cGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAgg                                                                              >  1:1513347/1‑91 (MQ=255)
             cATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCa                                                                            <  1:3762487/91‑1 (MQ=255)
              aTTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGcct                                                                           <  1:2733151/91‑3 (MQ=255)
              aTTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCaa                                                                           <  1:2635792/91‑1 (MQ=255)
               ttCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAAt                                                                          <  1:601002/91‑1 (MQ=255)
                    aaaaTATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAaat                                                                     <  1:959196/91‑1 (MQ=255)
                     aaaTATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAaata                                                                    >  1:1909414/1‑91 (MQ=255)
                     aaaTATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAaata                                                                    <  1:730633/91‑1 (MQ=255)
                        tatCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAataata                                                                 >  1:2574414/1‑91 (MQ=255)
                         atCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAAtaataa                                                                <  1:1263183/91‑1 (MQ=255)
                          tCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATaaa                                                               >  1:577932/1‑91 (MQ=255)
                            taaataaaTAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATg                                                             >  1:2720620/1‑91 (MQ=255)
                            taaataaaTAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATg                                                             <  1:645973/91‑1 (MQ=255)
                             aaataaatAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATgg                                                            <  1:2930008/91‑1 (MQ=255)
                              aataaataACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGc                                                           >  1:2541850/1‑91 (MQ=255)
                                 aaataaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAg                                                        <  1:1169991/91‑1 (MQ=255)
                                 aaataaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAg                                                        <  1:1170679/91‑1 (MQ=255)
                                 aaataaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAg                                                        <  1:2999359/91‑1 (MQ=255)
                                   ataaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGc                                                      <  1:3035919/91‑1 (MQ=255)
                                   ataaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGc                                                      <  1:3582243/91‑1 (MQ=255)
                                    taaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCt                                                     <  1:1977341/91‑1 (MQ=255)
                                    taaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCt                                                     <  1:1977314/91‑1 (MQ=255)
                                     aaCACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTa                                                    <  1:2561415/91‑1 (MQ=255)
                                      acacGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTAt                                                   >  1:1820496/1‑91 (MQ=255)
                                      acacGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTAt                                                   >  1:1819980/1‑91 (MQ=255)
                                        acGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATta                                                 >  1:3162751/1‑91 (MQ=255)
                                          gCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATtaat                                               <  1:1193794/91‑1 (MQ=255)
                                          gCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATtaat                                               <  1:1193749/91‑1 (MQ=255)
                                          gCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATtaat                                               >  1:3867867/1‑91 (MQ=255)
                                            aaTAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATtaataa                                             <  1:3229314/91‑1 (MQ=255)
                                              tAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAg                                           >  1:803866/1‑91 (MQ=255)
                                                 gTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACt                                        >  1:4064959/1‑91 (MQ=255)
                                                  ttgattgaTGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACtt                                       <  1:1284556/91‑1 (MQ=255)
                                                   tgattgatGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACttt                                      >  1:140614/1‑91 (MQ=255)
                                                    gattgatGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACtttt                                     <  1:1740239/91‑1 (MQ=255)
                                                       tgatGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAGTGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGt                                  <  1:3599240/91‑1 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               <  1:2137722/91‑1 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               >  1:2699829/1‑91 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               >  1:2599860/1‑91 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               >  1:2298634/1‑91 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               <  1:2137493/91‑1 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               <  1:2137492/91‑1 (MQ=255)
                                                          tGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaat                               >  1:1189375/1‑91 (MQ=255)
                                                             aaCATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaattaa                            <  1:2992203/91‑1 (MQ=255)
                                                              aCATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaattaat                           >  1:2785465/1‑91 (MQ=255)
                                                              aCATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGtaattaat                           >  1:723143/1‑91 (MQ=255)
                                                               cATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTaattaatt                          <  1:1219861/91‑1 (MQ=255)
                                                                aTGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAAttt                         <  1:3700346/91‑1 (MQ=255)
                                                                aTGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAAttt                         <  1:3421548/91‑1 (MQ=255)
                                                                  gTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTaa                       >  1:1478203/1‑91 (MQ=255)
                                                                   tAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTaaa                      >  1:1868264/1‑91 (MQ=255)
                                                                      gagaATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTaaaaaa                   <  1:1142457/91‑1 (MQ=255)
                                                                      gagaATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTaaaaaa                   <  1:1143016/91‑1 (MQ=255)
                                                                       agaATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAag                  >  1:960268/1‑91 (MQ=255)
                                                                       agaATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAag                  >  1:1355887/1‑91 (MQ=255)
                                                                       agaATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAag                  <  1:250018/91‑1 (MQ=255)
                                                                       agaATGTTTTAAATAAACTTTATTTAAAAGGCAACGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAag                  >  1:1302653/1‑91 (MQ=255)
                                                                          aTGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAgg               <  1:1545413/91‑1 (MQ=255)
                                                                          aTGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAgg               <  1:4104245/91‑1 (MQ=255)
                                                                           tGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGt              >  1:3632175/1‑91 (MQ=255)
                                                                            gTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTa             >  1:2042642/1‑91 (MQ=255)
                                                                             ttttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAt            <  1:223167/91‑1 (MQ=255)
                                                                             ttttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAt            <  1:1372213/91‑1 (MQ=255)
                                                                             ttttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAt            >  1:3015647/1‑91 (MQ=255)
                                                                             ttttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAt            >  1:287288/1‑91 (MQ=255)
                                                                              tttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAtt           >  1:4068210/1‑91 (MQ=255)
                                                                              tttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAtt           >  1:1347611/1‑91 (MQ=255)
                                                                              tttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTAtt           <  1:1918689/91‑1 (MQ=255)
                                                                               ttAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTc          >  1:64866/1‑91 (MQ=255)
                                                                                taaataaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCt         >  1:3263030/1‑91 (MQ=255)
                                                                                 aaataaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTa        <  1:1984469/91‑1 (MQ=255)
                                                                                  aataaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTaa       <  1:3699785/91‑1 (MQ=255)
                                                                                  aataaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTaa       >  1:3696943/1‑91 (MQ=255)
                                                                                  aataaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTaa       >  1:1066210/1‑91 (MQ=255)
                                                                                    taaaGTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAAtt     >  1:3717672/1‑91 (MQ=255)
                                                                                    taaaGTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAAtt     >  1:3108818/1‑91 (MQ=255)
                                                                                    taaaGTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAAtt     >  1:3109435/1‑91 (MQ=255)
                                                                                     aaaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTa    <  1:1832978/91‑1 (MQ=255)
                                                                                      aaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTaa   >  1:2704261/1‑91 (MQ=255)
                                                                                      aaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTaa   >  1:4102963/1‑91 (MQ=255)
                                                                                      aaCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTaa   >  1:2695751/1‑91 (MQ=255)
                                                                                       aCTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTAAc  <  1:639739/91‑1 (MQ=255)
                                                                                        |                                                                                         
AATTTAGAATTCGCATTCAAAAAATATCTAAATAAATAACACGCAATAAGTTGATTGATGTAACATGTAAGAGAATGTTTTAAATAAACTTTATTTAAAAGGCAATGAAATAATAAATGGCAAGGCTATTAATAAAGACTTTTAGTAATTAATTTAAAAAAGAGGTATTCTAATTAAC  >  USA300TCH1516_ALE/1270372‑1270549

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: