Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,399,046 C→G 10.4% intergenic (+120/+52) 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,0460CG10.4% 733.3 / 55.5 220intergenic (+120/+52)fdhD/USA300HOU_RS12350Sulfurtransferase FdhD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (104/93);  new base G (11/12);  total (115/105)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 223 total aligned reads displayed.

GATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAA  >  USA300TCH1516_ALE/2398956‑2399135
                                                                                          |                                                                                         
gATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAACTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTAc                                                                                           >  2:492799/1‑91 (MQ=255)
 atatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACt                                                                                          <  2:835610/91‑1 (MQ=255)
 atatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACt                                                                                          <  2:1497338/91‑1 (MQ=255)
  tatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTc                                                                                         >  1:1469505/1‑91 (MQ=255)
  tatTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTc                                                                                         >  1:3489068/1‑91 (MQ=255)
  tatTCTACATGAATCAACAATTTAAACGCACGATTTTTATAAGTGTAAGTAGCGCAATACAATTTAAAAAGCAGTAAGTTGATTTCTACTc                                                                                         <  1:3094473/91‑1 (MQ=255)
   atTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCg                                                                                        >  1:277878/1‑91 (MQ=255)
     tCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATGAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaa                                                                                      <  2:1375229/91‑1 (MQ=255)
     tCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaa                                                                                      >  2:3871001/1‑91 (MQ=255)
      cTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAGCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaaa                                                                                     >  2:3840726/1‑91 (MQ=255)
      cTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGaaa                                                                                     >  1:1093719/1‑91 (MQ=255)
       tACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAc                                                                                    >  2:2722433/1‑91 (MQ=255)
        aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                   >  1:3388206/1‑91 (MQ=255)
        aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                   <  2:1162542/91‑1 (MQ=255)
        aCATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAAcc                                                                                   <  1:1114787/91‑1 (MQ=255)
         cATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCa                                                                                  >  2:1792445/1‑91 (MQ=255)
            gAATCAACAATTTAAACGCGCTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACt                                                                               <  2:2673948/91‑1 (MQ=21)
            gAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCt                                                                               >  1:3852856/1‑91 (MQ=255)
             aaTCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCtt                                                                              <  2:705755/91‑1 (MQ=255)
             aaTCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCtt                                                                              <  2:2911425/91‑1 (MQ=255)
               tCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTAc                                                                            >  1:1286520/1‑91 (MQ=255)
                 aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTg                                                                          >  2:3877586/1‑91 (MQ=255)
                 aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTg                                                                          <  1:644387/91‑1 (MQ=18)
                 aacaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTg                                                                          <  2:302944/91‑1 (MQ=18)
                  acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGc                                                                         >  1:1385991/1‑91 (MQ=255)
                  acaaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGc                                                                         <  1:3153532/91‑1 (MQ=18)
                   caaTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCt                                                                        <  1:3567723/91‑1 (MQ=17)
                     aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                      >  1:2986394/1‑91 (MQ=255)
                     aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                      <  1:2714923/91‑1 (MQ=255)
                     aTTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                      >  2:3513673/1‑91 (MQ=255)
                     aTTTAAACGCACCATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttt                                                                      >  1:2989994/1‑91 (MQ=255)
                      tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                     >  2:2653576/1‑91 (MQ=255)
                      tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                     >  2:2775459/1‑91 (MQ=255)
                      tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttt                                                                     >  2:1117045/1‑91 (MQ=255)
                      tttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCtttt                                                                     <  2:1574715/91‑1 (MQ=17)
                       ttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttttt                                                                    >  2:676172/1‑91 (MQ=255)
                       ttAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCttttt                                                                    >  1:3624296/1‑91 (MQ=255)
                        tAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCtttttt                                                                   <  1:2079446/91‑1 (MQ=255)
                             gCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCt                                                              <  1:689989/91‑1 (MQ=255)
                             gCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCt                                                              >  1:2493818/1‑91 (MQ=255)
                             gCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCt                                                              >  2:2297513/1‑91 (MQ=255)
                               aCTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTaa                                                            >  1:2016854/1‑91 (MQ=255)
                                 tATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                          <  1:105438/91‑1 (MQ=255)
                                  aTTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATAc                                                         <  1:530777/91‑1 (MQ=255)
                                  aTTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATAc                                                         <  1:339585/91‑1 (MQ=255)
                                   tttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACa                                                        <  1:1043003/91‑1 (MQ=255)
                                   tttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACa                                                        >  1:468907/1‑91 (MQ=255)
                                    ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                          >  2:3268741/1‑88 (MQ=255)
                                    ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACaa                                                       <  1:1798107/91‑1 (MQ=255)
                                    ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACaa                                                       >  2:2467882/1‑91 (MQ=255)
                                    ttttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTGCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATa                                                          <  1:3268741/88‑1 (MQ=255)
                                      ttATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAAc                                                     >  2:2346378/1‑91 (MQ=255)
                                        ataagtataagtAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCt                                                   <  1:233651/91‑1 (MQ=255)
                                          aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                                 >  2:106827/1‑91 (MQ=255)
                                          aagtataagtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTat                                                 >  1:2564254/1‑91 (MQ=255)
                                           agtataagtaGCTCAACACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTatc                                                <  2:3701072/91‑1 (MQ=255)
                                                  gtaGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATc                                         <  2:1740438/91‑1 (MQ=255)
                                                   taGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCg                                        >  2:2588067/1‑91 (MQ=255)
                                                    aGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGt                                       >  2:3400356/1‑91 (MQ=255)
                                                     gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                      >  1:1937046/1‑91 (MQ=255)
                                                     gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                      >  2:3790392/1‑91 (MQ=255)
                                                     gCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGtt                                      >  2:531277/1‑91 (MQ=255)
                                                      cTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTa                                     >  2:1837233/1‑91 (MQ=255)
                                                        cAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAAt                                   >  1:819511/1‑91 (MQ=255)
                                                         aaaaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAAtt                                  <  1:1959477/91‑1 (MQ=255)
                                                           aaCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTc                                >  2:1948708/1‑91 (MQ=255)
                                                            aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                               >  1:137794/1‑91 (MQ=255)
                                                            aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                               <  2:3619964/91‑1 (MQ=255)
                                                            aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                               >  2:546848/1‑91 (MQ=255)
                                                            aCAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTcc                               >  1:3208864/1‑91 (MQ=255)
                                                              aaTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCaa                             <  2:3206833/91‑1 (MQ=255)
                                                               aTTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAat                            >  1:2477404/1‑91 (MQ=255)
                                                                  tAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtc                         >  1:2745534/1‑91 (MQ=255)
                                                                   aaaaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtct                        <  2:3120553/91‑1 (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)
                                                                     aaaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTGTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttc                      <  1:3361546/91‑1 (MQ=255)
                                                                      aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                     <  1:2035742/91‑1 (MQ=255)
                                                                      aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                     <  1:1117045/91‑1 (MQ=255)
                                                                      aaGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttct                     >  1:2868047/1‑91 (MQ=255)
                                                                       aGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATcttctt                    <  2:897097/91‑1 (MQ=255)
                                                                        gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt                   >  1:1984430/1‑91 (MQ=18)
                                                                        gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCttt                   >  2:2212662/1‑91 (MQ=18)
                                                                        gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTTTATTATTCTAATACAAACCTATCATCAGTCGTTAATTTCCAATATCTTCttt                   >  2:3597029/1‑91 (MQ=12)
                                                                        gCAGTAAGCTGGTTTCGAGTAGAAATCAACTTACTGCTATAATATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATAtcttctat                   >  2:1975125/1‑89 (MQ=12)
                                                                            tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat               <  2:1013012/91‑1 (MQ=255)
                                                                            tAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAat               >  1:3303207/1‑91 (MQ=255)
                                                                              aGTTGATTTCTACTCGAAACCATCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAatat             <  2:3134752/91‑1 (MQ=255)
                                                                              aGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAatat             >  2:3222925/1‑91 (MQ=255)
                                                                               gTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATAtt            <  1:478703/91‑1 (MQ=255)
                                                                                ttGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATAttt           >  2:2222657/1‑91 (MQ=255)
                                                                                    tttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATa       >  2:1666781/1‑91 (MQ=255)
                                                                                     ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa      >  1:580655/1‑91 (MQ=255)
                                                                                     ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa      >  1:1578365/1‑91 (MQ=255)
                                                                                     ttCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATaa      >  1:2464074/1‑91 (MQ=255)
                                                                                      tCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAAc     <  2:3348089/91‑1 (MQ=255)
                                                                                       cTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACa    <  2:1828660/91‑1 (MQ=255)
                                                                                         aCTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACaaa  <  2:567031/91‑1 (MQ=255)
                                                                                          |                                                                                         
GATATTCTACATGAATCAACAATTTAAACGCACTATTTTTATAAGTATAAGTAGCTCAAAACAATTTAAAAAGCAGTAAGTTGATTTCTACTCGAAACCAGCTTACTGCTTTTTTATTCTAATACAAACCTATCATCAATCGTTAATTTCCAATATCTTCTTTTAATATTTGATAACAAA  >  USA300TCH1516_ALE/2398956‑2399135

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: