Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,332 A→G 9.3% intergenic (‑37/+122) gcvT ← / ← aroK Aminomethyltransferase/Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,3320AG 9.3% 623.5 / 41.1 215intergenic (‑37/+122)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base A (79/115);  new base G (11/9);  total (90/125)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 216 total aligned reads displayed.

ACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACC  >  USA300TCH1516_ALE/1653242‑1653421
                                                                                          |                                                                                         
aCCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCa                                                                                           >  1:1076796/1‑91 (MQ=255)
  ctctATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAAt                                                                                         <  2:178107/91‑1 (MQ=255)
     tATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGc                                                                                      <  1:3439879/91‑1 (MQ=255)
     tATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGc                                                                                      <  1:3392909/91‑1 (MQ=255)
         aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                                  <  2:744659/91‑3 (MQ=255)
         aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                                  <  1:1797288/91‑3 (MQ=255)
           cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                <  1:1250349/91‑1 (MQ=255)
           cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                >  1:1398890/1‑91 (MQ=255)
           cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                <  1:2855907/91‑1 (MQ=255)
            aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTg                                                                               <  1:728596/91‑1 (MQ=255)
            aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTg                                                                               <  1:2691721/91‑1 (MQ=255)
               aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                            <  1:3665012/91‑1 (MQ=255)
               aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                            <  1:290332/91‑1 (MQ=255)
                ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGtttt                                                                           <  2:416652/91‑2 (MQ=37)
                ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTta                                                                           >  1:2873006/1‑91 (MQ=255)
                 tttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttttt                                                                          <  2:45115/91‑3 (MQ=37)
                   tGATATAAAGGTGTTTGTTATAAATCACTTGGCATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGTTtttat                                                                        <  2:2600321/91‑2 (MQ=21)
                     atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                      <  2:158894/91‑1 (MQ=255)
                     atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                      <  1:2814374/91‑1 (MQ=255)
                      tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                     <  2:2277439/91‑1 (MQ=255)
                      tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                     <  2:3039809/91‑1 (MQ=255)
                      tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                     >  2:3235170/1‑91 (MQ=255)
                       ataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCa                                                                    <  1:1207389/91‑1 (MQ=255)
                        taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                   >  2:184749/1‑91 (MQ=255)
                        taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                   <  1:2564777/91‑1 (MQ=255)
                         aaaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAcc                                                                  <  2:1352545/91‑1 (MQ=255)
                            ggTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAt                                                               <  1:3537128/91‑1 (MQ=255)
                               gtttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTcg                                                            <  1:2138002/91‑2 (MQ=255)
                                  tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAGTCTGTTTATAAATCACCAATTCAGGa                                                         >  2:3354732/1‑91 (MQ=255)
                                  tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                         <  1:1842948/91‑1 (MQ=255)
                                   gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                        <  2:320320/91‑1 (MQ=255)
                                   gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                        >  1:2781441/1‑91 (MQ=255)
                                     tttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcg                                                      <  1:3061295/91‑1 (MQ=255)
                                        aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                   >  2:3018843/1‑91 (MQ=255)
                                         aaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGtt                                                  <  1:209502/91‑1 (MQ=255)
                                           tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                <  1:1467438/91‑1 (MQ=255)
                                           tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                <  1:2010971/91‑1 (MQ=255)
                                                 gACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatt                                          >  1:434129/1‑91 (MQ=255)
                                                  aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                         <  2:1726527/91‑1 (MQ=255)
                                                  aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                         <  2:710712/91‑1 (MQ=255)
                                                   cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                        <  2:2936010/91‑1 (MQ=255)
                                                    aTTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATc                                       <  2:2156268/91‑1 (MQ=255)
                                                     ttGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATcc                                      <  1:1501680/91‑1 (MQ=255)
                                                      tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                     <  1:1752070/91‑1 (MQ=255)
                                                      tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                     <  1:526449/91‑1 (MQ=255)
                                                        aaaTACCGCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCttt                                   <  2:3451226/91‑1 (MQ=255)
                                                         aaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtttt                                  >  1:3140892/1‑91 (MQ=255)
                                                          aTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTg                                 >  1:2549146/1‑91 (MQ=255)
                                                            ccccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTAcc                               >  2:3705507/2‑91 (MQ=255)
                                                            aCCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGcc                               >  1:379489/1‑91 (MQ=255)
                                                             cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                              >  2:1881392/1‑91 (MQ=255)
                                                             cccccTTAAAATTATAAAAACAGGATAGCAATAGCTACAATCTTTATAAATCACCAATTGATGAACGCGTTACAGTATTATCCTTTTGCCt                              >  2:2414271/1‑91 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  1:1298476/91‑1 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  2:2781441/91‑1 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  1:808405/91‑1 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  1:134164/91‑1 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  1:1298946/91‑1 (MQ=255)
                                                              ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                             <  1:2184463/91‑1 (MQ=255)
                                                                 cTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTgaga                          <  2:2367052/91‑1 (MQ=255)
                                                                   tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGt                        <  1:2364356/91‑1 (MQ=255)
                                                                   tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGt                        >  2:1445004/1‑91 (MQ=255)
                                                                     aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                      >  1:2780774/1‑91 (MQ=255)
                                                                     aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                      <  2:386545/91‑1 (MQ=255)
                                                                     aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                      <  1:1733762/91‑1 (MQ=255)
                                                                      aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                     <  1:843805/91‑1 (MQ=255)
                                                                        ttATAAAAACAGGATAGCCATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                   <  1:2521752/91‑1 (MQ=255)
                                                                        ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                   >  1:630409/1‑91 (MQ=255)
                                                                         tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                  <  1:3170204/91‑1 (MQ=255)
                                                                         tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                  >  1:2340442/1‑91 (MQ=255)
                                                                          ataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaa                 >  1:1275693/1‑91 (MQ=255)
                                                                           taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                >  1:1009392/1‑91 (MQ=255)
                                                                           taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                >  1:897482/1‑91 (MQ=255)
                                                                           taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                >  1:1667455/1‑91 (MQ=255)
                                                                            aaaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaata               <  1:2998387/91‑1 (MQ=255)
                                                                             taaaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa              >  2:2516290/2‑91 (MQ=37)
                                                                             aaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa              <  2:1133919/91‑1 (MQ=255)
                                                                             aaaaCAGGATAGCAAGCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa              <  2:3500812/91‑1 (MQ=255)
                                                                               aaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGCACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt            >  2:2164945/1‑91 (MQ=37)
                                                                               aaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt            >  1:3579694/1‑91 (MQ=38)
                                                                               aaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt            <  1:1359989/91‑1 (MQ=255)
                                                                                aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt           <  2:216851/91‑1 (MQ=38)
                                                                                aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt           >  1:777805/1‑91 (MQ=38)
                                                                                aCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt           <  2:795115/91‑1 (MQ=255)
                                                                                 cAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa          >  1:1912169/1‑91 (MQ=255)
                                                                                   ggATAGCAATCGCTACCCTGTTTATAAATCACCAGTGCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACt        <  2:149733/91‑1 (MQ=255)
                                                                                    ggtAGCGATTGCTATCCTGTTTATAAATCCCCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATGAATTACtt       >  2:3378554/3‑91 (MQ=255)
                                                                                    gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt       >  2:2093267/1‑91 (MQ=255)
                                                                                     aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg      <  1:2405989/91‑1 (MQ=255)
                                                                                     aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg      <  2:900672/91‑1 (MQ=255)
                                                                                     aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg      <  2:1250370/91‑1 (MQ=255)
                                                                                      tAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc     >  2:1206470/1‑91 (MQ=255)
                                                                                      tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc     >  2:2777913/1‑91 (MQ=255)
                                                                                      tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc     <  1:991289/91‑1 (MQ=255)
                                                                                      tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc     <  2:1333074/91‑1 (MQ=255)
                                                                                       aGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa    >  1:3394319/1‑91 (MQ=255)
                                                                                       aGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa    >  2:1354518/1‑91 (MQ=255)
                                                                                        gCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAc   <  2:691783/91‑1 (MQ=255)
                                                                                         cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc  >  1:3819373/1‑91 (MQ=255)
                                                                                         cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc  <  2:867168/91‑1 (MQ=255)
                                                                                         cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc  <  1:3290513/91‑1 (MQ=255)
                                                                                          |                                                                                         
ACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACC  >  USA300TCH1516_ALE/1653242‑1653421

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