Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,399,039 A→G 10.9% intergenic (+113/+59) fdhD → / ← USA300HOU_RS12350 Sulfurtransferase FdhD/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,399,0390AG10.9% 642.9 / 49.7 219intergenic (+113/+59)fdhD/USA300HOU_RS12350Sulfurtransferase FdhD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (111/84);  new base G (11/13);  total (122/97)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 225 total aligned reads displayed.

AATGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGA  >  USA300TCH1516_ALE/2398949‑2399128
                                                                                          |                                                                                         
aaTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGa                                                                                           <  2:2105295/91‑1 (MQ=255)
aaTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGa                                                                                           <  1:3106455/91‑1 (MQ=255)
aaTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGa                                                                                           >  1:2442935/1‑91 (MQ=255)
aaTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGa                                                                                           <  1:1804209/91‑1 (MQ=255)
aaTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGa                                                                                           <  2:1173495/91‑1 (MQ=255)
 aTGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGAt                                                                                          >  2:3235763/1‑91 (MQ=255)
    aaCGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTc                                                                                       <  1:699937/91‑1 (MQ=255)
     aCGCTATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATctct                                                                                      >  1:3348089/1‑91 (MQ=255)
     aCGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCt                                                                                      >  2:3025606/1‑91 (MQ=255)
     aCGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCt                                                                                      >  1:567031/1‑91 (MQ=255)
      cGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTa                                                                                     >  1:662852/1‑91 (MQ=255)
       gATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTAc                                                                                    >  1:1542971/1‑91 (MQ=255)
       gATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTcgag                                                                                    <  1:1310470/91‑4 (MQ=39)
        atatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACt                                                                                   <  2:835610/91‑1 (MQ=255)
        atatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACt                                                                                   <  2:1497338/91‑1 (MQ=255)
         tatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTc                                                                                  >  1:1469319/1‑91 (MQ=255)
         tatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTc                                                                                  >  2:541490/1‑91 (MQ=255)
          atTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCg                                                                                 <  2:2664685/91‑1 (MQ=255)
          atTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCg                                                                                 >  1:3619964/1‑91 (MQ=255)
             cTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaaa                                                                              >  1:1093719/1‑91 (MQ=255)
             cTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaaa                                                                              >  2:1797403/1‑91 (MQ=255)
              tACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAc                                                                             <  2:3741378/91‑1 (MQ=255)
              tACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAc                                                                             >  2:2722433/1‑91 (MQ=255)
               aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                            <  1:1114787/91‑1 (MQ=255)
               aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                            >  1:3388371/1‑91 (MQ=255)
                  tGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGc                                                                         >  2:2226375/1‑91 (MQ=255)
                   gAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCt                                                                        >  1:3852856/1‑91 (MQ=255)
                    aaTCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCtt                                                                       <  2:705755/91‑1 (MQ=255)
                      tCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTAc                                                                     >  1:1286520/1‑91 (MQ=255)
                       caacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACt                                                                    <  1:671846/91‑1 (MQ=255)
                       caacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACt                                                                    >  2:2894911/1‑91 (MQ=255)
                        aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTg                                                                   <  2:302944/91‑1 (MQ=18)
                         acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGc                                                                  >  1:1385991/1‑91 (MQ=255)
                         acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGc                                                                  <  1:3153532/91‑1 (MQ=18)
                          caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCt                                                                 <  2:3800888/91‑1 (MQ=255)
                          caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCt                                                                 <  2:1418417/91‑1 (MQ=255)
                          caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCt                                                                 <  2:492256/91‑1 (MQ=255)
                           aaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtt                                                                <  2:2326884/91‑1 (MQ=255)
                           aaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtt                                                                <  2:3434660/91‑1 (MQ=255)
                           aaTTTAAACGCACTATTTTTATAAGTATAAGGAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTGAACTCGAAACCAGCTTACTGCtt                                                                >  2:2955887/1‑91 (MQ=255)
                            aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                               >  2:1738425/1‑91 (MQ=255)
                            aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                               >  1:2986394/1‑91 (MQ=255)
                            aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                               <  1:2714923/91‑1 (MQ=255)
                            aTTTAAACGCACCATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                               >  1:2989994/1‑91 (MQ=255)
                             tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                              >  2:2653535/1‑91 (MQ=255)
                             tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                              >  2:2775459/1‑91 (MQ=255)
                             tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                              <  2:387483/91‑1 (MQ=255)
                             tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCtttt                                                              <  2:1574715/91‑1 (MQ=17)
                              ttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttttt                                                             >  2:676172/1‑91 (MQ=255)
                              ttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttttt                                                             <  1:2002790/91‑1 (MQ=255)
                                aaaCGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCtttttaa                                                           <  2:519740/91‑3 (MQ=14)
                                 aaCGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGATGATTTCTACTCGAAACCAGCTTACTGCTTTTTTac                                                          >  2:3842275/1‑90 (MQ=255)
                                  aCGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTAtt                                                         >  1:3640882/1‑91 (MQ=255)
                                    gCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCt                                                       >  1:2498627/1‑91 (MQ=255)
                                        tATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                   <  1:105438/91‑1 (MQ=255)
                                          tttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACa                                                 <  1:1043003/91‑1 (MQ=255)
                                           ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                   >  2:3268741/1‑88 (MQ=255)
                                           ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACaa                                                >  2:2467882/1‑91 (MQ=255)
                                             ttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAAc                                              >  2:2346378/1‑91 (MQ=255)
                                              tataagtataagTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAAcc                                             >  2:3629819/1‑91 (MQ=255)
                                               ataagtataagtAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCt                                            <  1:541490/91‑1 (MQ=255)
                                               ataagtataagtAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCt                                            <  1:233651/91‑1 (MQ=255)
                                                taagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTa                                           >  1:3368662/1‑91 (MQ=255)
                                                 aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                          >  2:1742668/1‑91 (MQ=255)
                                                 aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                          >  1:2564254/1‑91 (MQ=255)
                                                 aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                          <  2:1052394/91‑1 (MQ=255)
                                                         gtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATc                                  <  1:2784085/91‑1 (MQ=255)
                                                         gtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATc                                  <  2:1328526/91‑1 (MQ=255)
                                                          taGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCg                                 >  2:2588067/1‑91 (MQ=255)
                                                             cTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTa                              >  2:1837233/1‑91 (MQ=255)
                                                             cTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTa                              >  1:3226475/1‑91 (MQ=255)
                                                              tCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTaa                             >  1:3207424/1‑91 (MQ=255)
                                                                aaaaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAAtt                           <  1:1959477/91‑1 (MQ=255)
                                                                  gataaTTTAAAAAGCAGTAAGTTGATTGCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTc                         <  2:1149921/88‑1 (MQ=255)
                                                                  aaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTc                         >  2:1948708/1‑91 (MQ=255)
                                                                   aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                        <  2:562833/91‑1 (MQ=255)
                                                                   aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                        <  1:289008/91‑1 (MQ=255)
                                                                    cAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCa                       <  1:47339/91‑1 (MQ=255)
                                                                     aaTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCaa                      <  2:3206833/91‑1 (MQ=255)
                                                                      aTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAat                     >  1:2477404/1‑91 (MQ=255)
                                                                        taaaaaaGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAatat                   >  2:2511058/3‑91 (MQ=14)
                                                                         tAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtc                  >  2:1102194/1‑91 (MQ=255)
                                                                         tAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtc                  >  1:2745534/1‑91 (MQ=255)
                                                                          aaaaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtct                 <  2:3120553/91‑1 (MQ=255)
                                                                           aaaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtctt                <  1:1277561/91‑1 (MQ=255)
                                                                              aGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATcttctt             <  2:897097/91‑1 (MQ=255)
                                                                              aGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATcttctt             <  1:1980664/91‑1 (MQ=255)
                                                                              aGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATcttctt             >  1:2022054/1‑91 (MQ=255)
                                                                               gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt            >  1:1199519/1‑91 (MQ=18)
                                                                               gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt            >  2:2212662/1‑91 (MQ=18)
                                                                               gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTATTATTCTAATACAAACCTATCATCAGTCGTTAATTTCCAATATCTTCttt            >  2:3597740/1‑91 (MQ=12)
                                                                               gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTATAATATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttctat            >  2:1975125/1‑89 (MQ=12)
                                                                                cAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCtttt           <  1:2113885/91‑1 (MQ=255)
                                                                                   tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat        >  1:3303207/1‑91 (MQ=255)
                                                                                   tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat        <  2:3564217/91‑1 (MQ=255)
                                                                                   tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat        <  1:505952/91‑1 (MQ=255)
                                                                                    aaGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAata       <  1:3556271/91‑1 (MQ=255)
                                                                                     aGTTGATTTCTACTCGAAACCATCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAatat      <  2:3134752/91‑1 (MQ=255)
                                                                                     aGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAatat      >  2:3222925/1‑91 (MQ=255)
                                                                                       ttGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATAttt    >  1:991990/1‑91 (MQ=255)
                                                                                          |                                                                                         
AATGAACGATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGA  >  USA300TCH1516_ALE/2398949‑2399128

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