Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,046,831 A→G 5.8% intergenic (+39/+50) putP → / ← USA300HOU_RS10335 Sodium/proline symporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,046,8310AG 5.8% 407.1 / 10.4 121intergenic (+39/+50)putP/USA300HOU_RS10335Sodium/proline symporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (63/51);  new base G (3/4);  total (66/55)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 126 total aligned reads displayed.

ATTTGTTGAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCATCT  >  USA300TCH1516_ALE/2046741‑2046920
                                                                                          |                                                                                         
aTTTGTTGAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTaa                                                                                           <  1:3255258/91‑1 (MQ=255)
aTTTGTTGAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTaa                                                                                           <  1:3253776/91‑1 (MQ=255)
 tttGTTGAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAAt                                                                                          <  1:34008/91‑1 (MQ=255)
   tgttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCa                                                                                        >  1:3608927/1‑91 (MQ=255)
   tgttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCa                                                                                        >  1:4069078/1‑91 (MQ=255)
   tgttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCa                                                                                        >  1:3610610/1‑91 (MQ=255)
    gttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCaa                                                                                       >  1:264345/1‑91 (MQ=255)
    gttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCaa                                                                                       <  1:499098/91‑1 (MQ=255)
     ttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCaaa                                                                                      >  1:1758536/1‑91 (MQ=255)
     ttgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCaaa                                                                                      >  1:3650652/1‑91 (MQ=255)
      tgAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGACTTAACAAATTAAAAAGGTACTAATGTTAATCaaaa                                                                                     >  1:1566021/1‑91 (MQ=255)
       gAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAAt                                                                                    >  1:3367320/1‑91 (MQ=255)
        aaaCTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAAtt                                                                                   >  1:375146/1‑91 (MQ=255)
         aaCTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTa                                                                                  <  1:14581/91‑1 (MQ=255)
         aaCTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTa                                                                                  >  1:883264/1‑91 (MQ=255)
         aaCTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTa                                                                                  <  1:175042/91‑1 (MQ=255)
         aaCTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTa                                                                                  <  1:1206185/91‑1 (MQ=255)
           cTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATg                                                                                >  1:293213/1‑91 (MQ=255)
           cTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATg                                                                                >  1:293353/1‑91 (MQ=255)
           cTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATg                                                                                <  1:3930771/91‑1 (MQ=255)
             gACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGAc                                                                              >  1:1698533/1‑91 (MQ=255)
              aCTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACt                                                                             >  1:1768084/1‑91 (MQ=255)
              aCTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACt                                                                             <  1:2415385/91‑1 (MQ=255)
              aCTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACt                                                                             <  1:2415406/91‑1 (MQ=255)
               cTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTa                                                                            >  1:3067761/1‑91 (MQ=255)
                ttAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTaa                                                                           >  1:3992918/1‑91 (MQ=255)
                ttAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTaa                                                                           <  1:2439541/91‑1 (MQ=255)
                 tAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAAc                                                                          <  1:497199/91‑1 (MQ=255)
                  aaaCAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACa                                                                         >  1:2068876/1‑91 (MQ=255)
                  aaaCAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACa                                                                         <  1:2916674/91‑1 (MQ=255)
                   aaCAAAGTTCGTGACATCGTTAGAGAAAAATAATTCACAAGTCTTAACAAATTAAAAAGGTACTATTGTTAATCAAAATTATGACAAACAt                                                                        >  1:3860786/1‑91 (MQ=255)
                    aCAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACAtt                                                                       <  1:3161572/91‑1 (MQ=255)
                     cAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAGTCATAATTTTGATTAACAtta                                                                      <  1:180184/91‑2 (MQ=25)
                      aaaGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTgg                                                                     <  1:944748/91‑1 (MQ=255)
                      aaaGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTgg                                                                     <  1:1394730/91‑1 (MQ=255)
                       aaGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGt                                                                    <  1:2112330/91‑1 (MQ=255)
                       aaGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGt                                                                    <  1:2110735/91‑1 (MQ=255)
                         gTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTAc                                                                  >  1:3721783/1‑91 (MQ=255)
                         gTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTAc                                                                  <  1:2832743/91‑1 (MQ=255)
                         gTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTAc                                                                  >  1:1274317/1‑91 (MQ=255)
                           tCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACCAATGTTAGTCATAATTTTGATTAACATTAGTACCt                                                                <  1:896496/91‑1 (MQ=17)
                            cGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCtt                                                               >  1:468327/1‑91 (MQ=255)
                             gTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCttt                                                              <  1:2798554/91‑1 (MQ=255)
                             gTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACCAATGTTAGTCATAATTTTGATTAACATTAGTACCttt                                                              <  1:1619260/91‑1 (MQ=17)
                               gACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCttttt                                                            >  1:1202090/1‑91 (MQ=255)
                                  aTCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTttatt                                                         <  1:3678456/91‑1 (MQ=255)
                                   tCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTttatta                                                        <  1:3624257/91‑1 (MQ=255)
                                    cGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTtattat                                                       <  1:1050682/91‑1 (MQ=255)
                                     gTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATc                                                      <  1:3358248/91‑1 (MQ=255)
                                     gTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATc                                                      >  1:3062498/1‑91 (MQ=255)
                                      ttAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCt                                                     <  1:3242632/91‑1 (MQ=255)
                                        agagAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCttt                                                   >  1:3956504/1‑91 (MQ=255)
                                           gaAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAAt                                                >  1:1376852/1‑91 (MQ=255)
                                           gaAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAAt                                                >  1:3528025/1‑91 (MQ=255)
                                             aaaaTAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAAttt                                              <  1:877043/91‑1 (MQ=255)
                                                ataatTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACa                                           >  1:2982824/1‑91 (MQ=255)
                                                ataatTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACa                                           >  1:2982584/1‑91 (MQ=255)
                                                 taatTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaa                                          <  1:2720249/91‑1 (MQ=255)
                                                  aatTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaat                                         >  1:435233/1‑91 (MQ=255)
                                                   atTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaatt                                        >  1:3393846/1‑91 (MQ=255)
                                                     tCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaattaa                                      >  1:3144711/1‑91 (MQ=255)
                                                     tCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaattaa                                      >  1:4032880/1‑91 (MQ=255)
                                                      cATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaattaat                                     <  1:3749862/91‑1 (MQ=255)
                                                       aTAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaattaatt                                    <  1:3934502/91‑1 (MQ=255)
                                                       aTAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACaattaatt                                    >  1:2048753/1‑91 (MQ=255)
                                                        tAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAattaatta                                   <  1:475197/91‑1 (MQ=255)
                                                           gTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTg                                >  1:1331397/1‑91 (MQ=255)
                                                            tCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGt                               <  1:907273/91‑1 (MQ=255)
                                                            tCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGt                               <  1:907088/91‑1 (MQ=255)
                                                             cTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTa                              >  1:440071/1‑91 (MQ=255)
                                                               tAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTaaa                            <  1:816114/91‑1 (MQ=255)
                                                                aaCAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAat                           <  1:601140/91‑1 (MQ=255)
                                                                aaCAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAat                           >  1:767525/1‑91 (MQ=255)
                                                                aaCAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAat                           >  1:767146/1‑91 (MQ=255)
                                                                 aCAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAata                          <  1:384299/91‑1 (MQ=255)
                                                                    taatAAAAAGGTACCAATGTTAGTCATAATTTTGATTAACATTAGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGa                       >  1:2520247/4‑91 (MQ=14)
                                                                    aaTTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGa                       >  1:3108635/1‑91 (MQ=255)
                                                                       tAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACt                    >  1:625849/1‑91 (MQ=255)
                                                                         aaaaGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTg                  <  1:868413/91‑1 (MQ=255)
                                                                         aaaaGGTACCAATGTTAGTCATAATTTTGATTAACATTAGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTg                  >  1:1335606/1‑91 (MQ=17)
                                                                          aaaGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGc                 >  1:1591344/1‑91 (MQ=255)
                                                                           aaGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCg                <  1:1423183/91‑1 (MQ=255)
                                                                              gTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGtt             >  1:2308059/1‑91 (MQ=255)
                                                                                 cgaaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCtt          <  1:2246897/89‑1 (MQ=255)
                                                                                 cTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCtt          >  1:2280401/1‑91 (MQ=255)
                                                                                 cTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCtt          <  1:4008247/91‑1 (MQ=255)
                                                                                 cTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCtt          <  1:1776423/91‑1 (MQ=255)
                                                                                  tAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCttt         >  1:1880476/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATTAACTTGCGGTTCTTTg        >  1:2698802/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        <  1:1763107/91‑1 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        >  1:2992617/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        >  1:1989974/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        >  1:1990565/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        >  1:2698303/1‑91 (MQ=255)
                                                                                   aaTGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTg        >  1:1075737/1‑91 (MQ=255)
                                                                                     tGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGtt      >  1:760481/1‑90 (MQ=255)
                                                                                      gTTAGTCATAATTTTGATTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCa     >  1:3689525/1‑91 (MQ=37)
                                                                                         aaTCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCATCt  <  1:2750932/91‑1 (MQ=255)
                                                                                         aaTCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCATCt  <  1:2916122/91‑1 (MQ=255)
                                                                                         aaTCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCATCt  <  1:2750786/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATTTGTTGAAACTGACTTAAACAAAGTTCGTGACATCGTTAGAGAAAAATAATTCATAAGTCTTAACAAATTAAAAAGGTACTAATGTTAATCAAAATTATGACTAACATTGGTACCTTTTTATTATCTTTAATTTACAATTAATTACTGTAAATATGAACTTGCGGTTCTTTGTCATCT  >  USA300TCH1516_ALE/2046741‑2046920

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: