Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,211,684 T→C 5.9% intergenic (+677/+136) USA300HOU_RS11275 → / ← yidC hypothetical protein/Membrane protein insertase YidC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,211,6840TC 5.9% 712.6 / 22.5 221intergenic (+677/+136)USA300HOU_RS11275/yidChypothetical protein/Membrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base T (104/104);  new base C (11/2);  total (115/106)
Fisher's exact test for biased strand distribution p-value = 2.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

Only 100 of 230 total aligned reads displayed.

TGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGTTAAAA  >  USA300TCH1516_ALE/2211594‑2211775
                                                                                           |                                                                                           
tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCt                                                                                             >  1:2491637/1‑91 (MQ=255)
  ttattaACTACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTa                                                                                           >  2:188402/1‑91 (MQ=255)
  ttattaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTa                                                                                           >  2:1549981/1‑91 (MQ=255)
  ttattaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTa                                                                                           <  1:1836769/91‑1 (MQ=255)
   tatttACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAc                                                                                          <  1:3224853/91‑1 (MQ=255)
   tattaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAc                                                                                          <  2:3165582/91‑1 (MQ=255)
   tattaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTAc                                                                                          <  2:1628750/91‑1 (MQ=17)
    attaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAct                                                                                         >  2:922140/1‑91 (MQ=255)
    attaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAct                                                                                         <  1:3534978/91‑1 (MQ=255)
    attaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAct                                                                                         >  2:257036/1‑91 (MQ=255)
    attaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAct                                                                                         >  2:1169236/1‑91 (MQ=255)
      taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                                       <  2:1136553/91‑1 (MQ=255)
      taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                                       <  2:2934057/91‑1 (MQ=255)
      taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                                       >  1:131560/1‑91 (MQ=255)
       aaCCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgct                                                                                      >  2:1801311/1‑91 (MQ=255)
           aCTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttt                                                                                  >  1:769527/1‑91 (MQ=255)
             tttttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCt                                                                                >  2:3306525/1‑91 (MQ=255)
             tttttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCt                                                                                <  1:3631700/91‑1 (MQ=255)
                ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                                             >  1:588790/1‑91 (MQ=255)
                 tACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAttt                                                                            >  1:2303361/1‑91 (MQ=255)
                    aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                                                         <  2:2681417/91‑1 (MQ=255)
                           aaaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATAtt                                                                  <  1:3699014/91‑1 (MQ=255)
                           aaaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATAtt                                                                  <  1:3756957/91‑1 (MQ=255)
                             aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGACATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                                                >  2:3057821/1‑91 (MQ=255)
                             aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                                                >  1:2736495/1‑91 (MQ=255)
                              cTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCGCATATTGta                                                               >  1:3710461/1‑91 (MQ=255)
                               tCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtat                                                              >  1:2674967/1‑91 (MQ=255)
                                cGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtata                                                             <  2:2389856/91‑1 (MQ=255)
                                cGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtata                                                             <  2:1700117/91‑1 (MQ=255)
                                   tAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtg                                                          >  2:2342349/1‑91 (MQ=255)
                                    aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCAGTTTCTATTTATCACATATTGTATAAtgt                                                         >  2:1744300/1‑91 (MQ=255)
                                      tattatTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGa                                                       <  2:2897865/91‑1 (MQ=255)
                                       attattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGac                                                      <  2:1369528/91‑1 (MQ=255)
                                        ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                                                     <  1:3286367/91‑1 (MQ=255)
                                         tattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGacac                                                    <  1:3040398/91‑1 (MQ=255)
                                         tattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGacac                                                    <  2:585783/91‑1 (MQ=255)
                                         tattTCAAATAAAAACAGCAGTAGGATGACTCTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGacac                                                    >  2:2947812/1‑91 (MQ=255)
                                          attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                                                   >  2:1302670/1‑91 (MQ=255)
                                          attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                                                   >  1:3828715/1‑91 (MQ=255)
                                          attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                                                   <  2:2529370/91‑1 (MQ=255)
                                           tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                                                  >  1:1427222/1‑91 (MQ=255)
                                           tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                                                  >  2:2804710/1‑91 (MQ=255)
                                            gataaaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                                                 <  2:2033378/88‑1 (MQ=255)
                                              cAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGt                                               >  2:1651630/1‑91 (MQ=255)
                                                aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                                             >  1:531734/1‑91 (MQ=255)
                                                  tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                                           <  1:1405335/91‑1 (MQ=255)
                                                  tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                                           <  2:1811815/91‑1 (MQ=255)
                                                   agaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGc                                          >  1:1456074/3‑91 (MQ=14)
                                                    aaaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt                                         >  2:2861850/1‑91 (MQ=255)
                                                       aCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCttca                                      <  1:1262119/91‑4 (MQ=255)
                                                       aCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTAtt                                      >  1:1362336/1‑91 (MQ=255)
                                                        cagcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCttcaa                                     <  1:1019678/91‑5 (MQ=255)
                                                        cagcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTg                                     >  1:719758/1‑91 (MQ=255)
                                                         agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa                                    >  2:1534553/1‑91 (MQ=255)
                                                          gcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa                                   >  2:3635406/1‑91 (MQ=255)
                                                          gcagTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa                                   >  2:3834178/1‑91 (MQ=17)
                                                            agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc                                 <  1:999941/91‑1 (MQ=255)
                                                            agTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc                                 >  1:193415/1‑91 (MQ=17)
                                                             gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg                                <  1:623452/91‑1 (MQ=255)
                                                              tAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGa                               >  1:3142535/1‑91 (MQ=255)
                                                               aGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaa                              <  1:930295/91‑1 (MQ=255)
                                                               aGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaa                              <  1:1614827/91‑1 (MQ=255)
                                                                 gATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAA                            <  2:769527‑M1/91‑1 (MQ=255)
                                                                 gATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAA                            >  1:2832324‑M1/1‑91 (MQ=255)
                                                                  aTGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAA                           <  1:773236‑M1/91‑1 (MQ=255)
                                                                  aTGACTTCCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAA                           >  2:644147‑M1/1‑91 (MQ=255)
                                                                    tgatttcaaATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAAT                          >  2:1070488‑M1/4‑91 (MQ=255)
                                                                    tgatttcaaATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAAT                          >  2:1067793‑M1/4‑91 (MQ=255)
                                                                    gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA                         <  1:220150‑M1/91‑1 (MQ=255)
                                                                    gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA                         <  2:650039‑M1/91‑1 (MQ=255)
                                                                    gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA                         <  2:3239185‑M1/91‑1 (MQ=255)
                                                                    gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA                         <  1:1759980‑M1/91‑1 (MQ=255)
                                                                    gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA                         <  1:831611‑M1/91‑1 (MQ=255)
                                                                     aCTTTCACATTTGAAATCATCTTACTTCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAA                        <  1:3637879‑M1/91‑1 (MQ=255)
                                                                      cTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAAT                       <  1:923909‑M1/91‑1 (MQ=255)
                                                                      cTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAAT                       <  1:188402‑M1/91‑1 (MQ=255)
                                                                        ttCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGT                     <  2:273687‑M1/91‑1 (MQ=255)
                                                                         tCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTG                    >  2:1294429‑M1/1‑91 (MQ=255)
                                                                          cacaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGC                   >  2:3328595‑M1/1‑91 (MQ=255)
                                                                          cacaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGC                   >  2:1986581‑M1/1‑91 (MQ=255)
                                                                           acaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcg                  >  2:2481947/1‑91 (MQ=255)
                                                                           acaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcg                  >  1:2454445/1‑91 (MQ=255)
                                                                           acaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcg                  <  2:354648/91‑1 (MQ=255)
                                                                           acaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcg                  >  2:2483771/1‑91 (MQ=255)
                                                                             aTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGcc                >  1:1684599/1‑91 (MQ=255)
                                                                             aTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGcc                >  1:2054514/1‑91 (MQ=255)
                                                                               ttGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCt              >  1:2352478/1‑91 (MQ=255)
                                                                               ttGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCt              >  2:803142/1‑91 (MQ=255)
                                                                                  aaaTCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCtata           >  2:3269051/1‑91 (MQ=255)
                                                                                  aaaTCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCtata           >  2:3738722/1‑91 (MQ=255)
                                                                                   aaTCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATaa          >  1:1526392/1‑91 (MQ=255)
                                                                                    atcatcTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATaaa         >  1:3468013/1‑91 (MQ=255)
                                                                                       catcTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGt       <  1:3592113/91‑1 (MQ=255)
                                                                                       catcTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGt       >  1:577013/1‑91 (MQ=255)
                                                                                        atcTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGtt      <  2:3468013/91‑1 (MQ=255)
                                                                                        atcTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGtt      >  1:112045/1‑91 (MQ=255)
                                                                                          cTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGTTaa    >  1:223434/1‑91 (MQ=255)
                                                                                           ttACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGTTaaa   >  1:2002636/1‑91 (MQ=255)
                                                                                           ttACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGTTaaa   <  1:1169236/91‑1 (MQ=255)
                                                                                           ttACTGCTGTTTCTATTTATCACATATTGAATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATGAGTTaaaa  >  2:3373561/1‑91 (MQ=255)
                                                                                           |                                                                                           
TGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGCCCTATAAAGTTAAAA  >  USA300TCH1516_ALE/2211594‑2211775

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: