Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,093 T→A 16.8% 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,0930TA16.8% 599.3 / 102.5 220intergenic (+87/‑275)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (85/98);  new base A (24/13);  total (109/111)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 241 total aligned reads displayed.

CTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAT  >  USA300TCH1516_ALE/950003‑950183
                                                                                          |                                                                                          
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            <  1:3189086/91‑1 (MQ=255)
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            <  2:3554966/91‑1 (MQ=255)
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            <  2:2299998/91‑1 (MQ=255)
cTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                            <  2:1942872/91‑1 (MQ=255)
 tAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtt                                                                                           >  1:100958/1‑91 (MQ=255)
  aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAttt                                                                                          >  1:3325675/1‑91 (MQ=255)
   aaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAtttt                                                                                         <  1:447318/91‑1 (MQ=255)
        aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAg                                                                                    <  1:3749559/91‑1 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                >  2:3165656/1‑91 (MQ=255)
            aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                                <  1:2373653/91‑1 (MQ=37)
             tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                               <  2:3752731/91‑1 (MQ=255)
              aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTa                                                                              <  2:1666655/91‑1 (MQ=37)
               aaacGTACAGAAAATCATACAAATGATCTGAAATGAAAGTTTCTTAAATTATAATAGAGATTAGTATTTATTAAATTTTATTAGATACTaa                                                                             <  2:1215804/88‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             <  2:1336575/91‑1 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             >  1:912312/1‑91 (MQ=255)
               aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                             >  2:3342190/1‑91 (MQ=255)
               aaTCGTACAGAAAATTATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATATATTAAATTTTATTAGATATTaa                                                                             >  1:1225740/1‑91 (MQ=255)
                aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                            <  2:140971/91‑1 (MQ=255)
                aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAt                                                                            <  2:1476450/91‑1 (MQ=25)
                 tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                           <  1:3446993/91‑1 (MQ=255)
                 tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                           <  1:1452566/91‑1 (MQ=255)
                  cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGTCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTCAtct                                                                          <  1:3069923/91‑1 (MQ=255)
                   gTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtctc                                                                         <  1:737496/91‑1 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        <  1:1070822/91‑1 (MQ=255)
                    tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATctct                                                                        <  2:2248690/91‑1 (MQ=255)
                     aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAATCTCtt                                                                       <  1:837131/91‑1 (MQ=25)
                       cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                     >  1:2746067/1‑91 (MQ=255)
                        gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                                    <  1:255874/91‑1 (MQ=255)
                        gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                                    <  1:2512430/91‑1 (MQ=255)
                         aaaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTg                                                                   >  2:2487160/1‑91 (MQ=255)
                                tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                            >  1:3756644/1‑91 (MQ=255)
                                  cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                          >  2:935531/1‑91 (MQ=255)
                                        gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                                    >  1:3884510/1‑91 (MQ=255)
                                        gTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTaa                                                    <  2:2444637/91‑1 (MQ=255)
                                         tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                   >  2:1442164/1‑91 (MQ=255)
                                         tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                   >  2:1403318/1‑91 (MQ=255)
                                         tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                   >  1:3389538/1‑91 (MQ=255)
                                          cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                  >  2:2074456/1‑91 (MQ=255)
                                          cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                  <  2:3255051/91‑1 (MQ=255)
                                          cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                  <  1:826289/91‑1 (MQ=255)
                                           tgaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAatat                                                 >  2:2783870/1‑91 (MQ=255)
                                              aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                              <  2:3510342/91‑1 (MQ=255)
                                                tgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattg                                            <  1:1831514/91‑1 (MQ=255)
                                                 gaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattga                                           <  1:377702/91‑1 (MQ=255)
                                                  aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                          >  2:3861849/1‑91 (MQ=255)
                                                  aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                          >  1:582449/1‑91 (MQ=255)
                                                    aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                        >  2:1144003/1‑91 (MQ=255)
                                                     gTTTCTTAGACTATAGAATAGATTAGTATCTATAAAATTTTATTAGCTACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCt                                       <  1:1108127/91‑1 (MQ=38)
                                                      tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                      >  2:1888792/1‑91 (MQ=255)
                                                      tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                      >  1:2026416/1‑91 (MQ=255)
                                                      tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                      >  1:345534/1‑91 (MQ=255)
                                                      tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                      <  1:2210/91‑1 (MQ=255)
                                                        tCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAtt                                    <  2:2704982/91‑1 (MQ=255)
                                                         cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                                   >  2:378627/1‑91 (MQ=255)
                                                         cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                                   >  2:569462/1‑91 (MQ=255)
                                                           tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                                 >  1:2810474/1‑91 (MQ=255)
                                                           tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                                 >  2:2733154/1‑91 (MQ=255)
                                                            aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                <  2:1885901/91‑1 (MQ=255)
                                                             gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                               <  2:903049/91‑1 (MQ=255)
                                                             gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                               <  1:2023861/91‑1 (MQ=255)
                                                             gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                               <  2:1276306/91‑1 (MQ=255)
                                                              aCTATAAAAGATATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                              <  1:3523625/91‑1 (MQ=255)
                                                               cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                             <  2:2388503/91‑1 (MQ=255)
                                                                  taAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                          <  2:3649684/91‑1 (MQ=255)
                                                                  aaaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                          <  1:2559084/90‑1 (MQ=255)
                                                                   aaaaGAGCTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         <  1:1932225/91‑1 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         >  1:1665143/1‑91 (MQ=255)
                                                                   aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                         >  2:2021279/1‑91 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        <  1:1723956/91‑1 (MQ=255)
                                                                    aaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                        >  2:3442111/1‑91 (MQ=255)
                                                                     aaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTCCAAAt                       >  2:3581037/1‑91 (MQ=255)
                                                                      agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      <  2:1706747/91‑1 (MQ=255)
                                                                      agagATTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                      >  1:2363663/1‑91 (MQ=25)
                                                                        agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                    <  2:2424990/91‑1 (MQ=255)
                                                                        agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                    <  2:1767973/91‑1 (MQ=255)
                                                                        agaTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                    >  1:504123/1‑91 (MQ=25)
                                                                         gaTTATTATCTATTAAATTTTATTATATACTAATCTCTTTTTGTCTACGATAATGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   <  1:3608705/91‑1 (MQ=255)
                                                                         gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   >  2:2989785/1‑91 (MQ=255)
                                                                         gaTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGtt                   >  2:3063439/1‑91 (MQ=25)
                                                                          aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                  >  2:603632/1‑91 (MQ=25)
                                                                          aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTCTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                  >  1:3276058/1‑91 (MQ=25)
                                                                           ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                 <  1:2432850/91‑1 (MQ=255)
                                                                           ttAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                 <  1:322940/91‑1 (MQ=255)
                                                                             aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               <  2:1046604/91‑1 (MQ=255)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTTTATTTTTACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  2:282512/1‑91 (MQ=21)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  2:2120132/1‑91 (MQ=37)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  1:2261645/1‑91 (MQ=37)
                                                                             aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTCGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg               >  2:466861/1‑91 (MQ=25)
                                                                              gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg              <  1:2737981/91‑1 (MQ=255)
                                                                                aTCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg            >  1:499666/1‑91 (MQ=37)
                                                                                 tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa           <  1:1412575/91‑1 (MQ=255)
                                                                                     atAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata       >  1:1970028/2‑91 (MQ=255)
                                                                                     atAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata       >  2:2549628/2‑91 (MQ=255)
                                                                                      tAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACatat      >  1:957740/1‑91 (MQ=255)
                                                                                      tAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACatat      >  1:957918/1‑91 (MQ=255)
                                                                                       aaaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATc     >  2:816599/1‑91 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTGATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  1:3165656/91‑1 (MQ=255)
                                                                                        aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc    <  1:2874073/91‑1 (MQ=255)
                                                                                         aTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCa   <  1:1525326/91‑1 (MQ=255)
                                                                                          ttttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAt  <  1:682438/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: