Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,413 C→G 6.8% 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 6.8% 1185.2 / 60.0 333intergenic (+164/+28)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base C (163/147);  new base G (11/12);  total (174/159)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 339 total aligned reads displayed.

TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503
                                                                                          |                                                                                           
   cGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAtt                                                                                          >  1:1172651/1‑91 (MQ=255)
    ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                         >  1:2218921/1‑91 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                        <  1:4141708/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                        <  1:2075184/91‑1 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTgg                                                                                        >  2:1181807/1‑91 (MQ=255)
     gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTgg                                                                                        >  1:4443406/1‑91 (MQ=255)
      atatTTAGAGACATTCATAACATCGCGTTATGTTATTTCGCTTGTGAATTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                       <  2:4039128/91‑1 (MQ=255)
       tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggta                                                                                      <  1:275333/91‑3 (MQ=39)
       tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggta                                                                                      <  1:3329236/91‑3 (MQ=39)
       tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggta                                                                                      <  1:4781156/91‑3 (MQ=39)
       tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCa                                                                                      <  2:1524146/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                     <  1:5201019/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                     <  1:2323936/91‑1 (MQ=255)
        atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                     <  2:3765090/91‑1 (MQ=255)
        atTTAAAAAAATTCATAATATCTCTTTATATTATGTAGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                                     <  2:5075981/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                    <  2:313078/91‑1 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                    >  2:200919/1‑91 (MQ=255)
         tttAAAGACATTCATAATATCGCTTTATGTTATTCCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGt                                                                                    >  2:3363857/1‑91 (MQ=255)
          ttAAAGACTTTCATAATATCGCGTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTc                                                                                   <  2:3828427/91‑1 (MQ=37)
          ttAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTc                                                                                   <  1:2443522/91‑1 (MQ=255)
            aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                                >  1:5560807/1‑91 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGttttt                                                                            <  2:4876586/91‑1 (MQ=37)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  1:4350102/91‑1 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  1:4704967/91‑1 (MQ=255)
                aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                            <  2:3303797/91‑1 (MQ=255)
                  aTTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtg                                                                          <  2:1049594/91‑1 (MQ=255)
                   ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         >  1:3520847/1‑91 (MQ=255)
                   ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                         <  1:226265/91‑1 (MQ=255)
                     cATAATATCGTTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAGCGTCTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       <  2:4484397/91‑1 (MQ=255)
                     cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                       >  1:4847830/1‑91 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  1:2371353/91‑1 (MQ=255)
                      ataataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      >  1:2109681/1‑91 (MQ=255)
                      ataataTCGCTTTAGGTTATTTCGCTTCTTATTTTATCCTTACTTAAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTc                                                                      <  2:174518/91‑1 (MQ=255)
                       taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                                     <  1:2003219/91‑1 (MQ=255)
                        aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                                    >  2:2917597/1‑91 (MQ=255)
                          tatCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGa                                                                  <  2:2658978/91‑1 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  1:450257/1‑91 (MQ=255)
                           atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                                 >  2:1787588/1‑91 (MQ=255)
                            tCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAtt                                                                >  1:2952251/1‑91 (MQ=255)
                             cacTTTGTGTTATTTCGCATCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                               <  1:4733949/89‑1 (MQ=255)
                             cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                               >  1:5101461/1‑91 (MQ=255)
                               cTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTaa                                                             >  1:1333440/1‑91 (MQ=255)
                                tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAc                                                            >  2:1302287/1‑91 (MQ=255)
                                 ttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAcc                                                           <  2:5322626/91‑1 (MQ=255)
                                 ttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAcc                                                           >  2:823875/1‑91 (MQ=255)
                                  tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                          <  1:624197/91‑1 (MQ=255)
                                     gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                       <  1:4215229/91‑1 (MQ=255)
                                      ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                      >  1:2568062/1‑91 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                     >  1:948072/1‑91 (MQ=255)
                                       tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                                     <  1:4661332/91‑1 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    <  1:906256/91‑1 (MQ=255)
                                        aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                                    <  1:983007/91‑1 (MQ=255)
                                         tttCGCTTCTTATCTTATCCTTACTTAAATGAAAAAAACGTCTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGCTTTACCGATGAACaa                                                   >  2:555759/1‑91 (MQ=255)
                                            cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                                >  2:1479375/1‑91 (MQ=255)
                                             gCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAatat                                               <  2:1906326/91‑1 (MQ=255)
                                              cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                              >  2:2457320/1‑91 (MQ=255)
                                              cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                              <  2:4932761/91‑1 (MQ=255)
                                              cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                              >  2:5471575/1‑91 (MQ=255)
                                                 cttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAg                                           <  1:3906649/91‑1 (MQ=255)
                                                   tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAGCCATGAACAAATATGAAGTa                                         <  2:643692/91‑1 (MQ=255)
                                                     ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                             <  2:4174680/85‑1 (MQ=255)
                                                     ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                       >  2:392525/1‑91 (MQ=255)
                                                       ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                     >  1:4553429/1‑91 (MQ=255)
                                                       ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                                     <  1:4056769/91‑1 (MQ=255)
                                                            cTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaa                                <  2:5033494/91‑1 (MQ=255)
                                                             ttACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAgaatgaat                               <  2:4415027/91‑1 (MQ=255)
                                                              tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                              >  1:541997/1‑91 (MQ=255)
                                                                 ttCAAAGAAAAAAACGACTACCAGTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                           <  2:1142912/91‑1 (MQ=255)
                                                                   cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                         >  1:2835333/1‑91 (MQ=255)
                                                                   cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                         <  1:1215268/91‑1 (MQ=255)
                                                                     aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                       >  2:3486713/1‑91 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      <  1:3215727/91‑1 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  1:3276579/1‑91 (MQ=255)
                                                                      aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                      >  1:1045437/1‑91 (MQ=255)
                                                                       gAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTa                     <  2:3030622/91‑1 (MQ=255)
                                                                         aaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAat                   >  2:956325/1‑91 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                <  1:1873003/91‑1 (MQ=255)
                                                                            aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca                <  1:4900719/91‑1 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               <  2:1333440/91‑1 (MQ=255)
                                                                             aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat               >  2:412632/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  2:1572758/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  1:736820/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              <  1:3364287/91‑1 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              >  1:3751818/1‑91 (MQ=255)
                                                                              aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata              <  1:3452217/91‑1 (MQ=255)
                                                                               cGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  1:1601100/1‑91 (MQ=255)
                                                                                gACTGCCAATGATTGGTAGTCGTTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacatac             >  2:3720154/1‑91 (MQ=37)
                                                                                gACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca            >  1:3130122/1‑91 (MQ=38)
                                                                                gACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca            <  1:1432563/91‑1 (MQ=255)
                                                                                  ctgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGACTAGTACCATTAATACATACAtt          >  2:2839670/4‑91 (MQ=38)
                                                                                  cTACCAATCATTGGCAGTCGTTTTTTTTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATacatacata          >  1:2966301/1‑90 (MQ=255)
                                                                                  cTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt          >  1:1381028/1‑91 (MQ=255)
                                                                                   tgccAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa         >  1:2606212/3‑91 (MQ=39)
                                                                                    aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt        >  1:697448/1‑91 (MQ=255)
                                                                                     ccAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATc       >  2:1664769/1‑91 (MQ=255)
                                                                                       aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATCGTACCATTAATACATACATTATCgg     >  2:1452210/1‑91 (MQ=255)
                                                                                       aaTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCgg     <  1:836784/91‑1 (MQ=255)
                                                                                        aTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa    <  2:3635530/91‑1 (MQ=255)
                                                                                        aTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGa    <  2:2961307/91‑1 (MQ=255)
                                                                                          cATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATg  <  2:2048438/91‑1 (MQ=255)
                                                                                          |                                                                                           
TATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTATCGGATG  >  USA300TCH1516_ALE/2372323‑2372503

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