Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,419 C→T 19.8% intergenic (+170/+22) 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,4190CT19.8% 558.9 / 118.1 252intergenic (+170/+22)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (99/103);  new base T (22/28);  total (121/131)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 261 total aligned reads displayed.

ATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTC  >  USA300TCH1516_ALE/2372329‑2372509
                                                                                          |                                                                                          
atatTTAAAGACATTCATAATATCGCTTTATGTTCTTTCGCTTCTTATTTTCTCCTTACTTCAAAGAAAAAAACGCCTACCAATCATTGGc                                                                                            >  2:1426498/1‑91 (MQ=255)
atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                            <  1:124344/91‑1 (MQ=255)
atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                            <  1:2774719/91‑1 (MQ=255)
atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                            <  1:3820410/91‑1 (MQ=255)
  atTTAAAGGCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                          <  2:279526/91‑1 (MQ=255)
   tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGt                                                                                         <  1:2092759/91‑1 (MQ=38)
    ttAATGACAAACATAATATCGCTTTATGTTATTTCGCTTCGTATATTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                        <  2:1279032/91‑1 (MQ=255)
    ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGTCAATGATTGGTAGTc                                                                                        <  2:1771364/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                       <  1:496021/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                       <  1:2235118/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                       <  1:2038917/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                       <  1:3315182/91‑1 (MQ=38)
     tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCg                                                                                       <  2:2475670/91‑1 (MQ=38)
      aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAGGAAAAAAACGACTGCCAATCATTGGCAGTCGt                                                                                      <  2:2234557/91‑1 (MQ=255)
      aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGt                                                                                      <  1:2925174/91‑1 (MQ=38)
      aaaGACATTCATAATATCGCTTTATGTTATTCCGCTTCTTATTTGATCCTTACGTCACAGAAAAAAACGACTGCCAATGATTGGTAGTCGt                                                                                      <  1:3160114/91‑1 (MQ=25)
       aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtt                                                                                     >  1:1120982/1‑91 (MQ=255)
          aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                                  <  1:77295/91‑1 (MQ=255)
           cATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtttttt                                                                                 >  1:410414/1‑91 (MQ=255)
             ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                               <  1:2176276/91‑1 (MQ=255)
             ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                               >  1:1274173/1‑91 (MQ=255)
              tCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgtg                                                                              >  1:2274071/1‑91 (MQ=255)
               cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                             >  2:277571/1‑91 (MQ=255)
                ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                            >  2:1991198/1‑91 (MQ=255)
                 taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                           >  2:1619292/1‑91 (MQ=255)
                 taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                           <  2:670581/91‑1 (MQ=255)
                 taagaTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                           <  1:3058539/91‑1 (MQ=255)
                  aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                          >  1:914689/1‑91 (MQ=255)
                    tatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGa                                                                        <  2:721930/91‑1 (MQ=255)
                    tatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGa                                                                        <  1:3356404/91‑1 (MQ=255)
                     atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                       <  1:2344217/91‑1 (MQ=255)
                     atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGAATACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                       <  1:1527573/91‑1 (MQ=255)
                       cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTCTTTTGTGTCTGGAttt                                                                     >  2:3790456/1‑91 (MQ=255)
                          tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAc                                                                  <  1:358993/91‑1 (MQ=255)
                            tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                                <  1:1785879/91‑1 (MQ=255)
                             aTGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCAt                                                               >  1:3896266/1‑91 (MQ=255)
                              tGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATg                                                              >  2:2342871/1‑91 (MQ=255)
                               gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                             >  2:430284/1‑91 (MQ=255)
                                 tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                           <  1:3645961/91‑1 (MQ=255)
                                 tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                           <  1:2008019/91‑1 (MQ=255)
                                  aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                          >  1:1036596/1‑91 (MQ=255)
                                   tttCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACaa                                                         >  2:8654/1‑91 (MQ=255)
                                      cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                                      >  2:2072161/1‑91 (MQ=255)
                                        cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                                    <  2:1939030/91‑1 (MQ=255)
                                        cttcttATTTTATCCATACTTTAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTCTTTTGTGTCTGGATATAACCATCACGAAATATg                                                    >  2:1359359/1‑91 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               <  2:2952730/91‑1 (MQ=255)
                                             tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                               >  1:772210/1‑91 (MQ=255)
                                              aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                              >  1:3853066/1‑91 (MQ=255)
                                              aTTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAc                                              <  1:2517653/91‑1 (MQ=255)
                                                tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                            >  2:3220928/1‑91 (MQ=255)
                                                 ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTTTGTCTGGATTTCACCATGAACAAATATGAAGTACAga                                           >  2:566810/1‑91 (MQ=255)
                                                   aTCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaat                                         >  1:576879/1‑91 (MQ=255)
                                                    tCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatg                                        >  1:2295809/1‑91 (MQ=255)
                                                     ccTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAATCGTTTTTTGCGTCTGGATGTAACAATGAACAAATCTGAAGTACAgaatga                                       >  2:3422315/1‑91 (MQ=255)
                                                          cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                                  >  1:466450/1‑91 (MQ=255)
                                                          cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                                  <  2:3689189/91‑1 (MQ=255)
                                                           ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                                 >  1:2850862/1‑91 (MQ=255)
                                                            tCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTAGAGAATGAATAGTAc                                <  2:3738017/91‑1 (MQ=255)
                                                             cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                               >  1:866215/1‑91 (MQ=255)
                                                              aaaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                              <  1:3811819/91‑1 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACACATATGAAGTACAGAATGAATAGTACCac                             >  2:1040257/1‑90 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                             >  2:2578450/1‑91 (MQ=255)
                                                               aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                             <  1:1355021/91‑1 (MQ=255)
                                                                aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                            >  1:2155986/1‑91 (MQ=255)
                                                                aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                            <  2:1205531/91‑1 (MQ=255)
                                                                aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                            >  2:1920429/1‑91 (MQ=255)
                                                                    aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata                        <  1:1853638/91‑1 (MQ=255)
                                                                      aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                      >  1:3276268/1‑91 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  2:961944/1‑91 (MQ=255)
                                                                       aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat                     >  2:3632457/1‑91 (MQ=255)
                                                                        aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata                    >  1:1177595/1‑91 (MQ=255)
                                                                        aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata                    <  2:536592/91‑1 (MQ=255)
                                                                        aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata                    <  1:2679098/91‑1 (MQ=255)
                                                                        aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata                    >  2:2120116/1‑91 (MQ=255)
                                                                         cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac                   <  1:3443171/91‑1 (MQ=255)
                                                                         cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac                   <  1:3201505/91‑1 (MQ=255)
                                                                          gACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca                  >  2:3827283/1‑91 (MQ=38)
                                                                           aCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat                 >  2:298036/1‑91 (MQ=255)
                                                                           aCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacat                 >  1:3472723/1‑91 (MQ=255)
                                                                            ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt                >  1:584375/4‑91 (MQ=39)
                                                                            ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt                >  1:521660/4‑91 (MQ=39)
                                                                            ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt                >  1:1657976/4‑91 (MQ=39)
                                                                            ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt                >  2:2697119/4‑91 (MQ=39)
                                                                            cTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGGATAGTACCATTAATACATACAtt                <  2:1879895/91‑1 (MQ=255)
                                                                            cNACCAATCATTGGCAGTCGTTTTTTGTGTCTGGNTTTAACCATGAACANATATGAAGTACAGAATGAATAGTACCATTAATANATACAtt                >  2:2361714/1‑91 (MQ=255)
                                                                             tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa               >  2:3277533/1‑91 (MQ=255)
                                                                              gccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGCATAGTACCATTAATACATACATTAt              >  2:2588186/2‑91 (MQ=39)
                                                                              aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt              <  1:1450956/91‑1 (MQ=255)
                                                                              aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt              <  1:2579732/91‑1 (MQ=255)
                                                                                cAATCATTGGCAGTCGTTTTTTGTGCCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCg            <  1:3095156/91‑1 (MQ=255)
                                                                                 aatgaTTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg           >  1:946672/5‑91 (MQ=255)
                                                                                 aaTGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGCATGAATAGTACCATTAATACATACATTATCgg           >  2:3876120/1‑91 (MQ=255)
                                                                                 aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAACTATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg           >  2:89442/1‑91 (MQ=255)
                                                                                   tgaTTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt         >  2:1154917/3‑91 (MQ=255)
                                                                                   tCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt         <  2:3959100/91‑1 (MQ=255)
                                                                                   tCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGAt         <  1:65029/91‑1 (MQ=255)
                                                                                    cATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGGATAGTACCATTAATACATACATTATCGGATg        <  2:3017942/91‑1 (MQ=255)
                                                                                      ttGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTa      >  1:1908507/1‑91 (MQ=255)
                                                                                         gCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActt   <  1:2405112/91‑1 (MQ=255)
                                                                                          cAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTActtc  >  1:3664384/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATGTACTTC  >  USA300TCH1516_ALE/2372329‑2372509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: