Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 594,888 T→A 5.7% intergenic (+42/‑81) rpsG → / → fusA 30S ribosomal protein S7/Elongation factor G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE594,8880TA 5.7% 689.8 / 21.6 204intergenic (+42/‑81)rpsG/fusA30S ribosomal protein S7/Elongation factor G
Reads supporting (aligned to +/- strand):  ref base T (99/93);  new base A (5/7);  total (104/100)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CACAAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTAGAG  >  USA300TCH1516_ALE/594799‑594978
                                                                                         |                                                                                          
cacaAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAttt                                                                                           <  2:3228634/91‑1 (MQ=255)
cacaAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAttt                                                                                           >  1:987587/1‑91 (MQ=255)
cacaAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAttt                                                                                           >  1:2742208/1‑91 (MQ=255)
 acaAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAtttt                                                                                          >  2:2689552/1‑91 (MQ=255)
 acaAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAtttt                                                                                          >  2:1911753/1‑91 (MQ=255)
  caAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTc                                                                                         >  2:747635/1‑91 (MQ=255)
  caAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTc                                                                                         <  2:2070107/91‑1 (MQ=255)
   aaaaTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCa                                                                                        <  2:1702764/91‑1 (MQ=255)
   aaaaTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAAAGAATTTTCa                                                                                        >  1:3709105/1‑91 (MQ=255)
    aaaTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCAt                                                                                       <  2:1961139/91‑1 (MQ=255)
     aaTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTATCAtt                                                                                      >  2:3146736/1‑91 (MQ=255)
      aTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTa                                                                                     <  2:199530/91‑1 (MQ=255)
      aTGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTa                                                                                     <  2:3349734/91‑1 (MQ=255)
       tGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAAATTCATTaa                                                                                    <  1:3762667/91‑1 (MQ=255)
       tGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATAAATTTTCATTaa                                                                                    <  1:2993355/91‑1 (MQ=255)
        ggCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAg                                                                                   >  1:2947206/1‑91 (MQ=255)
        ggCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAg                                                                                   >  2:6706/1‑91 (MQ=255)
        ggCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAg                                                                                   >  2:3022292/1‑91 (MQ=255)
         gCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGc                                                                                  >  1:511493/1‑91 (MQ=255)
         gCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAAAATTCATTaata                                                                                  <  1:1638977/91‑3 (MQ=255)
           tGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGt                                                                                <  2:2456742/91‑1 (MQ=255)
            gAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGtt                                                                               >  2:3028089/1‑91 (MQ=255)
                caaacaaaGCATTTGCTCACTACCGTTGGTAAGGTAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATg                                                                           <  1:3450676/91‑1 (MQ=255)
                   acaaaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCtt                                                                        <  1:228346/91‑1 (MQ=255)
                   acaaaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCtt                                                                        <  1:2005128/91‑1 (MQ=255)
                    caaaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTa                                                                       >  1:776265/1‑91 (MQ=255)
                     aaaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAg                                                                      <  1:453603/91‑1 (MQ=255)
                     aaaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAg                                                                      <  1:1329069/91‑1 (MQ=255)
                      aaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAgg                                                                     <  1:1305948/91‑1 (MQ=255)
                      aaGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAgg                                                                     >  2:3529913/1‑91 (MQ=255)
                        gCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGc                                                                   >  1:1175523/1‑91 (MQ=255)
                        gCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGc                                                                   >  1:108993/1‑91 (MQ=255)
                         cATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCa                                                                  >  2:1583223/1‑91 (MQ=255)
                         cATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCa                                                                  <  2:3277545/91‑1 (MQ=255)
                          aTTTGCTCACTACCGCTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAACGAATTTTCATTAAGCGTTCATGCTTAGGGCAt                                                                 >  1:3302963/1‑91 (MQ=255)
                           tttGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATc                                                                >  1:1159175/1‑91 (MQ=255)
                              gCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGcc                                                             >  2:1316925/1‑91 (MQ=255)
                              gCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGcc                                                             <  2:1768943/91‑1 (MQ=255)
                              gCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGcc                                                             >  1:1456132/1‑91 (MQ=255)
                                  aCTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCatat                                                         >  2:512132/1‑91 (MQ=255)
                                  aCTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCatat                                                         <  2:3486399/91‑1 (MQ=255)
                                     aCCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCAtatcta                                                      >  2:380780/1‑91 (MQ=255)
                                       cGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCTTCGCCAtatctatc                                                    <  2:170400/91‑1 (MQ=255)
                                       cGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCAtatctatc                                                    >  2:54977/1‑91 (MQ=255)
                                       cGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCAtatctatc                                                    <  2:2610827/91‑1 (MQ=255)
                                       cGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCAtatctatc                                                    <  1:2359963/91‑1 (MQ=255)
                                         ttGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGt                                                  <  1:2704107/91‑1 (MQ=255)
                                           ggTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGtat                                                >  2:3745325/1‑91 (MQ=255)
                                           ggTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGtat                                                >  1:2137182/1‑91 (MQ=255)
                                             tAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGtattt                                              <  2:2241233/91‑1 (MQ=255)
                                              aaGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGtattta                                             >  1:3423573/1‑91 (MQ=255)
                                                gATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGtatttatt                                           >  1:3024897/1‑91 (MQ=255)
                                                      aGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTaa                                     <  2:2486384/91‑1 (MQ=255)
                                                       gCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAat                                    <  2:1808114/91‑1 (MQ=255)
                                                          tttACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAAtata                                 >  1:1122670/1‑91 (MQ=255)
                                                          tttACCCTGAGTATGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTCGGGCATCGCCATATCTATCGTATTTATTCAGTAAtata                                 >  1:799150/1‑91 (MQ=255)
                                                           ttACCCTGAGTGTGCTCTATATTAATGAACTTTCTCGCAGCGGTCATTCTTAGGGCATCGCTATATCTATCGTAATTATTCAGTAatgtta                                >  1:1770796/1‑87 (MQ=255)
                                                             aCCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAAc                              <  2:2699463/91‑1 (MQ=255)
                                                             aCCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAAc                              <  1:1100954/91‑1 (MQ=255)
                                                              cccTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACt                             >  1:2309577/1‑91 (MQ=255)
                                                              cccTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACt                             <  2:2604002/91‑1 (MQ=255)
                                                                cTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTgg                           <  1:3114985/91‑1 (MQ=255)
                                                                   aGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAg                        <  1:475687/91‑1 (MQ=255)
                                                                   aGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAg                        <  1:474560/91‑1 (MQ=255)
                                                                   aGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAg                        >  2:203210/1‑91 (MQ=255)
                                                                   aGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAg                        <  2:2917636/91‑1 (MQ=255)
                                                                   aGTATGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAg                        <  1:3490803/91‑1 (MQ=255)
                                                                    gtgtgtTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTgg                           >  2:3839931/1‑87 (MQ=255)
                                                                        gtTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaa                   >  1:2625923/1‑91 (MQ=255)
                                                                        gtTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaa                   >  1:3554033/1‑91 (MQ=255)
                                                                          tCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaaaa                 <  1:2427724/91‑1 (MQ=255)
                                                                          tCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaaaa                 <  2:1074226/91‑1 (MQ=255)
                                                                           cTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaaaaa                <  2:3342250/91‑1 (MQ=255)
                                                                            tataTTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaaaaaa               >  1:1048191/1‑91 (MQ=255)
                                                                            tataTTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGaaaaaa               <  1:3022292/91‑1 (MQ=255)
                                                                             atatTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAAt              <  1:1671665/91‑1 (MQ=255)
                                                                             atatTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAAt              <  1:1563130/91‑1 (MQ=255)
                                                                               atTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATAc            <  1:2418269/91‑1 (MQ=255)
                                                                                ttAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACa           >  2:2190691/1‑91 (MQ=255)
                                                                                ttAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACa           >  1:3805400/1‑91 (MQ=255)
                                                                                 tAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACAt          >  2:116628/1‑91 (MQ=255)
                                                                                 tAATGAAAATTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTCATTCCGTAATATAAACTGGAAGGAGAAAAAATACAt          >  2:3053509/1‑91 (MQ=255)
                                                                                 tAATGAAAATTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTCATTCAGTAATATAAACTGGAAGGAGAAAAAATACAt          >  2:3053748/1‑91 (MQ=255)
                                                                                  aacgAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATg         >  2:2320584/4‑91 (MQ=255)
                                                                                  aaTGAAAATTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATg         >  1:558196/1‑91 (MQ=255)
                                                                                    tGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGc       >  2:2604330/1‑91 (MQ=255)
                                                                                    tGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGc       >  2:163104/1‑91 (MQ=255)
                                                                                     gAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCt      >  1:2292969/1‑91 (MQ=255)
                                                                                     gAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCt      >  2:1132335/1‑91 (MQ=255)
                                                                                     gAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCt      >  2:1466776/1‑91 (MQ=255)
                                                                                      aaTTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTa     <  1:3369964/91‑1 (MQ=255)
                                                                                      aaTTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTa     >  1:2367857/1‑91 (MQ=255)
                                                                                      aaTTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTa     <  2:447091/91‑1 (MQ=255)
                                                                                       aTTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTag    >  1:3846121/1‑91 (MQ=255)
                                                                                        ttttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTaga   >  2:1169935/1‑91 (MQ=255)
                                                                                         tttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTagag  <  2:2586501/91‑1 (MQ=255)
                                                                                         tttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTagag  <  1:1224440/91‑1 (MQ=255)
                                                                                         tttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTagag  <  1:1221999/91‑1 (MQ=255)
                                                                                         tttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTagag  >  1:3853880/1‑91 (MQ=255)
                                                                                         tttCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTagag  >  2:1040947/1‑91 (MQ=255)
                                                                                         |                                                                                          
CACAAAATGGCTGAAGCAAACAAAGCATTTGCTCACTACCGTTGGTAAGATAAAAGCTTTTACCCTGAGTGTGTTCTATATTAATGAATTTTCATTAAGCGTTCATGCTTAGGGCATCGCCATATCTATCGTATTTATTCAGTAATATAAACTGGAAGGAGAAAAAATACATGGCTAGAG  >  USA300TCH1516_ALE/594799‑594978

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: