Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,399,053 C→T 11.2% intergenic (+127/+45) 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,0530CT11.2% 653.2 / 48.1 213intergenic (+127/+45)fdhD/USA300HOU_RS12350Sulfurtransferase FdhD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (98/91);  new base T (11/13);  total (109/104)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGATTGA  >  USA300TCH1516_ALE/2398963‑2399143
                                                                                          |                                                                                          
tACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAc                                                                                            >  2:2722433/1‑91 (MQ=255)
 aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                           >  1:3388206/1‑91 (MQ=255)
 aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                           <  2:1162542/91‑1 (MQ=255)
 aCAGGAATCAACGATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAATCAATTTAATGAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                           <  1:1114982/91‑1 (MQ=255)
  cATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCa                                                                                          >  2:853406/1‑91 (MQ=255)
    tGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGc                                                                                        >  2:2226375/1‑91 (MQ=255)
     gAATCAACAATTTAAACGCGCTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACt                                                                                       <  2:2673948/91‑1 (MQ=21)
     gAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCt                                                                                       >  1:3852856/1‑91 (MQ=255)
      aaTCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCtt                                                                                      <  2:2911425/91‑1 (MQ=255)
      aaTCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCtt                                                                                      <  2:705755/91‑1 (MQ=255)
        tCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTAc                                                                                    >  1:1286520/1‑91 (MQ=255)
         caacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACt                                                                                   <  1:671846/91‑1 (MQ=255)
          aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTg                                                                                  >  2:3877586/1‑91 (MQ=255)
          aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTg                                                                                  <  2:3239285/91‑1 (MQ=18)
           acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGc                                                                                 >  1:1385912/1‑91 (MQ=255)
           acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGc                                                                                 <  1:3153532/91‑1 (MQ=18)
            caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCt                                                                                <  2:492256/91‑1 (MQ=255)
            caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCt                                                                                <  1:3567723/91‑1 (MQ=17)
             aaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtt                                                                               <  2:2326884/91‑1 (MQ=255)
              aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                              <  1:2714761/91‑1 (MQ=255)
              aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                              >  2:1738425/1‑91 (MQ=255)
              aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                              >  2:3513673/1‑91 (MQ=255)
              aTTTAAACGCACCATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                              >  1:2989994/1‑91 (MQ=255)
               tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                             >  2:2653535/1‑91 (MQ=255)
               tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                             >  2:1117045/1‑91 (MQ=255)
               tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                             >  2:2775459/1‑91 (MQ=255)
               tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCtttt                                                                             <  2:1574715/91‑1 (MQ=17)
                 tAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttttt                                                                           <  1:2079446/91‑1 (MQ=255)
                  aaaCGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTa                                                                          <  2:629680/91‑1 (MQ=255)
                    aCGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTAtt                                                                        >  1:3640882/1‑91 (MQ=255)
                      gCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCt                                                                      >  1:2493818/1‑91 (MQ=255)
                        aCTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTaa                                                                    >  1:2016854/1‑91 (MQ=255)
                          tATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                                  <  1:105438/91‑1 (MQ=255)
                            tttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACa                                                                <  1:1043003/91‑1 (MQ=255)
                             ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACaa                                                               <  1:1798107/91‑1 (MQ=255)
                             ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACaa                                                               >  2:2467882/1‑91 (MQ=255)
                               ttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAAc                                                             >  2:2346378/1‑91 (MQ=255)
                                 ataagtataagtAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCt                                                           <  1:541490/91‑1 (MQ=255)
                                  taagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTa                                                          >  1:3368662/1‑91 (MQ=255)
                                   aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                                         >  2:1742668/1‑91 (MQ=255)
                                   aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                                         >  1:2564151/1‑91 (MQ=255)
                                    agtataagtaGCTCAACACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTatc                                                        <  2:3701072/91‑1 (MQ=255)
                                      tataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTatcat                                                      <  2:3204993/91‑1 (MQ=255)
                                         aagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCaa                                                   <  2:1870714/91‑1 (MQ=255)
                                           gtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATc                                                 <  2:1328526/91‑1 (MQ=255)
                                             aGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGt                                               >  2:3400356/1‑91 (MQ=255)
                                             aGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGt                                               >  1:1437871/1‑91 (MQ=255)
                                              gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                              >  2:3790392/1‑91 (MQ=255)
                                              gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                              >  2:531277/1‑91 (MQ=255)
                                              gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                              >  1:1937046/1‑91 (MQ=255)
                                                tCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTaa                                            >  1:3207424/1‑91 (MQ=255)
                                                  aaaaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAAtt                                          <  1:1959477/91‑1 (MQ=255)
                                                   aaaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAAttt                                         <  1:2775459/91‑1 (MQ=255)
                                                    aaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTc                                        >  2:1948708/1‑91 (MQ=255)
                                                     aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                                       >  1:137794/1‑91 (MQ=255)
                                                     aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                                       <  2:562833/91‑1 (MQ=255)
                                                      cAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCa                                      <  1:47339/91‑1 (MQ=255)
                                                       aaTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCaa                                     <  2:3207424/91‑1 (MQ=255)
                                                       aaTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCaa                                     <  1:2950871/91‑1 (MQ=255)
                                                        aTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAat                                    >  1:2477404/1‑91 (MQ=255)
                                                          ttAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAatat                                  >  1:2744323/1‑91 (MQ=255)
                                                             aaaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtctt                               >  2:2499646/1‑91 (MQ=255)
                                                             aaaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtctt                               <  1:1277561/91‑1 (MQ=255)
                                                             aaaaGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtctt                               >  2:1271088/1‑91 (MQ=17)
                                                             aaaaACAGTAAGTTGATTTCTACTCGAAACCCGCTTACTGCTTTTTTATTCTAATACAAACCTATCAGCAATCGTTAATGTCCAATAtctt                               >  2:2968454/1‑91 (MQ=255)
                                                               aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                             <  1:2035742/91‑1 (MQ=255)
                                                               aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                             >  1:2868047/1‑91 (MQ=255)
                                                               aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTCTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                             <  1:2188061/91‑1 (MQ=255)
                                                                aGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATcttctt                            <  2:897097/91‑1 (MQ=255)
                                                                 gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt                           >  2:2212662/1‑91 (MQ=18)
                                                                 gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt                           >  1:1984430/1‑91 (MQ=18)
                                                                 gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTATTATTCTAATACAAACCTATCATCAGTCGTTAATTTCCAATATCTTCttt                           >  2:3597029/1‑91 (MQ=12)
                                                                 gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTATAATATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttctat                           >  2:1975125/1‑89 (MQ=12)
                                                                   aGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTa                         >  1:1096575/1‑91 (MQ=21)
                                                                    gTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTaa                        >  2:3045897/1‑91 (MQ=255)
                                                                     tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCGTCAATCGTTAATTTCCAATATCTTCTTTTAat                       >  1:3302679/1‑91 (MQ=255)
                                                                     tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat                       <  2:1013012/91‑1 (MQ=255)
                                                                     tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat                       <  1:505952/91‑1 (MQ=255)
                                                                       aGTTGATTTCTACTCGAAACCATCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAatat                     <  2:3134752/91‑1 (MQ=255)
                                                                         ttGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATAttt                   >  1:991990/1‑91 (MQ=255)
                                                                             tttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATa               <  1:755085/91‑1 (MQ=255)
                                                                             tttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATa               >  2:1666781/1‑91 (MQ=255)
                                                                             tttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATa               >  1:3472655/1‑91 (MQ=255)
                                                                              ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa              >  1:1578365/1‑91 (MQ=255)
                                                                              ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa              >  1:580655/1‑91 (MQ=255)
                                                                              ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa              >  1:436261/1‑91 (MQ=255)
                                                                              ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa              >  1:2464074/1‑91 (MQ=255)
                                                                              ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa              <  1:3358708/91‑1 (MQ=255)
                                                                               tCTACTCGCAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATGTCCAATATCTTCTTTTAATATTTGATAAc             <  2:3380596/91‑1 (MQ=255)
                                                                                cTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACa            <  2:1828660/91‑1 (MQ=255)
                                                                                cTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACa            >  2:36432/1‑91 (MQ=255)
                                                                                    tCGAAACCAGCTTACTGCTTTTTTATTCTAACACAAACCTATCATCAATCGTTAATTTCCAATATCTACTTTTAATATTCGATAACAAAAc        >  2:2964832/1‑91 (MQ=255)
                                                                                     cGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACg       <  2:509784/91‑1 (MQ=255)
                                                                                     cGAAACCAGCTTACTGCTTTTTTACTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACg       <  2:3611953/91‑1 (MQ=255)
                                                                                      gAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGa      <  2:1082690/91‑1 (MQ=255)
                                                                                       aaaCCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGAt     >  1:3254364/1‑91 (MQ=255)
                                                                                        aaCCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAATACGAtt    <  1:3021919/91‑1 (MQ=255)
                                                                                         aCCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGATTg   >  2:3025282/1‑91 (MQ=255)
                                                                                         aCCAGCTTACTGCTTTTTTATTCTAATACAAAACTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGATTg   >  2:3806039/1‑91 (MQ=255)
                                                                                          ccAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGATTGa  >  1:1148566/1‑91 (MQ=255)
                                                                                          |                                                                                          
TACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAAACGATTGA  >  USA300TCH1516_ALE/2398963‑2399143

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: