Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,044,663 T→C 9.6% intergenic (+12/+34) yfmC_1 → / ← USA300HOU_RS05165 Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,044,6630TC 9.6% 738.7 / 48.2 247intergenic (+12/+34)yfmC_1/USA300HOU_RS05165Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (124/99);  new base C (11/13);  total (135/112)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

GTGGTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAT  >  USA300TCH1516_ALE/1044573‑1044753
                                                                                          |                                                                                          
gTGGTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCt                                                                                            <  1:2954779/91‑1 (MQ=255)
gTGGTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCt                                                                                            <  1:721396/91‑1 (MQ=255)
gTGGTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCt                                                                                            >  2:366703/1‑91 (MQ=255)
   gTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACt                                                                                         <  2:1088245/91‑1 (MQ=255)
   gTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACt                                                                                         >  1:2239220/1‑91 (MQ=255)
   gTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACt                                                                                         >  2:431408/1‑91 (MQ=255)
   gTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACt                                                                                         <  2:2418207/91‑1 (MQ=255)
   gTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACt                                                                                         >  1:1798981/1‑91 (MQ=255)
    gtGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTc                                                                                        <  2:2540618/90‑1 (MQ=255)
     tGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCg                                                                                       >  1:2593198/1‑91 (MQ=255)
      gAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGt                                                                                      <  2:266088/91‑1 (MQ=255)
      gAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGt                                                                                      <  2:3554295/91‑1 (MQ=255)
      gAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGt                                                                                      <  1:2846313/91‑1 (MQ=255)
         aTCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTaaa                                                                                   >  1:22929/1‑91 (MQ=255)
          tCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAAt                                                                                  >  2:3877663/1‑91 (MQ=255)
          tCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAAt                                                                                  >  2:3458767/1‑91 (MQ=255)
             aGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCat                                                                               >  2:3082155/1‑91 (MQ=255)
                 ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                           >  2:2750235/1‑91 (MQ=255)
                 ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                           >  1:1272041/1‑91 (MQ=255)
                   aTTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAaga                                                                         <  2:882543/91‑1 (MQ=255)
                     tATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGt                                                                       <  1:1102836/91‑1 (MQ=25)
                       tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                     >  1:2878686/1‑91 (MQ=255)
                        cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                    >  1:1780156/1‑91 (MQ=255)
                        cGCAAGAGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGtagc                                                                    <  2:156485/91‑4 (MQ=18)
                         gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                   >  2:941040/1‑91 (MQ=255)
                         gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                   <  1:858864/91‑1 (MQ=255)
                           aaGTGAAAGTATGGCAGAAGATTTAGAAACAATGGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                 <  1:483634/91‑1 (MQ=255)
                             gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGt                                                               <  2:836467/91‑1 (MQ=17)
                             gTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGt                                                               <  2:122784/91‑1 (MQ=17)
                              tGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTa                                                              <  1:2890280/91‑1 (MQ=255)
                               gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                             >  1:317578/1‑91 (MQ=255)
                                aaaGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtt                                                            >  1:3543564/1‑91 (MQ=255)
                                 aaGTATGGCAGAAGACTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGCGGCGCTGTAttt                                                           >  2:3645666/1‑91 (MQ=255)
                                  aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                          >  2:2695378/1‑91 (MQ=255)
                                  aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                          >  1:1665485/1‑91 (MQ=255)
                                  aGTATGGAAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                          <  2:3106545/91‑1 (MQ=255)
                                   gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                         >  2:2103518/1‑91 (MQ=255)
                                   gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                         >  2:2225589/1‑91 (MQ=255)
                                    tATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTa                                                        <  1:1406438/91‑1 (MQ=255)
                                      tGGCAGAAGATTTAGAAAAAATTGCGGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                      <  2:3829762/91‑1 (MQ=255)
                                         cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                   >  1:1998588/1‑91 (MQ=255)
                                         cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                   <  2:1720304/91‑1 (MQ=255)
                                         cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                   <  2:1860672/91‑1 (MQ=255)
                                          agaagaGTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATg                                                  <  2:3264871/91‑1 (MQ=255)
                                              gaTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTtta                                              <  1:1461100/91‑1 (MQ=255)
                                               aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                             <  1:3331038/91‑1 (MQ=255)
                                               aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                             >  2:3869249/1‑91 (MQ=255)
                                                 ttAGAAAAAATTGCTGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                           <  1:1181649/91‑1 (MQ=255)
                                                 ttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                           <  1:813329/91‑1 (MQ=255)
                                                  tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                          <  1:999331/91‑1 (MQ=255)
                                                  tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                          <  1:3056194/91‑1 (MQ=255)
                                                  tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                          >  1:77490/1‑91 (MQ=255)
                                                  tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                          >  2:280505/1‑91 (MQ=255)
                                                   aGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATc                                         >  2:1009729/1‑91 (MQ=255)
                                                    gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                        <  2:3736134/91‑1 (MQ=255)
                                                    gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                        <  2:2257837/91‑1 (MQ=255)
                                                     aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                       <  1:1002222/91‑1 (MQ=255)
                                                      aaaaaTTGCATAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGt                                      <  2:1202913/91‑1 (MQ=255)
                                                       aaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATAATATAAGAGTGGCGCTGTATTTTTAATTATGCTTATTATCTGTc                                     >  1:1118615/1‑91 (MQ=255)
                                                         aaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCgg                                   <  1:1321493/91‑1 (MQ=255)
                                                          aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                  <  1:3583245/91‑1 (MQ=255)
                                                          aTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGa                                  <  1:1726087/91‑1 (MQ=255)
                                                            tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGCAtg                                >  1:3496234/1‑91 (MQ=255)
                                                              cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                              >  1:1072183/1‑91 (MQ=255)
                                                              cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                              >  2:582500/1‑91 (MQ=255)
                                                              cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                              <  1:3656924/91‑1 (MQ=255)
                                                              cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                              >  1:23701/1‑91 (MQ=255)
                                                              cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGAAGCGCTGTATTTTTAATTATGTTTATTATATGTCGGAtgtg                              >  2:20330/1‑91 (MQ=255)
                                                               aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                             >  2:2481024/1‑91 (MQ=255)
                                                               aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                             >  1:1601271/1‑91 (MQ=255)
                                                                gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                            >  1:2634149/1‑91 (MQ=255)
                                                                gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                            <  1:520091/91‑1 (MQ=255)
                                                                    aaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                        <  2:789093/91‑1 (MQ=255)
                                                                    aaGCAAAATAAAAATAAAGCGCTACTCGTAAATCATATAAGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                        <  2:3474808/91‑1 (MQ=255)
                                                                      gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                      >  2:1495103/1‑91 (MQ=255)
                                                                       cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                     >  1:1535134/1‑91 (MQ=255)
                                                                       cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                     >  2:1837850/1‑91 (MQ=255)
                                                                       cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                     >  1:37880/1‑91 (MQ=255)
                                                                        aaaataaaaataCAGTGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTGTTATCTGTCGGATGTGATGATTTAcc                    >  1:1487927/1‑91 (MQ=255)
                                                                        aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                    >  1:2655127/1‑91 (MQ=255)
                                                                        aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                    <  2:2960378/91‑1 (MQ=255)
                                                                        aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                    >  2:163974/1‑91 (MQ=255)
                                                                         aattaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                   >  1:2957015/4‑91 (MQ=11)
                                                                         aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                   >  1:3379387/1‑91 (MQ=255)
                                                                         aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                   >  2:1920118/1‑91 (MQ=255)
                                                                          aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                  >  1:2166687/1‑91 (MQ=255)
                                                                          aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                  >  1:560801/1‑91 (MQ=255)
                                                                            taaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaa                <  1:2410113/91‑1 (MQ=255)
                                                                             aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa               <  2:3230649/91‑1 (MQ=255)
                                                                              aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg              <  1:3086149/91‑1 (MQ=255)
                                                                              aaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAg              <  2:1998588/91‑1 (MQ=255)
                                                                               aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt             >  2:2864852/1‑91 (MQ=255)
                                                                               aaataCAGCGCTACTCGTAAATCATATAAGAGCGGCGCTGTACCTCTCATTATGTTTATTCCCTGTCGGATGTGATGATTTACCTGAAAgc             >  2:2052286/1‑90 (MQ=255)
                                                                                aataCAGCGCTACTCGTAAATCATGTAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt            >  1:2767664/1‑91 (MQ=255)
                                                                                 ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt           >  2:490980/1‑91 (MQ=255)
                                                                                  taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTa          >  1:1824279/1‑91 (MQ=17)
                                                                                   aCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAt         >  1:401063/1‑91 (MQ=255)
                                                                                       cgcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCg     >  2:3062353/1‑91 (MQ=255)
                                                                                        gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa    <  1:923581/91‑1 (MQ=255)
                                                                                        gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa    <  2:1576102/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTGGTTGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAT  >  USA300TCH1516_ALE/1044573‑1044753

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: