Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,093 T→A 19.6% intergenic (+87/‑275) yugI_2 → / → USA300HOU_RS04740 General stress protein 13/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE950,0930TA19.6% 387.3 / 87.5 148intergenic (+87/‑275)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (49/70);  new base A (14/15);  total (63/85)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

Only 100 of 156 total aligned reads displayed.

CTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAT  >  USA300TCH1516_ALE/950003‑950183
                                                                                          |                                                                                          
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            <  1:1790144/91‑1 (MQ=255)
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            >  1:3827330/1‑91 (MQ=255)
  aaaGGAACAGATAGATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                          >  1:3141567/1‑91 (MQ=255)
  aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                          <  1:2573180/91‑1 (MQ=255)
    aGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTtta                                                                                        >  1:3559106/1‑91 (MQ=255)
       aaCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATttaata                                                                                     <  1:1951101/91‑4 (MQ=38)
        aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAg                                                                                    >  1:2705436/1‑91 (MQ=255)
        aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                    <  1:1327351/91‑1 (MQ=38)
        aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                    <  1:506923/91‑1 (MQ=38)
           gATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATa                                                                                 <  1:776276/91‑1 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                >  1:1123979/1‑91 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                >  1:1456225/1‑91 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                >  1:1456394/1‑91 (MQ=255)
             tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                               <  1:2633108/91‑1 (MQ=255)
              aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTa                                                                              <  1:2810610/91‑1 (MQ=255)
              aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTa                                                                              <  1:3731722/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             >  1:3674747/1‑91 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  1:1430157/91‑1 (MQ=255)
                aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                            >  1:1039448/1‑91 (MQ=255)
                 tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                           >  1:2819603/1‑91 (MQ=255)
                 tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                           <  1:1546167/91‑1 (MQ=255)
                  cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtct                                                                          <  1:1410216/91‑1 (MQ=25)
                   gTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtctc                                                                         <  1:1864260/91‑1 (MQ=25)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        >  1:3451376/1‑91 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        <  1:3903179/91‑1 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        <  1:2320094/91‑1 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        <  1:2319964/91‑1 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATctct                                                                        <  1:1925366/91‑1 (MQ=25)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                      >  1:1823617/1‑91 (MQ=255)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                      >  1:3387485/1‑91 (MQ=255)
                      ccGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATCTCttt                                                                      <  1:1836118/91‑1 (MQ=21)
                       cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                     <  1:3713297/91‑1 (MQ=255)
                       cGAAAATCATACAAAGGGTCTGAAATGAAAGGTTCTTAGACTATGAAAGAGATTAGTGTCTATTAAATTTTATTAGATACTAATCTCtttt                                                                     <  1:3310122/91‑1 (MQ=255)
                           aaTCATACAACGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTc                                                                 <  1:1898588/91‑1 (MQ=255)
                            aTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCt                                                                >  1:1088536/1‑91 (MQ=255)
                              cATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTAc                                                              <  1:479753/91‑1 (MQ=255)
                               aTACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACg                                                             <  1:3984164/91‑1 (MQ=255)
                                 aCAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGAt                                                           <  1:2415672/91‑1 (MQ=255)
                                    aaGGGTCTGAAGTGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                        <  1:3752489/91‑1 (MQ=255)
                                    aaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAAc                                                        >  1:459778/1‑91 (MQ=255)
                                      gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                      <  1:1892388/91‑1 (MQ=255)
                                          cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                  >  1:3871007/1‑91 (MQ=255)
                                           tgaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAatat                                                 <  1:2653472/91‑1 (MQ=255)
                                             aaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtga                                               >  1:1020030/1‑91 (MQ=255)
                                              aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                              >  1:2090366/1‑91 (MQ=255)
                                                tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                            <  1:3942952/91‑1 (MQ=255)
                                                tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                            <  1:621606/91‑1 (MQ=255)
                                                tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTGTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                            <  1:3302629/91‑1 (MQ=255)
                                                  aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                          <  1:1256949/91‑1 (MQ=255)
                                                  aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                          <  1:307382/91‑1 (MQ=255)
                                                   aaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATgattgatt                                         <  1:590255/91‑1 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        <  1:3026607/91‑1 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        <  1:3026594/91‑1 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        >  1:2716737/1‑91 (MQ=255)
                                                        tCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAtt                                    >  1:3357754/1‑91 (MQ=255)
                                                        tCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAtt                                    <  1:308619/91‑1 (MQ=255)
                                                           tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                                 <  1:1115607/91‑1 (MQ=255)
                                                              aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:2607350/91‑1 (MQ=255)
                                                              aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:535744/91‑1 (MQ=255)
                                                              aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:2857340/91‑1 (MQ=255)
                                                              aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:2852314/91‑1 (MQ=255)
                                                              aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:2518975/91‑1 (MQ=255)
                                                               cTATAAAAGTGATTAGTATCTATTAAATTTTATTAGATACTATTCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                             <  1:2513383/91‑1 (MQ=255)
                                                               cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                             >  1:1322987/1‑91 (MQ=255)
                                                                 ataAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTAc                           <  1:1507494/91‑1 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTACTCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         >  1:1379125/1‑91 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         <  1:570860/91‑1 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         <  1:3430016/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  1:2378833/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  1:1079546/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  1:2378418/91‑1 (MQ=255)
                                                                       gagaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATgg                     <  1:1497065/91‑1 (MQ=255)
                                                                        agaTTAGTAACTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTAGAAATGGt                    >  1:2080496/1‑91 (MQ=255)
                                                                         gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   <  1:1508598/91‑1 (MQ=255)
                                                                         gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   >  1:133083/1‑91 (MQ=255)
                                                                         gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   <  1:1234597/91‑1 (MQ=255)
                                                                          aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                  >  1:3650824/1‑91 (MQ=25)
                                                                          aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                  >  1:218184/1‑91 (MQ=25)
                                                                            tAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTaa                >  1:2355485/1‑91 (MQ=255)
                                                                             aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:3912994/1‑91 (MQ=255)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:2529825/1‑91 (MQ=37)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:3764901/1‑91 (MQ=37)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:3845382/1‑91 (MQ=37)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:462284/1‑91 (MQ=37)
                                                                               tATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGt             <  1:2393935/91‑1 (MQ=255)
                                                                                aTCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            <  1:4044651/91‑1 (MQ=255)
                                                                                 tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           <  1:1455320/91‑1 (MQ=255)
                                                                                 tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           <  1:329629/91‑1 (MQ=255)
                                                                                 tCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           >  1:2568361/1‑91 (MQ=37)
                                                                                  cTATGAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGAc          >  1:3280748/1‑91 (MQ=255)
                                                                                   tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa         >  1:1407739/1‑91 (MQ=255)
                                                                                   tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa         >  1:3389478/1‑91 (MQ=255)
                                                                                    aTTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACat        >  1:1962470/1‑91 (MQ=255)
                                                                                    aTTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACat        <  1:2894763/91‑1 (MQ=255)
                                                                                       aaaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATc     >  1:611364/1‑91 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  1:2713659/91‑1 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  1:2713830/91‑1 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAAGATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  1:4142480/91‑1 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  1:2823676/91‑1 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  1:3572456/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAT  >  USA300TCH1516_ALE/950003‑950183

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: