Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,327 A→G 8.3% 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.3% 442.8 / 23.3 144intergenic (‑32/+127)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base A (66/66);  new base G (3/9);  total (69/75)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 148 total aligned reads displayed.

TTTGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTG  >  USA300TCH1516_ALE/1653237‑1653417
                                                                                          |                                                                                          
tttGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGa                                                                                            <  1:3065631/91‑1 (MQ=255)
tttGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGa                                                                                            >  1:2170665/1‑91 (MQ=255)
 ttGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGAt                                                                                           <  1:2718349/91‑1 (MQ=255)
   gCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAg                                                                                         >  1:48326/1‑91 (MQ=255)
   gCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAg                                                                                         >  1:2301570/1‑91 (MQ=255)
   gCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAg                                                                                         >  1:1868159/1‑91 (MQ=255)
    cACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGc                                                                                        <  1:1716148/91‑1 (MQ=255)
     aCCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCa                                                                                       >  1:2224233/1‑91 (MQ=255)
          tATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGc                                                                                  <  1:1047517/91‑1 (MQ=255)
           aTCAACATAATTTTGATATAATGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCt                                                                                 <  1:2578268/91‑1 (MQ=255)
           aTCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCt                                                                                 <  1:2442434/91‑1 (MQ=255)
            tCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTa                                                                                <  1:527832/91‑1 (MQ=255)
             cAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTat                                                                               <  1:488188/91‑2 (MQ=255)
             cAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAc                                                                               <  1:873771/91‑1 (MQ=255)
             cAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAc                                                                               <  1:3642757/91‑1 (MQ=255)
             cAACATAATATTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTat                                                                               <  1:3724244/91‑2 (MQ=255)
              aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAcc                                                                              >  1:2920391/1‑91 (MQ=255)
               aCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAccc                                                                             <  1:2155272/91‑1 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:3664048/1‑91 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:2788250/1‑91 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:468537/1‑91 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:2009801/1‑91 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:3816958/1‑91 (MQ=255)
                cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                            >  1:587182/1‑91 (MQ=255)
                 aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTg                                                                           <  1:1963208/91‑1 (MQ=38)
                   aaTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGtt                                                                         <  1:2606434/91‑1 (MQ=255)
                    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        <  1:3464345/91‑1 (MQ=255)
                    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        >  1:2033117/1‑91 (MQ=255)
                    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        >  1:2033368/1‑91 (MQ=255)
                    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTAGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                        <  1:3402986/91‑1 (MQ=255)
                      tttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTtat                                                                      <  1:2767880/91‑1 (MQ=255)
                       ttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTtata                                                                     <  1:1601802/91‑1 (MQ=255)
                       ttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTtata                                                                     >  1:2839161/1‑91 (MQ=255)
                        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGtttttat                                                                    <  1:1713917/91‑5 (MQ=25)
                        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                    >  1:1910062/1‑91 (MQ=255)
                        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                    >  1:1909843/1‑91 (MQ=255)
                         gATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaaa                                                                   >  1:1751045/1‑91 (MQ=255)
                          atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                  <  1:654819/91‑1 (MQ=255)
                           tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                 <  1:1626350/91‑1 (MQ=255)
                           tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                 <  1:1356832/91‑1 (MQ=255)
                            ataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATTAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCa                                                                <  1:3507595/91‑1 (MQ=255)
                             taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                               <  1:1796020/91‑1 (MQ=255)
                              aaaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAcc                                                              <  1:2454614/91‑1 (MQ=255)
                               aaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCa                                                             <  1:1126763/91‑1 (MQ=255)
                               aaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCa                                                             >  1:1829450/1‑91 (MQ=255)
                                aGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCaa                                                            >  1:1795288/1‑91 (MQ=255)
                                  gtgtTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAGCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAtt                                                          <  1:2697092/91‑1 (MQ=255)
                                  gtgtTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAtt                                                          <  1:2697116/91‑1 (MQ=255)
                                    gttagtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATACATCACCAATTCa                                                        <  1:3212109/91‑1 (MQ=255)
                                       tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                     >  1:2014224/1‑91 (MQ=255)
                                       tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                     >  1:3175456/1‑91 (MQ=255)
                                       tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                     >  1:4117803/1‑91 (MQ=255)
                                        gtttTAAATCACCTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCCATTCAGGaa                                                    <  1:3626759/91‑1 (MQ=255)
                                          tttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcg                                                  <  1:2080474/91‑1 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 <  1:1333678/91‑1 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 >  1:493278/1‑91 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 <  1:92558/91‑1 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 >  1:4047301/1‑91 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 >  1:589089/1‑91 (MQ=255)
                                           ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgc                                                 <  1:1152935/91‑1 (MQ=255)
                                              aaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGtt                                              >  1:2862993/1‑91 (MQ=255)
                                                   cTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGt                                         >  1:1496139/1‑91 (MQ=255)
                                                    ttGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGta                                        >  1:611283/1‑91 (MQ=255)
                                                    ttGACATTGAAATACCCCCTCAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTCTATAAATCACCAATTCAGGAACGCGTTACAGtt                                        >  1:2016705/1‑90 (MQ=255)
                                                     tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                       <  1:2087914/91‑1 (MQ=255)
                                                     tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                       <  1:817941/91‑1 (MQ=255)
                                                     tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                       <  1:832167/91‑1 (MQ=255)
                                                       aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                     <  1:3826168/91‑1 (MQ=255)
                                                           tGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                 <  1:3570457/91‑1 (MQ=255)
                                                            gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                >  1:2911654/1‑91 (MQ=255)
                                                            gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                >  1:4106708/1‑91 (MQ=255)
                                                            gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                >  1:2120156/1‑91 (MQ=255)
                                                               aTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTg                             >  1:4037653/1‑91 (MQ=255)
                                                               aTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTg                             >  1:2914576/1‑91 (MQ=255)
                                                                tACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGc                            >  1:775744/1‑91 (MQ=255)
                                                                  cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                          <  1:995491/91‑1 (MQ=255)
                                                                  cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                          >  1:828427/1‑91 (MQ=255)
                                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         >  1:3569011/1‑91 (MQ=255)
                                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         <  1:1538075/91‑1 (MQ=255)
                                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                         >  1:4010462/1‑91 (MQ=255)
                                                                    cccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTga                        >  1:3073565/1‑91 (MQ=255)
                                                                      cTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTgaga                      >  1:1965086/1‑91 (MQ=255)
                                                                       ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                     >  1:1853129/1‑91 (MQ=255)
                                                                       ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                     >  1:1852574/1‑91 (MQ=255)
                                                                       ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                     >  1:1366651/1‑91 (MQ=255)
                                                                        tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGt                    <  1:985691/91‑1 (MQ=255)
                                                                            aTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                >  1:3042725/1‑91 (MQ=255)
                                                                            aTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                >  1:623864/1‑91 (MQ=255)
                                                                             ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt               <  1:1392881/91‑1 (MQ=255)
                                                                                taAAAACAGGATAGCGATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat            >  1:3529543/1‑91 (MQ=255)
                                                                                 aaaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaata           <  1:2744654/91‑1 (MQ=255)
                                                                                  aaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa          <  1:2986074/91‑1 (MQ=255)
                                                                                    aaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt        <  1:3377819/91‑1 (MQ=255)
                                                                                      cAGGGTAGCGATTGCTATCCTGTTTTTAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa      >  1:612036/1‑91 (MQ=37)
                                                                                      cAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa      >  1:2339578/1‑91 (MQ=38)
                                                                                      cAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa      <  1:1074867/91‑1 (MQ=255)
                                                                                         gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt   >  1:1178665/1‑91 (MQ=255)
                                                                                         gATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt   >  1:3356106/1‑91 (MQ=255)
                                                                                          aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg  >  1:1778296/1‑91 (MQ=255)
                                                                                          aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg  >  1:469029/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTTGCACCTCTATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTG  >  USA300TCH1516_ALE/1653237‑1653417

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: