Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 419,955 C→T 7.7% intergenic (+173/+93) rpsR → / ← USA300HOU_RS01950 30S ribosomal protein S18/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE419,9550CT 7.7% 715.1 / 31.9 221intergenic (+173/+93)rpsR/USA300HOU_RS0195030S ribosomal protein S18/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (100/104);  new base T (15/2);  total (115/106)
Fisher's exact test for biased strand distribution p-value = 1.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 228 total aligned reads displayed.

ATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATA  >  USA300TCH1516_ALE/419865‑420045
                                                                                          |                                                                                          
atgatgTGAATGTTTTGTGTTCGGAGTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTAc                                                                                            <  1:3599269/91‑1 (MQ=255)
 tgatgTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACg                                                                                           >  1:1024036/1‑91 (MQ=255)
 tgatgTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACg                                                                                           <  1:3658082/91‑1 (MQ=255)
 tgatgTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACg                                                                                           >  1:966537/1‑91 (MQ=255)
     gtgAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCAtt                                                                                       >  2:3774046/1‑91 (MQ=255)
      tgAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCATAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCAttt                                                                                      <  2:2990437/91‑1 (MQ=255)
      tgAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCAttt                                                                                      <  2:1577218/91‑1 (MQ=255)
       gAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCAtttt                                                                                     <  2:1993122/91‑1 (MQ=255)
        aaTGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTta                                                                                    >  1:2060790/1‑91 (MQ=255)
        aaTGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTta                                                                                    <  1:1492078/91‑1 (MQ=255)
          tGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaa                                                                                  <  2:1249037/91‑1 (MQ=255)
           gTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtattt                                                                                 <  1:451293/91‑4 (MQ=255)
           gTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaat                                                                                 >  1:2387186/1‑91 (MQ=255)
            ttttGTGTTAGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtattta                                                                                <  2:3340096/91‑5 (MQ=255)
             tttGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTtaaataa                                                                               >  1:11291/1‑91 (MQ=255)
                gtgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAAt                                                                            >  2:3292558/1‑91 (MQ=255)
                gtgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAAt                                                                            >  2:936453/1‑91 (MQ=255)
                 tgtTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATg                                                                           <  2:1540755/91‑1 (MQ=255)
                 tgctCGGCATTTGCACAAAGATATGTTTATATTGAAAAAATAATATGAATTAAAAAAAAAAAAAAAAGAACTACGCATTTTAAATAAAATg                                                                           <  2:1542021/88‑1 (MQ=255)
                   ttCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTATGCATTTTATTTAAAATGcg                                                                         <  1:2270597/91‑2 (MQ=38)
                     cGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTATTTAAAATGcgta                                                                       <  1:706570/91‑4 (MQ=255)
                      ggAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAg                                                                      <  1:1565701/91‑1 (MQ=255)
                         aTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGctc                                                                   <  1:290617/91‑1 (MQ=255)
                         aTTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGctc                                                                   <  2:2764986/91‑1 (MQ=255)
                         aTTTGCACAAAGATATGTTGATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGctc                                                                   <  2:912536/91‑1 (MQ=255)
                            tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                <  1:1721418/91‑1 (MQ=255)
                            tGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttc                                                                >  2:800129/1‑91 (MQ=255)
                             gCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttct                                                               <  2:2884212/91‑1 (MQ=255)
                             gCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCtcttct                                                               <  1:1463245/91‑1 (MQ=255)
                              cacaAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTcttctt                                                              <  1:526215/91‑1 (MQ=255)
                               acaAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCttt                                                             >  2:3280251/1‑91 (MQ=255)
                                   aGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCt                                                         >  1:276972/1‑91 (MQ=255)
                                   aGATATGTTTATATTGCAAAAACAATACGAATTTCGATGCATACTCAAAGAACTACGCATTTTAAATAAAATACATAGCTCTTCTTTTTcc                                                         >  2:2202368/1‑90 (MQ=255)
                                      tatGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGc                                                      <  1:2737950/91‑1 (MQ=255)
                                      tatGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGc                                                      <  2:3072418/91‑1 (MQ=255)
                                       atGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCa                                                     >  2:2561189/1‑91 (MQ=255)
                                       atGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCa                                                     >  2:538084/1‑91 (MQ=255)
                                       atGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCa                                                     <  1:554080/91‑1 (MQ=255)
                                        tGTTTATATTGCCAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCAt                                                    >  1:2969273/1‑91 (MQ=255)
                                         gTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATa                                                   >  1:1812641/1‑91 (MQ=255)
                                          tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                  >  1:1876635/1‑91 (MQ=255)
                                          tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                  >  2:1952076/1‑91 (MQ=255)
                                          tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                  <  1:626672/91‑1 (MQ=255)
                                          tttATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATAc                                                  <  1:1270050/91‑1 (MQ=255)
                                             atatTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGaa                                               <  1:2281122/91‑1 (MQ=255)
                                                  gCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaa                                          >  1:2570405/1‑91 (MQ=255)
                                                  gCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaa                                          >  2:462409/1‑91 (MQ=255)
                                                   cAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaaa                                         >  1:3456574/1‑91 (MQ=255)
                                                   cAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTaaaa                                         >  2:1133741/1‑91 (MQ=255)
                                                      aaaTAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATaa                                      <  2:635286/91‑1 (MQ=255)
                                                       aataatATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAAc                                     <  2:731709/91‑1 (MQ=255)
                                                       aataatATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAAc                                     <  1:3149263/91‑1 (MQ=255)
                                                       aataatATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAAc                                     <  2:600661/91‑1 (MQ=255)
                                                           atatGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgc                                 >  1:3534952/1‑91 (MQ=255)
                                                           atatGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgc                                 <  1:2557464/91‑1 (MQ=255)
                                                            tatGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                >  1:1438406/1‑91 (MQ=255)
                                                             atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                               >  2:1044239/1‑91 (MQ=255)
                                                             atGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCga                               >  2:856943/1‑91 (MQ=255)
                                                                aaTATAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAAGAACTCGCGAGAc                            <  1:3001129/91‑1 (MQ=255)
                                                                 aTTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGAcc                           >  1:1811260/1‑91 (MQ=255)
                                                                 aTTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGAcc                           >  2:2428402/1‑91 (MQ=255)
                                                                  tttAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                             <  1:2547630/88‑1 (MQ=255)
                                                                  tttAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCgaga                             >  2:2547630/1‑88 (MQ=255)
                                                                  tttAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCt                          <  1:2095263/91‑1 (MQ=255)
                                                                   ttAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCta                         <  2:189799/91‑1 (MQ=255)
                                                                      gATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATaa                      <  1:782419/91‑1 (MQ=255)
                                                                       aTGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAg                     <  2:645666/91‑1 (MQ=255)
                                                                        tGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGt                    <  2:309792/91‑1 (MQ=255)
                                                                        tGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGt                    <  2:614620/91‑1 (MQ=255)
                                                                        tGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGt                    >  2:3026900/1‑91 (MQ=255)
                                                                        tGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGt                    >  1:1638152/1‑91 (MQ=255)
                                                                          cATAAATAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtct                  >  1:3015079/1‑91 (MQ=255)
                                                                          cATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtct                  <  2:3529846/91‑1 (MQ=255)
                                                                           aTAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTcgcg                                <  1:59523/76‑1 (MQ=255)
                                                                           aTAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtctc                 <  2:922765/91‑1 (MQ=255)
                                                                            tAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                <  1:385855/91‑1 (MQ=255)
                                                                            tAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                <  2:2900764/91‑1 (MQ=255)
                                                                            aaaaagaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                >  1:490005/2‑91 (MQ=14)
                                                                            aaaaagaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTctct                >  2:3521828/2‑91 (MQ=14)
                                                                             aaaagaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCtt               >  2:2196781/1‑91 (MQ=17)
                                                                             aaaaaaaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCtt               <  2:1788690/91‑1 (MQ=255)
                                                                              aaaaaaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCttt              <  1:3447531/91‑1 (MQ=255)
                                                                                 gaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGCATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt           >  1:1493460/2‑91 (MQ=18)
                                                                                 gaaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCt           >  1:1533188/2‑91 (MQ=21)
                                                                                  aaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc          >  1:1159846/1‑91 (MQ=18)
                                                                                  aaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTTTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc          >  1:1159881/1‑91 (MQ=18)
                                                                                  aaGAGCTATGCATTTTATTTAAAATGCGTAGTTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc          >  1:1376987/1‑91 (MQ=21)
                                                                                  aaGAGCTATGCATTTTATTTAAAATGCGTAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGCGACCTATAAGTCTCTTTCCTc          >  2:3660004/1‑91 (MQ=25)
                                                                                  aaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc          >  2:1118061/1‑91 (MQ=255)
                                                                                  aaGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTc          >  1:806001/1‑91 (MQ=255)
                                                                                    gAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtc                   >  1:3200539/1‑80 (MQ=255)
                                                                                    gAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGtc                   <  2:3200539/80‑1 (MQ=255)
                                                                                    gAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCAc        >  1:1419505/1‑91 (MQ=255)
                                                                                     aaCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACt       <  2:2535930/91‑1 (MQ=255)
                                                                                     aaCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACt       <  1:3292558/91‑1 (MQ=255)
                                                                                      aCTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTa      >  1:506399/1‑91 (MQ=255)
                                                                                       cTACGCATTTTAAATAAAATGCATAGTTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAg     >  2:2219298/1‑91 (MQ=255)
                                                                                         aCGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCCTACGAATTAAAATAACTCGCGAGACCGATAAGTCTCTTTCCTCACTaact   >  2:3603747/1‑88 (MQ=255)
                                                                                          cGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATa  >  1:2126060/1‑91 (MQ=255)
                                                                                          cGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATa  >  1:2126198/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAACTACGCATTTTAAATAAAATGCATAGCTCTTCTTTTTCTTGCATACGAATTAAAATAACTCGCGAGACCTATAAGTCTCTTTCCTCACTAGATA  >  USA300TCH1516_ALE/419865‑420045

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: