Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,787 T→C 5.3% intergenic (+441/+62) 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,7870TC 5.3% 752.3 / 21.8 227intergenic (+441/+62)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (93/122);  new base C (3/9);  total (96/131)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 232 total aligned reads displayed.

GGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAG  >  USA300TCH1516_ALE/954697‑954877
                                                                                          |                                                                                          
ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                            <  1:1042262/91‑1 (MQ=255)
 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                           <  1:695138/91‑1 (MQ=255)
 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                           >  1:246989/1‑91 (MQ=255)
   gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                         <  2:3359682/91‑1 (MQ=255)
   gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                         <  1:1602072/91‑1 (MQ=255)
      acgacgAAATAAATTTTGCGAAGATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                                      <  2:3809838/91‑1 (MQ=255)
      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                                      <  2:3428295/91‑1 (MQ=255)
      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                                      <  2:1248178/91‑1 (MQ=17)
      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                                      <  1:1406997/91‑1 (MQ=17)
       cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                                                     <  2:830512/91‑1 (MQ=255)
       cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgct                                                                                     <  1:3233157/91‑1 (MQ=17)
        gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                                                    <  1:1116217/91‑1 (MQ=14)
         acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                                   <  2:612840/91‑1 (MQ=255)
         acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                                   <  2:1386820/91‑1 (MQ=255)
         acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                                   <  1:3562936/91‑1 (MQ=255)
         acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                                   <  2:3645055/91‑1 (MQ=255)
          cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                                  >  2:39414/1‑91 (MQ=255)
           gaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttt                                                                                 >  1:1567309/1‑91 (MQ=255)
             aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                               <  1:1616603/91‑1 (MQ=255)
             aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                               <  2:3356542/91‑1 (MQ=255)
              ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                              >  1:2543958/1‑91 (MQ=255)
              ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                              >  1:2541472/1‑91 (MQ=255)
                aaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtat                                                                            <  2:2208651/91‑1 (MQ=255)
                 aatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtata                                                                           >  2:2165894/1‑91 (MQ=255)
                 aatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtata                                                                           <  1:1202988/91‑1 (MQ=255)
                      tGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAAt                                                                      <  1:796148/91‑1 (MQ=255)
                       gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                     >  1:2776896/1‑91 (MQ=255)
                       gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGAAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                     <  2:641521/91‑1 (MQ=255)
                       gCGAAAATATCATTTCTGGCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                     <  2:3632779/91‑1 (MQ=255)
                         gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                                   <  2:317614/91‑1 (MQ=255)
                         gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                                   >  1:3394345/1‑91 (MQ=255)
                           aaaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCta                                                                 <  2:2064840/91‑1 (MQ=255)
                             atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                               <  1:2456668/91‑1 (MQ=255)
                              tatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTatat                                                              <  2:1602776/91‑1 (MQ=255)
                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                          >  2:2374176/1‑91 (MQ=255)
                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                          >  2:2425043/1‑91 (MQ=255)
                                    ttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATCTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGt                                                        <  1:3868577/91‑1 (MQ=255)
                                       tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                                     <  2:1007042/91‑1 (MQ=255)
                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                                   <  1:3018947/91‑1 (MQ=255)
                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                                   >  1:1969373/1‑91 (MQ=255)
                                           ccACTCCCAAAAAGACTGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                                 <  2:3727908/91‑1 (MQ=255)
                                            cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                                >  1:1168468/1‑91 (MQ=255)
                                                cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                            >  2:1790864/1‑91 (MQ=255)
                                                cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                            <  1:942276/91‑1 (MQ=255)
                                                 ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                                           >  1:1701717/1‑91 (MQ=255)
                                                 ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                                           <  2:1980982/91‑1 (MQ=255)
                                                  cAAAAAGACAGCAGTAGGATACTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                          >  2:1546322/1‑91 (MQ=255)
                                                  cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                          >  2:3101893/1‑91 (MQ=255)
                                                  cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                          >  1:2220698/1‑91 (MQ=255)
                                                  ataaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                          >  1:2937791/3‑91 (MQ=2)
                                                   aaaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                                         >  1:3039689/1‑91 (MQ=255)
                                                    aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                        <  2:1567419/91‑1 (MQ=255)
                                                     aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                                       >  1:2071612/1‑91 (MQ=255)
                                                      aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                                      <  2:2830173/91‑1 (MQ=255)
                                                      aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                                      <  2:2249663/91‑1 (MQ=255)
                                                       aGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTa                                     <  1:239192/91‑1 (MQ=255)
                                                        gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                                    <  2:1859851/91‑1 (MQ=255)
                                                         aCAGCAGTAGGCTAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                                   <  2:2964001/91‑1 (MQ=255)
                                                          cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                                  >  1:1657947/1‑91 (MQ=255)
                                                           agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                                 >  2:714914/1‑91 (MQ=255)
                                                           agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                                 >  2:2670256/1‑91 (MQ=255)
                                                           agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                                 <  2:927728/91‑1 (MQ=255)
                                                            gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                                >  2:1825992/1‑91 (MQ=255)
                                                            gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                                >  2:1917310/1‑91 (MQ=255)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               <  1:969902/91‑1 (MQ=255)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               >  1:634501/1‑91 (MQ=255)
                                                              agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                              >  1:256231/1‑91 (MQ=255)
                                                                  ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                          <  1:1123536/91‑1 (MQ=255)
                                                                  ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                          >  2:2942857/1‑91 (MQ=255)
                                                                   gATAATTTCAATTTGGAAATATGTTACTGCTGTTTTTTATAGTAATAGCGATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACtt                         <  2:2952050/91‑1 (MQ=255)
                                                                    aTAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTa                        <  2:480333/91‑1 (MQ=255)
                                                                    aTAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTa                        >  1:1585043/1‑91 (MQ=255)
                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt                      >  1:563744/1‑91 (MQ=255)
                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt                      <  1:2165894/91‑1 (MQ=255)
                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt                      <  2:256211/91‑1 (MQ=255)
                                                                       aTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATg                     <  2:1331681/91‑1 (MQ=255)
                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa                    >  2:2636501/1‑91 (MQ=255)
                                                                         ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc                   >  1:2898895/1‑91 (MQ=255)
                                                                         ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc                   <  2:3203689/91‑1 (MQ=255)
                                                                           cAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActt                 <  1:2067133/91‑1 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                >  1:1808469/1‑91 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                >  1:3076747/1‑91 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                >  1:2660817/1‑91 (MQ=255)
                                                                             aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct               >  1:3503792/1‑91 (MQ=255)
                                                                              tttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtctt              <  2:1302633/91‑1 (MQ=255)
                                                                               ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTTGTTTAACATTACTACTTAATGACTTCttt             <  2:3381746/91‑1 (MQ=255)
                                                                               ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCttt             >  1:3765888/1‑91 (MQ=255)
                                                                               ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCttt             >  1:1242340/1‑91 (MQ=255)
                                                                                tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta            <  2:789940/91‑1 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           >  2:3394144/1‑91 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           >  2:2742868/1‑91 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           <  1:3835510/91‑1 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           <  2:2028801/91‑1 (MQ=255)
                                                                                  gAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtata          >  1:1177380/1‑91 (MQ=255)
                                                                                   aaaTATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTatat         >  1:1716724/1‑91 (MQ=255)
                                                                                    aaTATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATAtt        >  1:2730709/1‑91 (MQ=255)
                                                                                     atatCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATAttt       >  1:2114572/1‑91 (MQ=255)
                                                                                      tatCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTa      >  1:1214972/1‑91 (MQ=255)
                                                                                         cTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCa   >  2:552516/1‑91 (MQ=255)
                                                                                          ttACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAg  >  1:650586/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAG  >  USA300TCH1516_ALE/954697‑954877

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