Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,327 A→G 8.8% intergenic (‑32/+127) gcvT ← / ← aroK Aminomethyltransferase/Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,3270AG 8.8% 600.2 / 37.1 205intergenic (‑32/+127)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base A (80/106);  new base G (8/10);  total (89/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 210 total aligned reads displayed.

TTGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTG  >  USA300TCH1516_ALE/1653238‑1653417
                                                                                         |                                                                                          
ttGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGAt                                                                                           <  2:754869/91‑1 (MQ=255)
 tGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATa                                                                                          >  1:1988861/1‑91 (MQ=255)
    aCCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCa                                                                                       >  1:1076796/1‑91 (MQ=255)
             aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                              <  2:744659/91‑3 (MQ=255)
             aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                              <  1:2881830/91‑3 (MQ=255)
             aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTatc                                                                              <  1:1797288/91‑3 (MQ=255)
               cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            <  1:2855836/91‑1 (MQ=255)
               cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            <  1:1250349/91‑1 (MQ=255)
               cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            <  1:2855938/91‑1 (MQ=255)
                aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTg                                                                           <  1:728596/91‑1 (MQ=255)
                aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTg                                                                           <  1:2691721/91‑1 (MQ=255)
                   aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttt                                                                        <  1:2392659/91‑1 (MQ=38)
                   aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        >  2:1749841/1‑91 (MQ=255)
                   aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        <  1:3665012/91‑1 (MQ=255)
                   aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        <  1:2343081/91‑1 (MQ=255)
                   aTTGTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        <  1:3134727/91‑1 (MQ=255)
                    ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGtttt                                                                       <  2:416652/91‑2 (MQ=37)
                    ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATGATAAAAACAGGATAGCAATCGCTACCCTGTTta                                                                       >  1:3180012/1‑91 (MQ=255)
                     tttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttttt                                                                      <  2:45115/91‑3 (MQ=37)
                      ttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTtata                                                                     >  1:3058110/1‑91 (MQ=255)
                       tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                    >  1:2990009/1‑91 (MQ=255)
                       tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                    >  2:1919316/1‑91 (MQ=255)
                       tGATATAAAGGTGTTTGTTATAAATCACTTGGCATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGTTtttat                                                                    <  2:2600321/91‑2 (MQ=21)
                        gATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaaa                                                                   <  1:2135402/91‑1 (MQ=255)
                         atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                  <  1:2814374/91‑1 (MQ=255)
                          tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                 >  2:3235170/1‑91 (MQ=255)
                          tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                 <  2:2277439/91‑1 (MQ=255)
                           ataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCa                                                                <  1:1207389/91‑1 (MQ=255)
                            taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                               <  1:2564777/91‑1 (MQ=255)
                                   gtttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCa                                                        >  2:3074127/1‑91 (MQ=255)
                                      tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAGTCTGTTTATAAATCACCAATTCAGGa                                                     >  2:3354732/1‑91 (MQ=255)
                                      tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                     <  1:1842948/91‑1 (MQ=255)
                                       gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                    >  1:2781441/1‑91 (MQ=255)
                                       gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                    <  2:320320/91‑1 (MQ=255)
                                          ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 >  2:797119/1‑91 (MQ=255)
                                           tAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgcg                                                >  1:3299578/1‑91 (MQ=255)
                                            aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                               <  1:16156/91‑1 (MQ=255)
                                               tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                            <  1:3534964/91‑1 (MQ=255)
                                               tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                            <  1:1470432/91‑1 (MQ=255)
                                               tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                            <  2:1857034/91‑1 (MQ=255)
                                                cACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACa                                           >  2:2359940/1‑91 (MQ=255)
                                                cACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACa                                           >  1:3547591/1‑91 (MQ=255)
                                                  cTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGt                                         <  1:3401468/91‑1 (MQ=255)
                                                    tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                       <  2:143650/91‑1 (MQ=255)
                                                    tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                       <  1:2423866/91‑1 (MQ=255)
                                                     gACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatt                                      <  1:2063322/91‑1 (MQ=255)
                                                     gACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatt                                      >  1:434129/1‑91 (MQ=255)
                                                      aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                     <  2:710712/91‑1 (MQ=255)
                                                      aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAGTCAGGAACGCGTTACAGtatta                                     >  2:2947329/1‑91 (MQ=255)
                                                       cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                    <  2:2936010/91‑1 (MQ=255)
                                                       cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                    >  2:657800/1‑91 (MQ=255)
                                                        aTTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATc                                   <  2:2156268/91‑1 (MQ=255)
                                                        aTTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATc                                   <  1:1308702/91‑1 (MQ=255)
                                                          tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                 <  1:526984/91‑1 (MQ=255)
                                                          tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                 <  1:1752070/91‑1 (MQ=255)
                                                            aaaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCttt                               <  2:3451724/91‑1 (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                         <  2:2781441/91‑1 (MQ=255)
                                                                  ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         <  1:808405/91‑1 (MQ=255)
                                                                  ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         >  2:1807171/1‑91 (MQ=255)
                                                                  ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         <  1:134164/91‑1 (MQ=255)
                                                                  ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         <  1:1298946/91‑1 (MQ=255)
                                                                  ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         <  1:1298476/91‑1 (MQ=255)
                                                                     cTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTgaga                      <  1:2142837/91‑1 (MQ=255)
                                                                      ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                     <  2:3027661/91‑1 (MQ=255)
                                                                       tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGt                    >  2:1445004/1‑91 (MQ=255)
                                                                         aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                  <  2:386545/91‑1 (MQ=255)
                                                                         aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                  <  1:2332125/91‑1 (MQ=255)
                                                                         aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                  <  1:1733762/91‑1 (MQ=255)
                                                                         aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                  >  1:2780774/1‑91 (MQ=255)
                                                                          aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                 <  1:843805/91‑1 (MQ=255)
                                                                           aTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGACAGTTTCg                >  2:3442741/1‑91 (MQ=255)
                                                                            ttATAAAAACAGGATAGCCATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt               <  1:2521752/91‑1 (MQ=255)
                                                                            ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt               >  1:630409/1‑91 (MQ=255)
                                                                            ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt               >  1:1453535/1‑91 (MQ=255)
                                                                             tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta              <  1:3169070/91‑1 (MQ=255)
                                                                              ataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaa             >  2:1660647/1‑91 (MQ=255)
                                                                              ataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaa             >  1:1698498/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)
                                                                                aaaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaata           <  2:2555886/91‑1 (MQ=255)
                                                                                 taaaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa          >  2:2516290/2‑91 (MQ=37)
                                                                                 aaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa          <  2:1133919/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)
                                                                                   aaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt        <  2:10207/91‑1 (MQ=255)
                                                                                    aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt       <  2:216851/91‑1 (MQ=38)
                                                                                    aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt       >  1:777805/1‑91 (MQ=38)
                                                                                     cAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa      >  1:1912169/1‑91 (MQ=255)
                                                                                     cAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa      <  2:1564755/91‑1 (MQ=255)
                                                                                       ggATAGCAATCGCTACCCTGTTTATAAATCACCAGTGCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACt    <  2:149733/91‑1 (MQ=255)
                                                                                        gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt   >  2:2093267/1‑91 (MQ=255)
                                                                                         aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg  <  2:1250370/91‑1 (MQ=255)
                                                                                         aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg  <  2:900672/91‑1 (MQ=255)
                                                                                         |                                                                                          
TTGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTG  >  USA300TCH1516_ALE/1653238‑1653417

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: