Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,768 A→T 7.9% intergenic (+422/+81) gluD → / ← glpQ_1 NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE954,7680AT 7.9% 780.1 / 41.4 240intergenic (+422/+81)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (88/132);  new base T (8/11);  total (97/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 249 total aligned reads displayed.

ACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATG  >  USA300TCH1516_ALE/954678‑954858
                                                                                          |                                                                                          
aCAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATaa                                                                                            <  1:20035/91‑1 (MQ=255)
aCAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATaa                                                                                            <  2:1925403/91‑1 (MQ=255)
    tAGTGTGCAAGTTGGGGTGGGACGACTAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTc                                                                                        >  2:2238639/1‑91 (MQ=255)
    tAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTATGATAATTTc                                                                                        >  2:1113/1‑91 (MQ=255)
      gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                      <  2:3459714/91‑1 (MQ=255)
      gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                      >  2:1691762/1‑91 (MQ=255)
       tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                     >  2:2731884/1‑91 (MQ=255)
       tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                     <  2:2651226/91‑1 (MQ=255)
       tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                     <  2:2650916/91‑1 (MQ=255)
        gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                    <  2:3862029/91‑1 (MQ=255)
        gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                    <  1:156522/91‑1 (MQ=255)
        gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                    <  2:3505339/91‑1 (MQ=255)
        gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                    <  1:3093457/91‑1 (MQ=255)
        gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                    >  1:852116/1‑91 (MQ=255)
         tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                                   >  1:2340697/1‑91 (MQ=255)
           cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                                 >  1:3200782/1‑91 (MQ=255)
             aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                               <  2:718058/91‑1 (MQ=255)
             aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                               >  1:3824862/1‑91 (MQ=255)
             aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTACCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                               <  2:1140862/91‑1 (MQ=255)
                tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGCCAGCAGTAGGATAATTTCAATTTGGACata                                                                            >  2:3562838/1‑91 (MQ=255)
                tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                            <  1:3887254/91‑1 (MQ=255)
                 ggggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAatat                                                                           <  1:803693/91‑1 (MQ=255)
                  gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                          >  2:1345614/1‑91 (MQ=255)
                  gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                          >  1:373744/1‑91 (MQ=255)
                  gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTAAAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                          <  2:3450174/91‑1 (MQ=21)
                   ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACCGCAGTAGGATAATTTCAATTAGGAAatatgt                                                                         >  1:2216962/1‑89 (MQ=255)
                   ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                         <  1:2608751/91‑1 (MQ=255)
                   ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                         >  2:558478/1‑91 (MQ=255)
                   ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                         <  1:929512/91‑2 (MQ=17)
                    gTGGGACGACGAAATAAATTTTGCGAAGATTTCATTTCTGTCCCACTCAAGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                        <  2:3666978/91‑1 (MQ=21)
                    gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                        <  2:2208700/91‑1 (MQ=255)
                    gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                        >  1:344150/1‑91 (MQ=255)
                    gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGCCCCACTCCCAAAAAGACAGCAGTAGGCTAATTTCAATTTGGAAATATccc                                                                        >  2:652858/1‑89 (MQ=255)
                      gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                      <  2:2895199/91‑1 (MQ=255)
                      gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                      <  1:1766161/91‑1 (MQ=255)
                       ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                     <  1:2262156/91‑1 (MQ=255)
                       ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                     >  1:2566455/1‑91 (MQ=255)
                        gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                    <  2:632348/91‑1 (MQ=255)
                         acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                   <  2:3212824/91‑1 (MQ=255)
                         acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                   <  2:3354671/91‑1 (MQ=14)
                         acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                   <  2:1400834/91‑1 (MQ=17)
                           gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgctg                                                                 <  1:1924023/91‑1 (MQ=255)
                           gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                                 <  1:1851459/91‑1 (MQ=14)
                           gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATCATCCTActgctg                                                                 <  1:3029696/91‑1 (MQ=9)
                             cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                               >  1:3486307/1‑91 (MQ=255)
                             cgaAATAAATGTTGCGAAAATATCATTTCTGTCCCACTCCAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                               <  2:2852727/91‑1 (MQ=255)
                               aaataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                             >  2:1639915/1‑91 (MQ=255)
                               aaataaatTTTGCGAAAATATCATTTCTGTACAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                             <  2:3367894/91‑1 (MQ=37)
                                 ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAGCGCTGtttttt                                                           >  2:3017930/1‑91 (MQ=255)
                                 ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                           >  1:13601/1‑91 (MQ=255)
                                  taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                          >  2:412323/1‑91 (MQ=255)
                                   aaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtat                                                         <  1:334020/91‑1 (MQ=255)
                                      ttttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGt                                                      >  1:3725168/1‑91 (MQ=255)
                                                   tCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                         >  2:1844426/1‑91 (MQ=255)
                                                      tttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTg                                      <  1:430937/91‑1 (MQ=255)
                                                        tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                    <  2:1797118/91‑1 (MQ=255)
                                                        tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                    <  1:314284/91‑1 (MQ=255)
                                                         cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                   <  2:1717916/91‑1 (MQ=255)
                                                         cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                   <  1:1638215/91‑1 (MQ=255)
                                                         cTGTCCCACTCCCAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                  <  2:3250474/91‑1 (MQ=255)
                                                          tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                  >  2:910085/1‑91 (MQ=255)
                                                          tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                  <  1:2032085/91‑1 (MQ=255)
                                                          tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                  <  2:3307025/91‑1 (MQ=255)
                                                           gTCCCACTCCCTAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGGTTTTTATAGTAATAGCTATATTTTTGTAttt                                 <  2:3024889/91‑1 (MQ=255)
                                                            tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                <  1:1319529/91‑1 (MQ=255)
                                                            tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                <  1:1071995/91‑1 (MQ=255)
                                                            tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                <  2:523719/91‑1 (MQ=255)
                                                             cccAGTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                               <  2:3637860/91‑1 (MQ=255)
                                                             cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                               <  2:764551/91‑1 (MQ=255)
                                                             cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                               <  1:3539409/91‑1 (MQ=255)
                                                              ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                              >  2:3717954/1‑91 (MQ=255)
                                                              ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                              <  1:2739316/91‑1 (MQ=255)
                                                               cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                             >  1:2696223/1‑91 (MQ=255)
                                                                  acacAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                          <  2:3583014/88‑1 (MQ=255)
                                                                  acaaaaaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                          <  2:2793902/87‑1 (MQ=255)
                                                                   cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                         <  1:144379/91‑1 (MQ=255)
                                                                      taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                      >  2:3386727/2‑91 (MQ=14)
                                                                      taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTCTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                      >  1:2747602/2‑91 (MQ=0)
                                                                       aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                     <  1:2402992/91‑1 (MQ=255)
                                                                       aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                     <  1:70232/91‑1 (MQ=255)
                                                                        aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                    >  1:2806627/1‑91 (MQ=255)
                                                                           gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTCATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                 >  2:2473552/1‑91 (MQ=255)
                                                                            aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                >  1:737700/1‑91 (MQ=255)
                                                                             cagcagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa               >  2:2414620/1‑91 (MQ=17)
                                                                              agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt              >  1:2371154/1‑91 (MQ=255)
                                                                               gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt             <  1:1249491/91‑1 (MQ=255)
                                                                               gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt             >  1:634420/1‑91 (MQ=255)
                                                                                cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta            <  1:1125318/91‑1 (MQ=255)
                                                                                cagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta            >  2:1462288/1‑91 (MQ=18)
                                                                                 agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac           <  1:3803159/91‑1 (MQ=255)
                                                                                 agTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac           >  2:4583/1‑91 (MQ=18)
                                                                                  gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact          >  2:966704/3‑91 (MQ=255)
                                                                                  gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact          <  2:1552644/91‑1 (MQ=255)
                                                                                   tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta         <  2:3118012/91‑1 (MQ=255)
                                                                                   tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta         >  1:3397173/1‑91 (MQ=255)
                                                                                   tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta         >  1:1588334/1‑91 (MQ=255)
                                                                                     ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact       <  2:2065995/91‑1 (MQ=255)
                                                                                        tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa    <  2:3154846/91‑1 (MQ=255)
                                                                                         aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt   >  2:3069087/1‑91 (MQ=255)
                                                                                          aTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATg  >  1:2713195/1‑91 (MQ=255)
                                                                                          |                                                                                          
ACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATG  >  USA300TCH1516_ALE/954678‑954858

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: