Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,372,407 A→G 5.1% intergenic (+158/+34) 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,4070AG 5.1% 1148.3 / 30.4 328intergenic (+158/+34)USA300TCH1516_02262/pbuGhypothetical protein/Guanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base A (156/155);  new base G (4/13);  total (160/168)
Fisher's exact test for biased strand distribution p-value = 4.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 338 total aligned reads displayed.

TATTTATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAT  >  USA300TCH1516_ALE/2372317‑2372497
                                                                                          |                                                                                          
tatTTATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTa                                                                                            <  2:1232969/91‑1 (MQ=255)
tatTTATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTa                                                                                            <  1:2855414/91‑1 (MQ=255)
tatTTATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTAAAAAAAAAAAAAAAACGACTa                                                                                            <  2:2685377/91‑1 (MQ=255)
  tttATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTAcc                                                                                          <  1:3138544/91‑1 (MQ=255)
  tttATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTAAAAGAAAAAAACGACTAcc                                                                                          <  2:202135/91‑1 (MQ=255)
   ttATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCa                                                                                         <  1:1283950/91‑1 (MQ=255)
    tataTCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCaa                                                                                        <  1:3861795/91‑1 (MQ=255)
    tataTCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCaa                                                                                        >  2:1414204/1‑91 (MQ=255)
      tatCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTTTCCTTACTTCAAAGAAAAAAACGACTACCaacc                                                                                      >  2:2231080/1‑89 (MQ=255)
        tCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCAt                                                                                    >  2:2772715/1‑91 (MQ=255)
          ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTg                                                                                  <  1:2530827/91‑1 (MQ=255)
          ggATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTg                                                                                  >  1:2148645/1‑91 (MQ=255)
           gATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTgg                                                                                 <  2:1739070/91‑1 (MQ=255)
            atatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGc                                                                                <  2:306611/91‑1 (MQ=255)
             tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggta                                                                               <  1:4781156/91‑3 (MQ=39)
             tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTggta                                                                               <  1:275333/91‑3 (MQ=39)
             tatTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCa                                                                               <  2:3678377/91‑1 (MQ=255)
              atTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAg                                                                              <  1:5201019/91‑1 (MQ=255)
                 tAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCg                                                                           >  2:3972595/1‑91 (MQ=255)
                  aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                          >  1:5560807/1‑91 (MQ=255)
                  aaaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGt                                                                          >  2:351803/1‑91 (MQ=255)
                   aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTGCCAATGATTGGTAGTCGtt                                                                         <  1:295968/91‑1 (MQ=38)
                   aaGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGtt                                                                         >  1:1916933/1‑91 (MQ=255)
                    aGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttt                                                                        >  1:3296222/1‑91 (MQ=255)
                    aGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttt                                                                        <  1:342224/91‑1 (MQ=255)
                      aCATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGttttt                                                                      <  2:3303797/91‑1 (MQ=255)
                        aTTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtg                                                                    <  2:4747758/91‑1 (MQ=255)
                        aTTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtg                                                                    <  2:1049594/91‑1 (MQ=255)
                         ttCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTtgt                                                                   >  1:1168678/1‑91 (MQ=255)
                           cATAATATCGTTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAGCGTCTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                 <  2:4484397/91‑1 (MQ=255)
                           cATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTgtgt                                                                 >  2:909784/1‑91 (MQ=255)
                             taataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCt                                                               <  2:3857512/91‑1 (MQ=255)
                              aataTCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTg                                                              >  2:2917597/1‑91 (MQ=255)
                                 atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                           >  1:98299/1‑91 (MQ=255)
                                 atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                           >  2:1787588/1‑91 (MQ=255)
                                 atCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAt                                                           <  1:3335990/91‑1 (MQ=255)
                                   ctatttGTGTTATTTCGCCTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                         <  1:4733697/88‑1 (MQ=255)
                                   cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                         <  2:537982/91‑1 (MQ=255)
                                   cGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAttt                                                         >  1:5101461/1‑91 (MQ=255)
                                    gCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTa                                                        <  1:4360331/91‑1 (MQ=255)
                                     cTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTTTCTGGATTTaa                                                       >  2:5650075/1‑91 (MQ=255)
                                     cTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTaa                                                       >  1:1333440/1‑91 (MQ=255)
                                      tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAc                                                      <  1:2588763/91‑1 (MQ=255)
                                      tttATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAAc                                                      >  2:1302287/1‑91 (MQ=255)
                                        tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                    <  1:624197/91‑1 (MQ=255)
                                        tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                    <  2:4912400/91‑1 (MQ=255)
                                        tATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCa                                                    >  1:225158/1‑91 (MQ=255)
                                           gTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGa                                                 >  2:683476/1‑91 (MQ=255)
                                            ttATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTTGCAGTCGTTTTTTGTGTCTGGATTTAACCATGaa                                                >  2:3384100/1‑91 (MQ=255)
                                             tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                               <  1:4661332/91‑1 (MQ=255)
                                             tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                               >  1:948072/1‑91 (MQ=255)
                                             tATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                               >  1:2310222/1‑91 (MQ=255)
                                             tATTTCGCTTCTGATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAAc                                               <  2:3734080/91‑1 (MQ=255)
                                              aTTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACa                                              >  2:2573072/1‑91 (MQ=255)
                                                 tCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAat                                           >  2:5315228/1‑91 (MQ=255)
                                                  cGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAata                                          >  2:1479375/1‑91 (MQ=255)
                                                    cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                        >  1:857678/1‑91 (MQ=255)
                                                    cttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATg                                        >  1:4963449/1‑91 (MQ=255)
                                                     ttcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                       <  2:115839/91‑1 (MQ=255)
                                                      tcttATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGaa                                      <  2:2973199/91‑1 (MQ=255)
                                                         tATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTa                                   <  1:3577909/91‑1 (MQ=255)
                                                           ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGa                                       <  2:4174680/85‑1 (MQ=255)
                                                           ttttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACa                                 >  2:392525/1‑91 (MQ=255)
                                                            tttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAg                                <  2:3768807/91‑1 (MQ=255)
                                                             ttATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAga                               <  1:4056769/91‑1 (MQ=255)
                                                                    tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                        <  1:200919/91‑1 (MQ=255)
                                                                    tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                        <  2:4408405/91‑1 (MQ=255)
                                                                    tACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATa                        >  2:1651577/1‑91 (MQ=255)
                                                                      cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                      <  1:4883897/91‑1 (MQ=255)
                                                                      cTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGt                      <  1:4158564/91‑1 (MQ=255)
                                                                       ttCAAAGAAAAAAACGACTACCAGTCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                     <  2:1142912/91‑1 (MQ=255)
                                                                       ttCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTa                     <  1:4468722/91‑1 (MQ=255)
                                                                         cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                   <  2:1935157/91‑1 (MQ=255)
                                                                         cAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTAcc                   <  1:1215142/91‑1 (MQ=255)
                                                                          aaaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCa                  <  1:300637/91‑1 (MQ=255)
                                                                           aaGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAt                 <  2:1878053/91‑1 (MQ=255)
                                                                            aGAAACAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                >  2:5526664/1‑91 (MQ=255)
                                                                            aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                >  1:2399864/1‑91 (MQ=255)
                                                                            aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                >  1:3947781/1‑91 (MQ=255)
                                                                            aGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCAtt                >  1:1045437/1‑91 (MQ=255)
                                                                              aaaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTaa              >  1:2339356/1‑91 (MQ=255)
                                                                               aaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAat             >  2:956325/1‑91 (MQ=255)
                                                                               aaaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAat             <  1:351803/91‑1 (MQ=255)
                                                                                aaaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAata            >  1:4549751/1‑91 (MQ=255)
                                                                                 aaaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatac           >  2:5064631/1‑91 (MQ=255)
                                                                                  aaaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAataca          <  2:1530472/91‑1 (MQ=255)
                                                                                  aaaCGACTACCAATCATTGACAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGACGAACAGAATGAATAGTACCATTAataca          >  1:322548/1‑91 (MQ=255)
                                                                                   aaCGACTGCCAATGATTGGTAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         >  1:1577555/1‑91 (MQ=38)
                                                                                   aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         <  2:1333440/91‑1 (MQ=255)
                                                                                   aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         >  1:1643449/1‑91 (MQ=255)
                                                                                   aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         >  2:412632/1‑91 (MQ=255)
                                                                                   aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTCAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         <  1:5488455/91‑1 (MQ=255)
                                                                                   aaCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGAGTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacat         >  1:4698998/1‑91 (MQ=255)
                                                                                    aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata        >  1:736820/1‑91 (MQ=255)
                                                                                    aCGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAatacata        >  2:3178934/1‑91 (MQ=255)
                                                                                      gACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAAtacataca      >  2:4943962/1‑91 (MQ=255)
                                                                                        cTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACAtt    >  1:1381028/1‑91 (MQ=255)
                                                                                         tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGATCAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa   <  1:4251125/91‑1 (MQ=255)
                                                                                         tACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTa   <  1:2082683/91‑1 (MQ=255)
                                                                                          aCCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAt  >  1:697448/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATTTATATCGGATATTTAAAGACATTCATAATATCGCTTTATGTTATTTCGCTTCTTATTTTATCCTTACTTCAAAGAAAAAAACGACTACCAATCATTGGCAGTCGTTTTTTGTGTCTGGATTTAACCATGAACAAATATGAAGTACAGAATGAATAGTACCATTAATACATACATTAT  >  USA300TCH1516_ALE/2372317‑2372497

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: