Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,335 C→T 8.9% intergenic (‑40/+119) gcvT ← / ← aroK Aminomethyltransferase/Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,3350CT 8.9% 684.8 / 35.1 221intergenic (‑40/+119)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base C (84/116);  new base T (11/9);  total (96/125)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCG  >  USA300TCH1516_ALE/1653247‑1653425
                                                                                        |                                                                                          
tATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGc                                                                                          <  2:3372412/91‑1 (MQ=255)
  tCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTa                                                                                        <  2:2270032/91‑1 (MQ=255)
    aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                                      <  2:744659/91‑3 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                    >  1:1398890/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                    <  1:2855938/91‑1 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                    <  1:2855836/91‑1 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttt                                                                                <  1:2392659/91‑1 (MQ=38)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                                <  1:290332/91‑1 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                                >  2:1749841/1‑91 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                                <  1:2343081/91‑1 (MQ=255)
            tttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttttt                                                                              <  2:45115/91‑3 (MQ=37)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                            >  1:366596/1‑91 (MQ=255)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                            >  1:1250370/1‑91 (MQ=255)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                            <  2:2906846/91‑1 (MQ=255)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                            >  2:1919316/1‑91 (MQ=255)
              tGATATAAAGGTGTTTGTTATAAATCACTTGGCATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGTTtttat                                                                            <  2:2600321/91‑2 (MQ=21)
                atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                          <  1:2814374/91‑1 (MQ=255)
                atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                          <  2:158894/91‑1 (MQ=255)
                 tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                         <  2:3039809/91‑1 (MQ=255)
                   taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                       >  2:184749/1‑91 (MQ=255)
                   taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                       <  1:2568850/91‑1 (MQ=255)
                    aaaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAcc                                                                      <  2:1352545/91‑1 (MQ=255)
                       ggTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAt                                                                   <  1:3537128/91‑1 (MQ=255)
                          gtttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCa                                                                >  1:359019/1‑91 (MQ=255)
                          gtttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCa                                                                <  1:380014/91‑1 (MQ=255)
                           tttgtttTAAATCACTTGGCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAg                                                               <  1:2440828/91‑1 (MQ=255)
                             tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                             <  1:1843001/91‑1 (MQ=255)
                              gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                            >  1:2781441/1‑91 (MQ=255)
                                tttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcg                                                          <  1:3061295/91‑1 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                         <  1:2756997/91‑1 (MQ=255)
                                   aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                       <  1:16156/91‑1 (MQ=255)
                                   aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                       >  2:3868369/1‑91 (MQ=255)
                                   aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                       >  2:3018843/1‑91 (MQ=255)
                                    aaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGtt                                                      <  1:209502/91‑1 (MQ=255)
                                      tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                    <  1:1467438/91‑1 (MQ=255)
                                      tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                    >  2:692250/1‑91 (MQ=255)
                                      tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                    <  2:1857069/91‑1 (MQ=255)
                                         cTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGt                                                 <  1:3401468/91‑1 (MQ=255)
                                           tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                               <  2:143650/91‑1 (MQ=255)
                                            gACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatt                                              >  1:434129/1‑91 (MQ=255)
                                             aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                             <  1:1229351/91‑1 (MQ=255)
                                              cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                            <  2:2936064/91‑1 (MQ=255)
                                              cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                            >  2:657800/1‑91 (MQ=255)
                                               aTTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATc                                           <  1:1308702/91‑1 (MQ=255)
                                                 tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                         <  1:1752070/91‑1 (MQ=255)
                                                 tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                         <  1:526449/91‑1 (MQ=255)
                                                   aaaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCttt                                       <  2:3451724/91‑1 (MQ=255)
                                                    aaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtttt                                      >  1:3140892/1‑91 (MQ=255)
                                                    aaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtttt                                      >  2:3320798/1‑91 (MQ=255)
                                                       ccccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTAcc                                   >  2:3705507/2‑91 (MQ=255)
                                                        cccccTTAAAATTATAAAAACAGGATAGCAATAGCTACAATCTTTATAAATCACCAATTGATGAACGCGTTACAGTATTATCCTTTTGCCt                                  >  2:2414271/1‑91 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 <  1:1298946/91‑1 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 <  1:808405/91‑1 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 <  1:1298255/91‑1 (MQ=255)
                                                          cccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTga                                <  2:2317239/91‑1 (MQ=255)
                                                            cTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTgaga                              <  1:2142837/91‑1 (MQ=255)
                                                              tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGt                            >  2:1445004/1‑91 (MQ=255)
                                                                aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                          <  1:2332125/91‑1 (MQ=255)
                                                                aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                          <  1:1733762/91‑1 (MQ=255)
                                                                 aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                         <  1:843805/91‑1 (MQ=255)
                                                                 aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                         >  1:3571860/1‑91 (MQ=255)
                                                                  aTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                        >  1:2555886/1‑91 (MQ=255)
                                                                   ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                       <  1:2521778/91‑1 (MQ=255)
                                                                    tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                      >  1:1133919/1‑91 (MQ=255)
                                                                    tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                      >  1:2035014/1‑91 (MQ=255)
                                                                    tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                      <  1:3169070/91‑1 (MQ=255)
                                                                     ataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaa                     >  1:1698498/1‑91 (MQ=255)
                                                                      taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                    >  1:897482/1‑91 (MQ=255)
                                                                      taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                    >  1:1009392/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)
                                                                        aaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa                  >  2:2157385/1‑91 (MQ=255)
                                                                         aaaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATaaa                 >  1:1350604/1‑91 (MQ=38)
                                                                          aaCAGGTTAGCGATTGCTATCCTGTTTATAAAGCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGaaataaat                >  2:2165038/1‑91 (MQ=38)
                                                                          aaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt                >  1:1840146/1‑91 (MQ=38)
                                                                           aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt               >  1:777805/1‑91 (MQ=38)
                                                                           aCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt               <  2:795115/91‑1 (MQ=255)
                                                                           aCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt               <  1:56530/91‑1 (MQ=255)
                                                                             aGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTAc             >  1:702428/1‑91 (MQ=39)
                                                                              ggATAGCAATCGCTACCCTGTTTATAAATCACCAGTGCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACt            <  2:149733/91‑1 (MQ=255)
                                                                               gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt           >  2:2093267/1‑91 (MQ=255)
                                                                               gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt           <  1:1981778/91‑1 (MQ=255)
                                                                                aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          <  2:900672/91‑1 (MQ=255)
                                                                                aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          <  2:3229348/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         <  2:1333074/91‑1 (MQ=255)
                                                                                 tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc         <  1:991289/91‑1 (MQ=255)
                                                                                  aGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa        >  1:3394266/1‑91 (MQ=255)
                                                                                   gCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAc       >  2:2380724/1‑91 (MQ=255)
                                                                                   gCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAc       >  2:342903/1‑91 (MQ=255)
                                                                                    cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc      <  2:867168/91‑1 (MQ=255)
                                                                                    cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc      <  1:3290513/91‑1 (MQ=255)
                                                                                      aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGTACCtt    >  2:3145887/1‑91 (MQ=255)
                                                                                      aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt    >  2:430871/1‑91 (MQ=255)
                                                                                      aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt    >  2:3522197/1‑91 (MQ=255)
                                                                                       tCGCTACCCTGTTTCTAAATCACGGAGTCAGGAACGCGTTACAGTATTACTCTTGTGCCTGAGAGTTTCGTAATAAATTACTTCCACCTTc   <  2:3605434/91‑1 (MQ=255)
                                                                                        cGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATCAATTACTTGCACCTTCg  >  1:3322165/1‑91 (MQ=255)
                                                                                        |                                                                                          
TATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCG  >  USA300TCH1516_ALE/1653247‑1653425

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