Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 581,901 G→C 10.8% V13L (GTG→CTG)  nusG_2 → Transcription termination/antitermination protein NusG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE581,9010GC10.8% 329.9 / 25.4 102V13L (GTG→CTG) nusG_2Transcription termination/antitermination protein NusG
Reads supporting (aligned to +/- strand):  ref base G (39/52);  new base C (8/3);  total (47/55)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 105 total aligned reads displayed.

GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991
                                                                                          |                                                                                          
gATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAg                                                                                            >  1:2096348/1‑91 (MQ=255)
   ttAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCACTGc                                                                                         <  1:3823102/91‑1 (MQ=255)
   ttAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCACTGc                                                                                         <  1:2029596/91‑1 (MQ=255)
   ttAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCACTGc                                                                                         <  1:2031159/91‑1 (MQ=255)
     atgAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcat                                                                                       <  1:2800990/89‑1 (MQ=255)
      ggAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcata                                                                                      >  1:852551/1‑91 (MQ=255)
       gAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatac                                                                                     <  1:266367/91‑1 (MQ=255)
        aaTTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcataca                                                                                    <  1:2540655/91‑1 (MQ=255)
          ttACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGcatacata                                                                                  <  1:1880770/91‑1 (MQ=255)
           tACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACatat                                                                                 <  1:1364047/91‑1 (MQ=255)
           tACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACatat                                                                                 <  1:2122259/91‑1 (MQ=255)
              aGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCt                                                                              >  1:603948/1‑91 (MQ=255)
               gCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTg                                                                             <  1:2991286/91‑1 (MQ=255)
                  ttGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGat                                                                          <  1:1246121/91‑1 (MQ=255)
                        aaTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGATATGCAGTGCATACATATTCTGGATATGaaa                                                                    <  1:1215518/91‑1 (MQ=255)
                         aTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGaaaa                                                                   <  1:3139269/91‑1 (MQ=255)
                         aTTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGaaaa                                                                   <  1:1454170/91‑1 (MQ=255)
                          tttATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAAt                                                                  <  1:4127713/91‑1 (MQ=255)
                            tattatTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATaa                                                                >  1:2861650/1‑91 (MQ=255)
                            tattatTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATaa                                                                <  1:1113417/91‑1 (MQ=255)
                            tattatTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATaa                                                                <  1:1336443/91‑1 (MQ=255)
                             attattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATaaa                                                               <  1:3625277/91‑1 (MQ=255)
                              ttattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAg                                                              <  1:3122798/91‑1 (MQ=255)
                              ttattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAg                                                              >  1:2034716/1‑91 (MQ=255)
                               tattTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGt                                                             <  1:338524/91‑1 (MQ=255)
                                attTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGtt                                                            >  1:1156886/1‑91 (MQ=255)
                                 tttCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTa                                                           >  1:1032502/1‑91 (MQ=255)
                                   tCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGCATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaa                                                         <  1:3592923/91‑1 (MQ=255)
                                    cGTTAGAGGAGTGAAGACATGTCTGAGGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTaaaa                                                        <  1:2681811/91‑1 (MQ=255)
                                      ttAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAg                                                      <  1:3892091/91‑1 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  1:3472432/91‑1 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  >  1:2873753/1‑91 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  >  1:2874181/1‑91 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  >  1:1660194/1‑91 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  1:318193/91‑1 (MQ=255)
                                          aggagTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAAtt                                                  <  1:1385256/91‑1 (MQ=255)
                                              gTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGa                                              <  1:2238288/91‑1 (MQ=255)
                                                 aaGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaa                                           >  1:70906/1‑91 (MQ=255)
                                                 aaGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaa                                           <  1:2686257/91‑1 (MQ=255)
                                                 aaGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaa                                           <  1:1823668/91‑1 (MQ=255)
                                                   gACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGaaaaaa                                         >  1:2923035/1‑91 (MQ=255)
                                                    aCATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAag                                        <  1:1752064/91‑1 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       >  1:2136541/1‑91 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       >  1:1492314/1‑91 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       >  1:436781/1‑91 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       >  1:255006/1‑91 (MQ=255)
                                                     cATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAaga                                       <  1:2062546/91‑1 (MQ=255)
                                                      aTGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAagag                                      <  1:3785837/91‑1 (MQ=255)
                                                      aTGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAagag                                      <  1:1762119/91‑1 (MQ=255)
                                                       tGTCTGATGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     <  1:2071046/91‑1 (MQ=255)
                                                       tGTCTGAGGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     <  1:3824851/91‑1 (MQ=255)
                                                       tGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGt                                     >  1:4050684/1‑91 (MQ=255)
                                                        gTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTa                                    >  1:1048452/1‑91 (MQ=255)
                                                         tCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAg                                   <  1:1975688/91‑1 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  <  1:892170/91‑1 (MQ=255)
                                                          cTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGa                                  >  1:2031773/1‑91 (MQ=255)
                                                              ctaagTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCt                              >  1:485215/3‑91 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           >  1:545981/1‑91 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           <  1:3154370/91‑1 (MQ=255)
                                                                 agTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATg                           >  1:2260319/1‑91 (MQ=255)
                                                                  gTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGa                          >  1:3939883/1‑91 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        >  1:3735540/1‑91 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        >  1:3576979/1‑91 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        >  1:3576887/1‑91 (MQ=255)
                                                                    tGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        <  1:2115375/91‑1 (MQ=255)
                                                                    tGGCGAAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAat                        <  1:2448987/91‑1 (MQ=255)
                                                                     ggCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAata                       >  1:3576587/1‑91 (MQ=255)
                                                                     ggCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAata                       >  1:3560876/1‑91 (MQ=255)
                                                                      gcgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAatat                      >  1:1940097/1‑91 (MQ=255)
                                                                       cgcAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATg                     >  1:1965406/1‑91 (MQ=255)
                                                                         cAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGAc                   <  1:1434490/91‑1 (MQ=255)
                                                                         cAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGAc                   <  1:1434526/91‑1 (MQ=255)
                                                                         cAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGAc                   <  1:276943/91‑1 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  1:1357005/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                <  1:3065318/91‑1 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  1:1144131/1‑91 (MQ=255)
                                                                            aGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGa                >  1:1144092/1‑91 (MQ=255)
                                                                             gCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGaa               <  1:592607/91‑1 (MQ=255)
                                                                             gCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGaa               <  1:2571952/91‑1 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  1:3368732/91‑1 (MQ=255)
                                                                              cGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAAc              <  1:881308/91‑1 (MQ=255)
                                                                                ttGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaa            <  1:2038670/91‑1 (MQ=255)
                                                                                tatgTATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACaa            >  1:3946015/4‑91 (MQ=255)
                                                                                  ggTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAAt          >  1:60485/1‑91 (MQ=255)
                                                                                   gTACGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATc         <  1:3136199/91‑1 (MQ=255)
                                                                                    tATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGGATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        <  1:901866/91‑1 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  1:2550084/1‑91 (MQ=255)
                                                                                    tATGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCt        >  1:2827340/1‑91 (MQ=255)
                                                                                     aTGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       <  1:2373835/91‑1 (MQ=255)
                                                                                     aTGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       >  1:1681830/1‑91 (MQ=255)
                                                                                     aTGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCtt       >  1:3567638/1‑91 (MQ=255)
                                                                                      tGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      >  1:33076/1‑91 (MQ=255)
                                                                                      tGCACTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCttt      >  1:2440527/1‑91 (MQ=255)
                                                                                       gCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     >  1:3512970/1‑91 (MQ=255)
                                                                                       gCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     <  1:316332/91‑1 (MQ=255)
                                                                                       gCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     <  1:2681561/91‑1 (MQ=255)
                                                                                       gCACTGCATCCATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTa     >  1:3054276/1‑91 (MQ=255)
                                                                                         aGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTaga   >  1:1208366/1‑91 (MQ=255)
                                                                                          gTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTagag  <  1:1878627/91‑1 (MQ=255)
                                                                                          gTGCATACATATTCCGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTagag  <  1:1878556/91‑1 (MQ=255)
                                                                                          |                                                                                          
GATTTAGGAATTACAGCATTGAAAAATTTATTATTTCGTTAGAGGAGTGAAGACATGTCTGAAGAAGTTGGCGCAAAGCGTTGGTATGCAGTGCATACATATTCTGGATATGAAAATAAAGTTAAAAAGAATTTAGAAAAAAGAGTAGAATCTATGAATATGACTGAACAAATCTTTAGAG  >  USA300TCH1516_ALE/581811‑581991

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