Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,413 C→G 7.3% intergenic (+164/+28) USA300TCH1516_02262 → / ← pbuG hypothetical protein/Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,372,4130CG 7.3% 1298.0 / 75.3 367intergenic (+164/+28)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (169/170);  new base G (15/12);  total (184/183)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 376 total aligned reads displayed.

TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503
                                                                                          |                                                                                          
tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                            >  2:1463468/1‑91 (MQ=255)
taacGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATc                                                                                            <  1:5118951/88‑1 (MQ=255)
 atCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCa                                                                                           <  1:714766/91‑1 (MQ=255)
 atCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTAAAAAAAAAAAAAAAACGACTACCAATCa                                                                                           <  2:1709525/91‑1 (MQ=255)
  tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                          >  1:1198679/1‑91 (MQ=255)
  tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                          >  1:2623598/1‑91 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                        >  2:4252294/1‑91 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                        >  1:2072585/1‑91 (MQ=255)
     gATATTTAAAGACATTCATAATATCGTTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTTAAATAAAAAAACGACTGCCAATGATTgg                                                                                       <  2:4717299/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTgg                                                                                       <  1:5047306/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTgg                                                                                       <  2:3799798/91‑1 (MQ=255)
      atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                      >  2:2843703/1‑91 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    >  1:4102539/1‑91 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  1:5549522/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  1:68118/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                    <  2:3459330/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGt                                                                                   <  1:2834494/91‑1 (MQ=38)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                   <  2:3229693/91‑1 (MQ=255)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                 <  1:4648909/91‑1 (MQ=38)
           tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATACTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCg                                                                                 >  1:5687603/1‑91 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                                <  2:5264229/91‑1 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                                <  2:1863197/91‑1 (MQ=255)
              aGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttt                                                                              >  1:2582450/1‑91 (MQ=255)
               gACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttt                                                                             <  2:5103854/91‑1 (MQ=255)
                 cATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttttt                                                                           <  2:2529291/91‑1 (MQ=255)
                  aTTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtg                                                                          <  1:4209535/91‑1 (MQ=255)
                   ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         >  2:2650665/1‑91 (MQ=255)
                    tCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgtg                                                                        <  2:1168259/91‑1 (MQ=255)
                    tCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgtg                                                                        <  1:3781621/91‑1 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  2:4541958/1‑91 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:100290/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:1711920/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                     <  2:4433508/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                     >  2:4083327/1‑91 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    <  2:3336216/91‑1 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    >  1:1891147/1‑91 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    >  1:3400073/1‑91 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  2:5400064/1‑91 (MQ=255)
                             cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                               <  2:2785923/91‑1 (MQ=255)
                             cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                               <  2:2921548/91‑1 (MQ=255)
                               cTTTATGTTATTTCGCTTCTTAATTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTaa                                                             >  2:4463836/1‑91 (MQ=255)
                                tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAc                                                            >  2:562484/1‑91 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                          <  2:3547261/91‑1 (MQ=255)
                                   aTGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCAt                                                         >  1:1300349/1‑91 (MQ=255)
                                    tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATg                                                        >  1:1630765/1‑91 (MQ=255)
                                     gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                       <  1:1390067/91‑1 (MQ=255)
                                     gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                       <  1:4750103/91‑1 (MQ=255)
                                      ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                      <  1:5151062/91‑1 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    <  1:5725696/91‑1 (MQ=255)
                                          ttCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACaaa                                                  >  1:369143/1‑91 (MQ=255)
                                           tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                                 <  2:2977968/91‑1 (MQ=255)
                                            cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                                >  2:5188628/1‑91 (MQ=255)
                                            cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                                <  2:3903236/91‑1 (MQ=255)
                                            cGCTTCTTATTGTATCCTTGCTTCAAAGAGAAAAACGACGACCAATCATTGGCAGTCGTGTTTTGTGTCTGGATTTAACCATGAACAAata                                                <  2:4925723/91‑1 (MQ=255)
                                                tcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            >  2:3685422/1‑59 (MQ=255)
                                                 cttATTTTATCCTTATTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                           <  2:4522131/91‑1 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                           >  2:1951275/1‑91 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                           <  1:2161215/91‑1 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                          <  1:407408/91‑1 (MQ=255)
                                                  ttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGt                                          <  2:1248167/91‑1 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         <  1:4954396/91‑1 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         <  1:3227418/91‑1 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         >  2:5726738/1‑91 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                         <  1:5288961/91‑1 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        >  1:1360537/1‑91 (MQ=255)
                                                    aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                        >  2:209790/1‑91 (MQ=255)
                                                      tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                      <  2:863573/91‑1 (MQ=255)
                                                      tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                      <  1:2650665/91‑1 (MQ=255)
                                                       ttATCCTTACTTCAAAGAAAAAAACGACTACCAGTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                     >  1:5835474/1‑91 (MQ=255)
                                                       ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                     >  1:3163312/1‑91 (MQ=255)
                                                        tATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaa                                    <  2:2262334/91‑1 (MQ=255)
                                                          tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                  <  1:3001479/91‑1 (MQ=255)
                                                          tCCTTACGTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGATTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                  <  2:4816627/91‑1 (MQ=255)
                                                           ccTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatga                                 >  2:4305668/1‑91 (MQ=255)
                                                            cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                <  1:795113/91‑1 (MQ=255)
                                                              tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                              >  2:4638374/1‑91 (MQ=255)
                                                              tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                              <  1:4305668/91‑1 (MQ=255)
                                                               aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGACTAg                             >  2:683487/1‑91 (MQ=255)
                                                               aCTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAg                             >  1:3511080/1‑91 (MQ=255)
                                                                cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                            >  1:5411169/1‑91 (MQ=255)
                                                                 ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                           <  1:4541958/91‑1 (MQ=255)
                                                                 ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                           <  1:1124481/91‑1 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      <  2:581429/91‑1 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  1:4393481/1‑91 (MQ=255)
                                                                          aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata                  >  1:913969/1‑91 (MQ=255)
                                                                           aaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatac                 <  1:3863020/91‑1 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                <  1:1056786/91‑1 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                <  2:4159142/91‑1 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                >  2:759899/1‑91 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               <  1:2425569/91‑1 (MQ=255)
                                                                              aCGACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  1:59202/1‑91 (MQ=38)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  2:2397243/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  1:1847329/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  2:401941/1‑91 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             <  2:2564988/91‑1 (MQ=255)
                                                                                gACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca            >  1:3904799/1‑91 (MQ=255)
                                                                                 aCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat           >  1:157803/1‑91 (MQ=255)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  1:5768905/4‑91 (MQ=39)
                                                                                      cAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg      >  2:1309110/1‑91 (MQ=255)
                                                                                          cATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATg  <  1:2234152/91‑1 (MQ=255)
                                                                                          |                                                                                          
TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503

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: