Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 582,522 C→G 7.2% intergenic (+109/‑72) nusG_2 → / → rplK Transcription termination/antitermination protein NusG/50S ribosomal protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE582,5220CG 7.2% 658.8 / 32.3 195intergenic (+109/‑72)nusG_2/rplKTranscription termination/antitermination protein NusG/50S ribosomal protein L11
Reads supporting (aligned to +/- strand):  ref base C (77/104);  new base G (9/5);  total (86/109)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 200 total aligned reads displayed.

AAACTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAGATA  >  USA300TCH1516_ALE/582432‑582612
                                                                                          |                                                                                          
aaaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTc                                                                                            <  1:530481/91‑1 (MQ=255)
aaaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTc                                                                                            >  1:2129925/1‑91 (MQ=255)
aaaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATCAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTc                                                                                            >  2:3251299/1‑91 (MQ=255)
 aaCTAACCGAAGGAGAAAAAAGAGGATTGATTTTTTAATTAGAAAAGTGTGAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           <  1:1022854/91‑1 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           <  1:3348275/91‑1 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           <  2:1438688/91‑1 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           >  1:2271097/1‑91 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           >  2:1603542/1‑91 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           <  1:1763261/91‑1 (MQ=255)
 aaCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCa                                                                                           >  1:240484/1‑91 (MQ=255)
  aCTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCAc                                                                                          <  1:334541/91‑1 (MQ=255)
   cTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTGACg                                                                                         <  2:2736846/91‑1 (MQ=255)
   cTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACg                                                                                         >  1:234808/1‑91 (MQ=255)
    tAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGa                                                                                        >  1:895047/1‑91 (MQ=255)
     aaCCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTGACGaa                                                                                       <  2:1458742/91‑1 (MQ=255)
     aaCCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGaa                                                                                       <  1:3774424/91‑1 (MQ=255)
      aCCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTGACGaaa                                                                                      <  2:2688031/91‑1 (MQ=255)
       ccAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAAt                                                                                     >  1:172088/1‑91 (MQ=255)
          aaaGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGtt                                                                                  >  1:1653915/1‑91 (MQ=255)
          aaaGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGtt                                                                                  >  1:3880416/1‑91 (MQ=255)
           aaGATAAAAAAGAGTATTGATTTTTTTTTTTTTTTAGTGTTAAAATTATGTGGTCGCGTTTTTAGAGCGCCCATTTCGTCACGAAATGTTa                                                                                 >  2:2183597/1‑91 (MQ=255)
           aaGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTa                                                                                 <  2:1556121/91‑1 (MQ=255)
             gATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGACATGTTAag                                                                               >  2:2888156/1‑91 (MQ=255)
             gATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAag                                                                               >  2:2282050/1‑91 (MQ=255)
             gATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAag                                                                               <  1:2612793/91‑1 (MQ=255)
             aaaaaaaaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAag                                                                               <  2:1642062/88‑1 (MQ=255)
                aaaaaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGt                                                                            <  2:174229/91‑1 (MQ=255)
                aaaaaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGt                                                                            <  1:2997955/91‑1 (MQ=255)
                  aaaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTgg                                                                          <  1:242747/91‑1 (MQ=255)
                  aaaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTgg                                                                          >  1:2545622/1‑91 (MQ=255)
                   aaaGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTggg                                                                         >  2:303588/1‑91 (MQ=255)
                    aaGAGTATTGATTTTTTAATTAGAAAAGTGTTATAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGa                                                                        <  1:3236014/91‑1 (MQ=255)
                      gagTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGCGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAgg                                                                      >  2:3120229/1‑91 (MQ=255)
                       agTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAggg                                                                     <  1:1349767/91‑1 (MQ=255)
                       agTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAggg                                                                     <  1:823117/91‑1 (MQ=255)
                       agTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAggg                                                                     <  1:1348354/91‑1 (MQ=255)
                        gTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGc                                                                    <  2:3656268/91‑1 (MQ=255)
                         tATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCa                                                                   <  1:3524052/91‑1 (MQ=255)
                          aTTGATTTTTTAATTATAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCaa                                                                  <  2:3034556/91‑1 (MQ=255)
                               ttttttAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTg                                                             <  2:2065700/91‑1 (MQ=255)
                                tttttAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGa                                                            <  2:792329/91‑1 (MQ=255)
                                 ttttAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAg                                                           <  1:656826/91‑1 (MQ=255)
                                 ttttAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAg                                                           <  1:1789776/91‑1 (MQ=255)
                                   ttAAATAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGACCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGcc                                                         <  1:2191140/91‑1 (MQ=255)
                                      aTTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCtg                                                      >  2:781727/1‑91 (MQ=255)
                                      aTTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCtg                                                      >  1:1432829/1‑91 (MQ=255)
                                       ttAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGCACTGTTAAGCGTGGGAGGGCAAAACTGAGCCCcgt                                                     >  2:475181/1‑91 (MQ=255)
                                       ttAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCtgt                                                     >  2:3188835/1‑91 (MQ=255)
                                           aaaaGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCGCGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGAcc                                                 <  2:500106/91‑1 (MQ=255)
                                           aaaaGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGAcc                                                 <  1:1112368/91‑1 (MQ=255)
                                           aaaaGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGAcc                                                 <  1:415896/91‑1 (MQ=255)
                                            aaaGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACca                                                <  1:1383391/91‑1 (MQ=255)
                                            aaaGTGTTAAAAGTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACca                                                >  1:3029839/1‑91 (MQ=255)
                                               gtgtTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGCCCACAt                                             >  2:2925778/1‑91 (MQ=255)
                                                tgtTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATc                                            >  1:125286/1‑91 (MQ=255)
                                                 gtTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGTAAAACTGAGCCCTGTGACCACATCa                                           >  1:3673139/1‑91 (MQ=255)
                                                 gtTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCa                                           <  1:1472408/91‑1 (MQ=255)
                                                  ttATAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCAc                                          <  2:2410101/91‑1 (MQ=255)
                                                      aaTTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGata                                      <  2:238090/91‑1 (MQ=255)
                                                       aTTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGatat                                     <  2:1057794/91‑1 (MQ=255)
                                                       aTTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGatat                                     <  1:1689138/91‑1 (MQ=255)
                                                        ttATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATc                                    <  2:1276938/91‑1 (MQ=255)
                                                        ttATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTGAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATc                                    <  1:3134424/91‑1 (MQ=255)
                                                           tgtgGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAag                                 <  2:2471588/91‑1 (MQ=255)
                                                             tgGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAagga                               <  2:2529922/91‑1 (MQ=255)
                                                             tgGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAagga                               >  2:3803122/1‑91 (MQ=255)
                                                             tgGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAagga                               >  1:186637/1‑91 (MQ=255)
                                                             tgGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAagga                               >  1:2283106/1‑91 (MQ=255)
                                                              ggTCGCGCTTACCGAGCGCCCATTTCGGCACGTAACGTCCAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACTTCTCGTTATCAaggag                              >  2:3583669/1‑91 (MQ=255)
                                                                  gcgcTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGc                          >  1:2528487/1‑91 (MQ=255)
                                                                  gcgcTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGc                          <  2:1092403/91‑1 (MQ=255)
                                                                   cgcTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGca                         >  1:2092056/1‑91 (MQ=255)
                                                                    gcTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGcac                        <  1:1576142/91‑1 (MQ=255)
                                                                       tttAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATc                     >  2:3847828/1‑91 (MQ=255)
                                                                         tAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGt                   >  2:2952866/1‑91 (MQ=255)
                                                                         tAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGt                   >  1:1307588/1‑91 (MQ=255)
                                                                           gagCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGCGCACATCGTgg                 >  1:3497904/1‑91 (MQ=255)
                                                                           gagCGCCCATTTCGTCACGAAAGGTTAAGGGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATAAGGGAGGTGCACATCGTgg                 <  2:642883/91‑1 (MQ=255)
                                                                            agCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGc                <  1:401302/91‑1 (MQ=255)
                                                                             gcgcCCATTTCGTCATGAAATGTTACTAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCt               <  2:3420375/91‑1 (MQ=255)
                                                                             gcgcCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCt               <  2:224987/91‑1 (MQ=255)
                                                                             gcgcCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCt               <  1:2378752/91‑1 (MQ=255)
                                                                             gcgcCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACGGAGCCCTGTGACCACATCACGATAGCCAGGAGGTGCACATCGTGGCt               <  2:2041295/91‑1 (MQ=255)
                                                                               gcCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTaa             >  1:3278431/1‑91 (MQ=255)
                                                                               gcCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTaa             <  2:1931376/91‑1 (MQ=255)
                                                                                cccATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTaaa            <  2:1369414/91‑1 (MQ=255)
                                                                                   aTTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTaaaaaa         >  1:1044813/1‑91 (MQ=255)
                                                                                   aTTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTaaaaaa         <  1:908106/91‑1 (MQ=255)
                                                                                    tttCGTGACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAg        >  1:574084/1‑91 (MQ=255)
                                                                                    tttCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAg        <  1:2156529/91‑1 (MQ=255)
                                                                                     ttCGTGACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGt       >  1:774566/1‑91 (MQ=255)
                                                                                     ttCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGCTATCAAGGAGGTGCACATCGTGGCTAAAAAAGt       >  1:2705041/1‑91 (MQ=255)
                                                                                     ttCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGt       >  2:103480/1‑91 (MQ=255)
                                                                                     ttCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGt       <  2:1537424/91‑1 (MQ=255)
                                                                                      tCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTa      <  2:455738/91‑1 (MQ=255)
                                                                                       cGTGACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTCAAAAAGTAg     >  1:804823/1‑91 (MQ=255)
                                                                                       cGTGACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAg     >  2:1251771/1‑91 (MQ=255)
                                                                                        gTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAGa    <  1:1311077/91‑1 (MQ=255)
                                                                                         tCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAGAt   >  1:1364458/1‑91 (MQ=255)
                                                                                          cACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAGATa  >  1:1800171/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAACTAACCAAAAGATAAAAAAGAGTATTGATTTTTTAATTAGAAAAGTGTTAAAATTATGTGGTCGCGCTTTTAGAGCGCCCATTTCGTCACGAAATGTTAAGAGTGGGAGGGCAAAACTGAGCCCTGTGACCACATCACGATATCAAGGAGGTGCACATCGTGGCTAAAAAAGTAGATA  >  USA300TCH1516_ALE/582432‑582612

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