Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,233,504 T→A 12.1% V470D (GTT→GAT)  prkC → Serine/threonine‑protein kinase PrkC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,233,5040TA12.1% 362.3 / 34.2 116V470D (GTT→GAT) prkCSerine/threonine‑protein kinase PrkC
Reads supporting (aligned to +/- strand):  ref base T (70/32);  new base A (10/4);  total (80/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 125 total aligned reads displayed.

AGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGA  >  USA300TCH1516_ALE/1233414‑1233594
                                                                                          |                                                                                          
aGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGt                                                                                            >  1:1556546/1‑91 (MQ=255)
 gggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGtt                                                                                           <  1:1101363/91‑1 (MQ=255)
 gggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGtt                                                                                           >  1:4009531/1‑91 (MQ=255)
 gggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGtt                                                                                           <  1:3004219/91‑1 (MQ=255)
  ggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTa                                                                                          >  1:347682/1‑91 (MQ=255)
   gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTAc                                                                                         >  1:1079238/1‑91 (MQ=255)
   gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTAc                                                                                         >  1:2236471/1‑91 (MQ=255)
     ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGa                                                                                       >  1:2843336/1‑91 (MQ=255)
     ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGa                                                                                       >  1:1341321/1‑91 (MQ=255)
     ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGa                                                                                       >  1:3115746/1‑91 (MQ=255)
      cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGAt                                                                                      >  1:1409216/1‑91 (MQ=255)
      cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGAt                                                                                      >  1:2312424/1‑91 (MQ=255)
      cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGAt                                                                                      <  1:3973096/91‑1 (MQ=255)
      cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGAt                                                                                      >  1:2312253/1‑91 (MQ=255)
      cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGAt                                                                                      >  1:3414349/1‑91 (MQ=255)
        gAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGCTTg                                                                                    >  1:3738136/1‑91 (MQ=255)
        gAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTg                                                                                    >  1:511064/1‑91 (MQ=255)
        gAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTg                                                                                    <  1:2506834/91‑1 (MQ=255)
        gAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTg                                                                                    >  1:3493657/1‑91 (MQ=255)
         aaaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                                   >  1:3546098/1‑91 (MQ=255)
         aaaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                                   >  1:2108313/1‑91 (MQ=255)
           aaGGTTAAATTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGaaa                                                                                 <  1:3974277/91‑1 (MQ=255)
           aaGGTTAAACTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGaaa                                                                                 <  1:3974945/91‑1 (MQ=255)
           aaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGaaa                                                                                 <  1:3380558/91‑1 (MQ=255)
           aaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaa                                                                                 >  1:698514/1‑91 (MQ=255)
                tAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGta                                                                            >  1:1965897/1‑91 (MQ=255)
                tAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGta                                                                            >  1:3007873/1‑91 (MQ=255)
                tAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGta                                                                            >  1:1257120/1‑91 (MQ=255)
                 aaaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGACAAAGtat                                                                           >  1:2068530/1‑91 (MQ=255)
                     tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATataa                                                                       <  1:1717319/91‑1 (MQ=255)
                       ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATataata                                                                     >  1:2802732/1‑91 (MQ=255)
                       ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataata                                                                     <  1:1886149/91‑1 (MQ=255)
                        cAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATAtaataa                                                                    >  1:892957/1‑91 (MQ=255)
                         aaaTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATaataat                                                                   >  1:3057896/1‑91 (MQ=255)
                         aaaTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATaataat                                                                   >  1:2979050/1‑91 (MQ=255)
                          aaTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATc                                                                  >  1:264387/1‑91 (MQ=255)
                           ttggcATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCa                                                                 >  1:895877/5‑91 (MQ=255)
                           aTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCa                                                                 >  1:723569/1‑91 (MQ=255)
                            tGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCaa                                                                >  1:2575162/1‑91 (MQ=255)
                            tGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCaa                                                                >  1:1174626/1‑91 (MQ=255)
                            tGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCaa                                                                >  1:2501813/1‑91 (MQ=255)
                               cATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAgcg                                                             >  1:2731810/1‑91 (MQ=255)
                                aTTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAgcgc                                                            >  1:1692799/1‑91 (MQ=255)
                                  tGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAAGCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCa                                                          >  1:3504305/1‑91 (MQ=255)
                                     tttACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCaaaa                                                       <  1:4017491/91‑1 (MQ=255)
                                       tACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAgg                                                     >  1:3529451/1‑91 (MQ=255)
                                       tACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAgg                                                     >  1:650825/1‑91 (MQ=255)
                                        aCCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGa                                                    <  1:2546482/91‑1 (MQ=255)
                                        aCCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGa                                                    <  1:2546621/91‑1 (MQ=255)
                                         ccTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                                   >  1:4082909/1‑91 (MQ=255)
                                         ccTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                                   >  1:4083000/1‑91 (MQ=255)
                                           tAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATGTAATAATCAAGCGCCAAAAGGATAc                                                 >  1:3852653/1‑91 (MQ=255)
                                           tAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATGTAATAATCAAGCGCCAAAAGGATAc                                                 >  1:3853881/1‑91 (MQ=255)
                                           tAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATCTAATAATCAAGCGCCAAAAGGATAc                                                 >  1:3853140/1‑91 (MQ=255)
                                              ggaggaAGCCTTGCAGAATTTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                              >  1:3142622/1‑91 (MQ=255)
                                               gaggaAGCCTTGCAGAAATTAAAAGCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                             <  1:3442683/91‑1 (MQ=255)
                                                 ggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCa                                           <  1:2621540/91‑1 (MQ=255)
                                                 ggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCa                                           <  1:2621246/91‑1 (MQ=255)
                                                   aaGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCaaa                                         >  1:2263702/1‑91 (MQ=255)
                                                     gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATc                                       <  1:3351904/91‑1 (MQ=255)
                                                     gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATc                                       >  1:4112649/1‑91 (MQ=255)
                                                     gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATc                                       >  1:1613217/1‑91 (MQ=255)
                                                        ttGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaaa                                    >  1:1981270/1‑91 (MQ=255)
                                                         tGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAg                                   <  1:3979582/91‑1 (MQ=255)
                                                             gAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTa                               >  1:1141268/1‑91 (MQ=255)
                                                              aaaTTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTaa                              <  1:1147193/91‑1 (MQ=255)
                                                               aaTTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                             >  1:2471339/1‑91 (MQ=255)
                                                               aaTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAACTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                             >  1:3590836/1‑91 (MQ=255)
                                                                 ttAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCg                           <  1:2654585/91‑1 (MQ=255)
                                                                   aaaaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                         <  1:3367337/91‑1 (MQ=255)
                                                                   aaaaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                         <  1:3904697/91‑1 (MQ=255)
                                                                    aaaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                        >  1:229486/1‑91 (MQ=255)
                                                                    aaaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                        >  1:230015/1‑91 (MQ=255)
                                                                     aaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaaa                       >  1:1497575/1‑91 (MQ=255)
                                                                     aaTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaaa                       >  1:3957499/1‑91 (MQ=255)
                                                                      aTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAAt                      <  1:2577546/91‑1 (MQ=255)
                                                                      aTCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAAt                      >  1:4096916/1‑91 (MQ=255)
                                                                       tCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATa                     >  1:1323855/1‑91 (MQ=255)
                                                                       tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATa                     >  1:317925/1‑91 (MQ=255)
                                                                        cGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATAc                    <  1:2368875/91‑1 (MQ=255)
                                                                         gTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTATCCGCAAATACt                   <  1:3582202/91‑1 (MQ=255)
                                                                         gTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                   <  1:941937/91‑1 (MQ=255)
                                                                         gTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt                   <  1:942200/91‑1 (MQ=255)
                                                                          ttAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                  >  1:2912817/1‑91 (MQ=255)
                                                                          ttAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                  >  1:3261382/1‑91 (MQ=255)
                                                                           tAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGa                 >  1:936972/1‑91 (MQ=255)
                                                                           tAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGa                 >  1:2384637/1‑91 (MQ=255)
                                                                            aGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaa                <  1:2054589/91‑1 (MQ=255)
                                                                             ggTCTTAAAGATGTTACGATTGAAAGAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaaa               <  1:3876441/91‑1 (MQ=255)
                                                                              gTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAAt              <  1:2627650/91‑1 (MQ=255)
                                                                                cTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg            >  1:3477646/1‑91 (MQ=255)
                                                                                cTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg            <  1:3146323/91‑1 (MQ=255)
                                                                                cTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg            >  1:385982/1‑91 (MQ=255)
                                                                                     aGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAtt       >  1:3544544/1‑91 (MQ=255)
                                                                                     aGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAtt       >  1:542398/1‑91 (MQ=255)
                                                                                      gatgatACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTc      <  1:3628182/91‑1 (MQ=255)
                                                                                      gATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTc      <  1:1678139/91‑1 (MQ=255)
                                                                                        tGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCAt    <  1:239514/91‑1 (MQ=255)
                                                                                          ttACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGa  >  1:2190997/1‑91 (MQ=255)
                                                                                          ttACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGa  <  1:3469208/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGA  >  USA300TCH1516_ALE/1233414‑1233594

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