Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,044,668 G→T 9.5% intergenic (+17/+29) 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,6680GT 9.5% 802.7 / 53.8 239intergenic (+17/+29)yfmC_1/USA300HOU_RS05165Fe(3+)‑citrate‑binding protein YfmC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (119/96);  new base T (11/12);  total (130/109)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

TGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAT  >  USA300TCH1516_ALE/1044578‑1044757
                                                                                          |                                                                                         
 gAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGt                                                                                          <  1:2846334/91‑1 (MQ=255)
 gAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGt                                                                                          <  2:266088/91‑1 (MQ=255)
   aaTCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTaa                                                                                        <  2:3168266/91‑1 (MQ=255)
   aaTCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTaa                                                                                        >  2:1858263/1‑91 (MQ=255)
    aTCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTaaa                                                                                       >  1:22929/1‑91 (MQ=255)
     tCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAAt                                                                                      >  2:3877663/1‑91 (MQ=255)
      cAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATc                                                                                     >  1:1781920/1‑91 (MQ=255)
        aGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCat                                                                                   >  2:3082155/1‑91 (MQ=255)
        aGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCat                                                                                   >  1:1461648/1‑91 (MQ=255)
            ggATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATaa                                                                               >  1:1272041/1‑91 (MQ=255)
              aTTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAaga                                                                             <  2:882543/91‑1 (MQ=255)
                tATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTGAATCATATAAGAGt                                                                           <  1:1713154/91‑1 (MQ=255)
                 aTCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGAAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTg                                                                          <  2:2454386/91‑1 (MQ=255)
                  tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                         >  2:2208438/1‑91 (MQ=255)
                  tCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTgg                                                                         <  2:2339811/91‑1 (MQ=255)
                   cGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgc                                                                        >  1:613546/1‑91 (MQ=255)
                    gCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGgcg                                                                       >  2:941040/1‑91 (MQ=255)
                      aaGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCt                                                                     <  1:3526952/91‑1 (MQ=255)
                       aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                    >  2:2970441/1‑91 (MQ=255)
                       aGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTg                                                                    >  2:3431937/1‑91 (MQ=255)
                         tGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTa                                                                  <  1:894054/91‑1 (MQ=255)
                          gAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAt                                                                 >  1:317578/1‑91 (MQ=255)
                             aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                              >  1:151632/1‑91 (MQ=255)
                             aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                              >  2:2695378/1‑91 (MQ=255)
                             aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                              >  1:1664090/1‑91 (MQ=255)
                             aGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTAtttt                                                              <  2:2076859/91‑1 (MQ=17)
                             aGTATGGAAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAtttt                                                              <  2:3106545/91‑1 (MQ=255)
                              gTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTAttttt                                                             <  2:2500338/91‑1 (MQ=255)
                                 tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                          <  1:3082155/91‑1 (MQ=255)
                                 tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAAt                                                          >  1:1928552/1‑91 (MQ=255)
                                 tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTatt                                                          <  1:358/91‑3 (MQ=14)
                                 tGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTatt                                                          <  1:2396622/91‑3 (MQ=14)
                                    cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                       >  1:2623089/1‑91 (MQ=255)
                                    cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                       >  1:1576102/1‑91 (MQ=255)
                                    cAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTAt                                                       <  2:1720304/91‑1 (MQ=255)
                                     agaagaGTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATg                                                      <  2:3264871/91‑1 (MQ=255)
                                          aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                 >  2:3867494/1‑91 (MQ=255)
                                          aTTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttat                                                 >  2:3869249/1‑91 (MQ=255)
                                           tttAGAAAAAATTGCAGATAAAGCAGGTTAAAAATACAGTGCTACTCGTAAATCATAGAAGTGTGGCGCTGTAATTTTCATTATGTttatt                                                >  2:1523630/1‑91 (MQ=255)
                                            ttAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTttatta                                               >  1:3773094/1‑91 (MQ=255)
                                             tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                              <  1:3056194/91‑1 (MQ=255)
                                             tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                              >  1:2752620/1‑91 (MQ=255)
                                             tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                              <  1:50221/91‑1 (MQ=255)
                                             tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                              >  1:1090456/1‑91 (MQ=255)
                                             tAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTtattat                                              <  1:999331/91‑1 (MQ=255)
                                              aGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATc                                             >  2:1009729/1‑91 (MQ=255)
                                               gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                            <  1:2640798/91‑1 (MQ=255)
                                               gAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCt                                            <  2:3738408/91‑1 (MQ=255)
                                                aaaaaaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTg                                           <  1:1002222/91‑1 (MQ=255)
                                                 aaaaaTTGCATAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGt                                          <  2:1202913/91‑1 (MQ=255)
                                                    aaTTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCgg                                       <  1:1321493/91‑1 (MQ=255)
                                                      ttGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAt                                     <  1:3458767/91‑1 (MQ=255)
                                                       tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                    <  1:1422922/91‑1 (MQ=255)
                                                       tGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCATAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtg                                    <  2:2239220/91‑1 (MQ=255)
                                                        gCAGAAAAAGCAAAATAAAAATACAGCGCTACACGTAAATCATATAAGAGTGGCGCTGTATTTCTAATTATGTATATTATCTGTCAGAtgt                                   >  2:1226836/1‑91 (MQ=39)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  >  2:582500/1‑91 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  >  1:1072183/1‑91 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  <  1:115770/91‑1 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  >  1:23701/1‑91 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  <  1:818381/91‑1 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGAtgtg                                  <  1:3594039/91‑1 (MQ=255)
                                                         cAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGAAGCGCTGTATTTTTAATTATGTTTATTATATGTCGGAtgtg                                  >  2:20330/1‑91 (MQ=255)
                                                          aGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtga                                 >  2:834063/1‑91 (MQ=255)
                                                           gAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGtgat                                <  1:520091/91‑1 (MQ=255)
                                                              aaaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTgatgat                             <  1:2238579/91‑1 (MQ=255)
                                                              aaaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTAAATTATGTTTATTATCTGTCGGATGTgatgat                             >  1:3319406/1‑91 (MQ=255)
                                                               aaGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                            <  2:789093/91‑1 (MQ=255)
                                                               aaGCAAAATAAAAATAAAGCGCTACTCGTAAATCATATAAGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGAtt                            <  2:3474808/91‑1 (MQ=255)
                                                                 gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                          >  2:1495103/1‑91 (MQ=255)
                                                                 gCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTa                          <  2:2752620/91‑1 (MQ=255)
                                                                  cAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAc                         >  2:1837850/1‑91 (MQ=255)
                                                                   aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                        >  2:163974/1‑91 (MQ=255)
                                                                   aaaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTAcc                        <  1:3801/91‑1 (MQ=255)
                                                                    aattaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                       >  1:2957015/4‑91 (MQ=11)
                                                                    aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                       >  1:3379387/1‑91 (MQ=255)
                                                                    aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                       >  2:1920118/1‑91 (MQ=255)
                                                                    aaataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCt                       >  2:3697686/1‑91 (MQ=255)
                                                                     aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                      >  1:2166687/1‑91 (MQ=255)
                                                                     aataaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTg                      >  2:1955016/1‑91 (MQ=255)
                                                                      ttaaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGa                     >  2:456600/2‑91 (MQ=17)
                                                                       taaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaa                    >  2:3512021/1‑91 (MQ=255)
                                                                        aaaaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                   <  1:982298/91‑1 (MQ=255)
                                                                        aaaaataCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGaaa                   >  2:2860290/1‑91 (MQ=14)
                                                                          aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                 <  1:779190/91‑1 (MQ=255)
                                                                          aaataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCAGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGt                 >  2:2052248/1‑91 (MQ=255)
                                                                           aataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGtt                >  1:2767901/1‑91 (MQ=255)
                                                                            ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt               <  1:1776503/91‑1 (MQ=255)
                                                                            ataCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGttt               >  2:491083/1‑91 (MQ=255)
                                                                             taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTa              >  1:1824279/1‑91 (MQ=17)
                                                                             taCAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATATTTAATTATGTTTATTATCTGTCGGGTGTGATGATTTACCTGAAAGTTTa              >  2:3597761/1‑91 (MQ=17)
                                                                              aCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAt             >  1:401063/1‑91 (MQ=255)
                                                                              aCAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAt             >  1:2884257/1‑91 (MQ=255)
                                                                               cAGCGCCACTCTTATATGATTTACGAGTAGCGCTGTATTCTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTAtt            >  1:772451/1‑91 (MQ=18)
                                                                                 gcgcCACTCTTATATGATTTACGAGTAGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTc          >  1:1117954/1‑91 (MQ=18)
                                                                                  cgcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCg         >  2:3062353/1‑91 (MQ=255)
                                                                                   gcTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGa        <  2:1576101/91‑1 (MQ=255)
                                                                                     tACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAAt      <  1:2894173/91‑1 (MQ=255)
                                                                                      aCTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATa     <  1:15473/91‑1 (MQ=255)
                                                                                        tCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATaaa   >  2:1683994/1‑91 (MQ=255)
                                                                                         cGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAt  <  2:2878686/91‑1 (MQ=255)
                                                                                          |                                                                                         
TGAAATCAAGGGGGATTATCGCAAGTGAAAGTATGGCAGAAGATTTAGAAAAAATTGCAGAAAAAGCAAAATAAAAATACAGCGCTACTCGTAAATCATATAAGAGTGGCGCTGTATTTTTAATTATGTTTATTATCTGTCGGATGTGATGATTTACCTGAAAGTTTATTTCGAATAAAT  >  USA300TCH1516_ALE/1044578‑1044757

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: