Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,211,670 C→A 7.9% intergenic (+663/+150) 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,6700CA 7.9% 671.7 / 30.7 193intergenic (+663/+150)USA300HOU_RS11275/yidChypothetical protein/Membrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base C (93/84);  new base A (9/7);  total (102/91)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

AGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGC  >  USA300TCH1516_ALE/2211581‑2211760
                                                                                         |                                                                                           
agaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTcaca                                                                                            <  2:205871/91‑1 (MQ=255)
 gaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAt                                                                                           <  1:3779285/91‑1 (MQ=255)
  aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                                          <  1:1573574/91‑1 (MQ=255)
  aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                                          >  2:2861185/1‑91 (MQ=255)
   aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                                         <  2:2392137/91‑1 (MQ=255)
   aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                                         <  2:2618816/91‑1 (MQ=255)
   aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                                         <  1:1990339/91‑1 (MQ=255)
    tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                                        >  1:514339/1‑91 (MQ=255)
     gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                                       >  2:1310114/1‑91 (MQ=255)
     gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                                       >  1:624288/1‑91 (MQ=255)
     gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                                       >  2:3355236/1‑91 (MQ=255)
     gACAACTCTGTTATTAACCACTTTTACCCTACTGAAAACGCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                                       >  1:1791217/1‑91 (MQ=255)
      aCAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaa                                                                                      <  1:2779262/91‑1 (MQ=255)
       cAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaaa                                                                                     >  1:3491773/1‑91 (MQ=255)
       cAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaaa                                                                                     >  1:3491724/1‑91 (MQ=255)
          ctctGTTTTTATCCACTTTTTACATACTGAAAACGCGTTAATATTATGTCATATAGAAACAGCAGTAAGATGATTTCAAATGTGaaagtca                                                                                 <  2:1000027/91‑5 (MQ=12)
          ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                                 >  2:2112277/1‑91 (MQ=255)
          ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                                 >  2:1519977/1‑91 (MQ=255)
          ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                                 >  2:531041/1‑91 (MQ=255)
            ctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGtcatc                                                                               <  2:1975395/91‑1 (MQ=21)
            ctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGtcatc                                                                               <  2:2332198/91‑1 (MQ=21)
             tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCt                                                                              <  2:1384482/91‑1 (MQ=255)
                   taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                        <  2:2782057/91‑1 (MQ=255)
                   taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                        <  1:3910473/91‑1 (MQ=255)
                   taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                        >  1:3886122/1‑91 (MQ=255)
                   taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTActgc                                                                        <  1:3071617/91‑1 (MQ=17)
                     aCCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgctg                                                                      <  2:3758978/91‑1 (MQ=255)
                      ccACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGt                                                                     >  1:1720526/1‑91 (MQ=255)
                       cACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGtt                                                                    <  2:683022/91‑1 (MQ=255)
                        aCTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttt                                                                   <  1:337939/91‑1 (MQ=255)
                         cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttta                                                                  <  2:501737/91‑2 (MQ=255)
                         cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                                  <  1:316110/91‑1 (MQ=255)
                         cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                                  <  2:1189698/91‑1 (MQ=255)
                         cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                                  <  1:985375/91‑1 (MQ=255)
                          tttttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCt                                                                 >  1:567655/1‑91 (MQ=255)
                             ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                              >  2:1680168/1‑91 (MQ=255)
                               aCATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTa                                                            <  1:2166019/91‑1 (MQ=255)
                                cATACTGAAAGCACGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTAt                                                           <  2:1237509/91‑1 (MQ=255)
                                cATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTAt                                                           >  1:1690566/1‑91 (MQ=255)
                                 aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                                          <  2:2407825/91‑1 (MQ=255)
                                 aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                                          <  2:708263/91‑1 (MQ=255)
                                  tACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATCTGAAATCATCTTACTGCTGTTTCTATTTATca                                                         >  1:3364436/1‑91 (MQ=255)
                                   aCTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcac                                                        >  2:3029883/1‑91 (MQ=255)
                                    cTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcaca                                                       <  2:2563653/91‑1 (MQ=255)
                                         aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGGCTTTCACATTTGAAATCATCTTACTGCTGTTTCTATATATCACATATTg                                                  <  2:3462046/91‑1 (MQ=255)
                                         aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTg                                                  <  2:2730562/91‑1 (MQ=255)
                                          aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                                 <  1:370999/91‑1 (MQ=255)
                                           cTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGta                                                <  1:1604078/91‑1 (MQ=255)
                                            tCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtat                                               >  1:2888759/1‑91 (MQ=255)
                                              gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                             >  2:1816971/1‑91 (MQ=255)
                                              gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                             >  2:146981/1‑91 (MQ=255)
                                              gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                             >  2:3490508/1‑91 (MQ=255)
                                               ttAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAt                                            >  1:1866684/1‑91 (MQ=255)
                                                 aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                          <  2:864767/91‑1 (MQ=255)
                                                 aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                          >  2:1167514/1‑91 (MQ=255)
                                                 aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                          <  1:3259578/91‑1 (MQ=255)
                                                   tattatTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGa                                        <  2:2533879/91‑1 (MQ=255)
                                                    attattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATCTTGTATAATGTGac                                       <  2:1857544/91‑1 (MQ=255)
                                                     ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                                      >  1:1053891/1‑91 (MQ=255)
                                                     ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                                      >  1:74612/1‑91 (MQ=255)
                                                       attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                                    >  1:3580974/1‑91 (MQ=255)
                                                        tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                                   <  2:3922191/91‑1 (MQ=255)
                                                        tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                                   <  2:2354220/91‑1 (MQ=255)
                                                         ttCAGATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                                  <  2:2689392/91‑1 (MQ=255)
                                                         ttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                                  <  1:1597155/91‑1 (MQ=255)
                                                           cAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGt                                <  2:1754574/91‑1 (MQ=255)
                                                           cAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGt                                <  2:3211635/91‑1 (MQ=255)
                                                             aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                              >  2:802669/1‑91 (MQ=255)
                                                             aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                              <  2:60826/91‑1 (MQ=255)
                                                             aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                              >  2:1179224/1‑91 (MQ=255)
                                                               tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                            >  1:2465573/1‑91 (MQ=255)
                                                                agaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGc                           >  1:2339208/3‑91 (MQ=11)
                                                                aaaaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGc                           >  1:2245449/1‑91 (MQ=255)
                                                                 gaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt                          >  1:1158714/2‑91 (MQ=11)
                                                                 aaaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt                          >  1:73374/1‑91 (MQ=255)
                                                                  aaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTa                         >  1:2633553/1‑91 (MQ=255)
                                                                    aCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTAtt                       >  2:2914965/1‑91 (MQ=14)
                                                                      agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa                     >  1:2074786/1‑91 (MQ=255)
                                                                      agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa                     >  1:1849712/1‑91 (MQ=255)
                                                                      agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa                     >  1:2558705/1‑91 (MQ=255)
                                                                       gcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa                    >  1:2894558/1‑91 (MQ=255)
                                                                       gcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa                    <  2:1613635/91‑1 (MQ=255)
                                                                       gcagTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa                    >  2:1667687/1‑91 (MQ=17)
                                                                        cagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAg                   >  2:2149931/1‑91 (MQ=255)
                                                                        cagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAg                   >  1:1237846/1‑91 (MQ=255)
                                                                         agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc                  >  1:3006682/1‑91 (MQ=255)
                                                                         agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc                  <  1:1040281/91‑1 (MQ=255)
                                                                         agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc                  >  2:1889431/1‑91 (MQ=255)
                                                                          gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg                 <  2:329944/91‑1 (MQ=255)
                                                                          gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg                 <  1:3572804/91‑1 (MQ=255)
                                                                                tGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTAATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAaat           <  1:931672/91‑1 (MQ=255)
                                                                                atgatttCAAATGTGAAAGTCATCCTACTGCTGTCTGTATTTATCACATACTGTCTAATGTGACACTAAGTTTCGCTCTTTAAGCGcaaaa            >  2:2496526/5‑91 (MQ=12)
                                                                                 gACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAaata          >  1:1351728/1‑91 (MQ=255)
                                                                                   atttCAAATGTGAAAGTCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAaataa         >  2:1430027/2‑91 (MQ=37)
                                                                                     ttCACATTGGGAATCATCGTACTGCTGTTTCTATTTATCACATATTCTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAAtgt      <  1:1305044/91‑1 (MQ=255)
                                                                                      tCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACAGTAAGTTTCGCTATTGAAGCGAAAAATAAtgtg     <  2:3886122/91‑1 (MQ=255)
                                                                                       cacaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATACTGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgc    <  1:3428194/91‑1 (MQ=255)
                                                                                       cacaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgc    >  1:1909017/1‑91 (MQ=255)
                                                                                        acaTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcg   <  2:3791553/91‑1 (MQ=255)
                                                                                         caTTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTgcgc  >  1:1248007/1‑91 (MQ=255)
                                                                                         |                                                                                           
AGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATAATGTGCGC  >  USA300TCH1516_ALE/2211581‑2211760

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