Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,504,989 T→C 5.7% intergenic (+293/+111) USA300HOU_RS12890 → / ← gltT hypothetical protein/Proton/sodium‑glutamate symport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,504,9890TC 5.7% 921.1 / 15.8 298intergenic (+293/+111)USA300HOU_RS12890/gltThypothetical protein/Proton/sodium‑glutamate symport protein
Reads supporting (aligned to +/- strand):  ref base T (160/116);  new base C (9/11);  total (170/128)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

TCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTT  >  USA300TCH1516_ALE/2504899‑2505079
                                                                                          |                                                                                          
tCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGt                                                                                            <  1:6315781/91‑1 (MQ=255)
ccGCCTCCGCTTTAGCGCATCCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGt                                                                                            <  2:5605635/90‑1 (MQ=255)
 ccccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTc                                                                                           <  2:3199113/89‑1 (MQ=255)
 cGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTc                                                                                           <  1:4674535/91‑1 (MQ=255)
   ccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTccc                                                                                         >  2:1151982/1‑91 (MQ=255)
      gcgcTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaa                                                                                      <  1:919956/91‑1 (MQ=255)
      gcgcGTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaa                                                                                      <  2:3500026/91‑1 (MQ=255)
          tttAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaacaac                                                                                  <  1:5429909/91‑1 (MQ=255)
          tttAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaacaac                                                                                  <  2:6286910/91‑1 (MQ=255)
            tAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCt                                                                                >  1:5958874/1‑91 (MQ=255)
             aGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATGAAAGTCCCTAACAACCTg                                                                               >  1:6879633/1‑91 (MQ=255)
             aGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTg                                                                               >  1:1222341/1‑91 (MQ=255)
              gcgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGa                                                                              >  2:2894628/1‑91 (MQ=255)
               cgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAg                                                                             <  1:4724526/91‑1 (MQ=255)
                gcatgcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAgg                                                                            <  1:2961130/91‑1 (MQ=255)
                gcatgcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAgg                                                                            <  2:2510553/91‑1 (MQ=21)
                  atgcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTc                                                                          >  1:6668513/1‑91 (MQ=255)
                      atAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTa                                                                      <  2:1986851/91‑1 (MQ=40)
                      atAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTa                                                                      <  1:6381332/91‑1 (MQ=255)
                        aaaaGCCCCTAACAACCTGAGGTCACTGCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACg                                                                    <  1:4233959/91‑1 (MQ=21)
                         aaaGCCCCTAACATCCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGc                                                                   <  1:3655151/91‑1 (MQ=255)
                         aaaGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGc                                                                   >  1:4985585/1‑91 (MQ=255)
                          aaGCCCCTAACAACCTGAGGTCACTACGCTTTGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCt                                                                  <  1:6482893/91‑1 (MQ=255)
                             ccccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTcc                                                               >  1:2981123/1‑91 (MQ=255)
                             ccccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTcc                                                               <  2:4985585/91‑1 (MQ=255)
                               ccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCCTGCATAAAAGTCCCTAACAACCAGAGGTCACTACGCTTCCGt                                                             >  2:3865579/1‑91 (MQ=255)
                               ccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGt                                                             >  1:2920838/1‑91 (MQ=255)
                                cTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGtt                                                            >  1:841580/1‑91 (MQ=255)
                                 tAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTCCCGTTc                                                           >  2:6058685/1‑91 (MQ=255)
                                 tAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCGCCCGTTc                                                           >  2:5298704/1‑91 (MQ=255)
                                  aacaacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCg                                                          <  1:4863539/91‑1 (MQ=255)
                                   acaacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGc                                                         >  1:2096061/1‑91 (MQ=255)
                                     aacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCt                                                       >  2:1549200/1‑91 (MQ=255)
                                         tGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgcgc                                                   <  1:3308523/91‑1 (MQ=255)
                                          gAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCt                                                  <  2:2722502/91‑1 (MQ=255)
                                            ggTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCttt                                                >  1:2683342/1‑91 (MQ=255)
                                             gTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTa                                               <  2:7003469/91‑1 (MQ=255)
                                              tCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCGCTACGCTTCCGTTCGCCTGCGCATTat                                              >  2:1503728/1‑90 (MQ=255)
                                              tCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAg                                              >  1:3898622/1‑91 (MQ=255)
                                                  tACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgca                                          <  1:2536670/91‑1 (MQ=255)
                                                  tACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgca                                          >  1:1463956/1‑91 (MQ=255)
                                                  caCGCCTCGGTTCGCCCGCGCGTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgca                                          <  2:5823914/90‑1 (MQ=14)
                                                    cGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatg                                        <  2:5950174/91‑1 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAATAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                      >  2:564584/1‑91 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                      <  2:3833679/91‑1 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                      >  1:5513151/1‑91 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                      <  1:4093308/91‑1 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                      <  1:6552803/91‑1 (MQ=255)
                                                      cTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTGTAGCgcatgca                                      >  2:2185736/1‑91 (MQ=255)
                                                        tCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATa                                    >  1:5325672/1‑91 (MQ=40)
                                                         cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaa                                   >  1:3445279/1‑91 (MQ=39)
                                                         cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaa                                   >  1:1510326/1‑91 (MQ=39)
                                                         cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaa                                   >  1:4804152/1‑91 (MQ=39)
                                                          ggTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaa                                  <  1:2222851/91‑1 (MQ=39)
                                                          ggTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaa                                  >  2:4858678/1‑91 (MQ=38)
                                                          ggTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaa                                  <  1:169196/91‑1 (MQ=38)
                                                           gTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaaa                                 >  1:3183890/1‑91 (MQ=39)
                                                              cGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcc                              <  1:742975/91‑1 (MQ=255)
                                                              cGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcc                              >  1:1031445/1‑91 (MQ=255)
                                                               gCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGccc                             <  2:4455756/91‑1 (MQ=255)
                                                                ccTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcccc                            >  1:3075664/1‑91 (MQ=255)
                                                                ccTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcccc                            >  1:1829208/1‑91 (MQ=255)
                                                                ccTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcccc                            >  2:2355332/1‑91 (MQ=255)
                                                                ccTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCGTCCGTTCGTCTGCGCTTTAGCGCATGCATAAAAGcccc                            >  1:6170604/1‑91 (MQ=255)
                                                                   gcgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaa                         >  1:2420873/1‑91 (MQ=255)
                                                                   gcgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaa                         >  2:3591206/1‑91 (MQ=255)
                                                                   gcgcTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaa                         >  2:4010083/1‑91 (MQ=18)
                                                                    cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGcc                                                        <  1:2841154/59‑1 (MQ=255)
                                                                    cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaac                        >  2:4728195/1‑91 (MQ=255)
                                                                    cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaac                        >  2:6196470/1‑91 (MQ=255)
                                                                    cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaac                        >  1:4782133/1‑91 (MQ=255)
                                                                    cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaac                        >  1:3437129/1‑91 (MQ=255)
                                                                    cgcCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTATGCTTCCGTTAGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaac                        <  2:6735931/91‑1 (MQ=11)
                                                                      cTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaa                      >  2:6377867/1‑91 (MQ=255)
                                                                       tttAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaac                     <  1:5112422/91‑1 (MQ=255)
                                                                        ttAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcc                    <  2:3183890/91‑1 (MQ=255)
                                                                        ttAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcc                    <  1:1385966/91‑1 (MQ=255)
                                                                        ttAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCCTCCGTTCGCCTGCGCTTTAGCGCACGCATAAAAGCCCCTAaccacc                    >  2:4602209/1‑91 (MQ=255)
                                                                         tAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCa                   >  1:6513596/1‑91 (MQ=255)
                                                                            cgcATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTGGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAg                <  1:5456215/91‑1 (MQ=255)
                                                                            cgcATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcctttg                <  2:4084326/91‑5 (MQ=255)
                                                                            cgcATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAg                <  1:1636015/91‑1 (MQ=255)
                                                                            cgcATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAg                <  1:1925506/91‑1 (MQ=255)
                                                                            cgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTTACAAcctttg                <  2:5169618/91‑5 (MQ=11)
                                                                              catgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGt              >  2:2945914/1‑91 (MQ=255)
                                                                              catgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGt              <  1:4479060/91‑1 (MQ=255)
                                                                              catgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGt              >  1:6819265/1‑91 (MQ=255)
                                                                                tgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTg            >  2:4285513/1‑91 (MQ=255)
                                                                                tgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTg            >  1:4455218/1‑91 (MQ=255)
                                                                                 gcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGt           <  2:5587556/91‑1 (MQ=255)
                                                                                  catAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTa          <  2:2562113/91‑1 (MQ=255)
                                                                                    tAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAg        <  2:3433449/91‑1 (MQ=255)
                                                                                     aaaaGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAgg       <  1:4006849/91‑1 (MQ=255)
                                                                                     aaaaGTCCCTAACAACCGGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTTACAACCTTTGGTTGTTAgg       <  1:3523551/91‑1 (MQ=14)
                                                                                      aaaGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAggg      >  1:4196650/1‑91 (MQ=255)
                                                                                       aaGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAgggg     >  1:5721933/1‑91 (MQ=255)
                                                                                        aGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc    <  1:2665558/91‑1 (MQ=255)
                                                                                        aGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc    >  2:6592248/1‑91 (MQ=255)
                                                                                        aGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc    >  2:6816118/1‑91 (MQ=255)
                                                                                          tCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCtt  >  2:690609/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTT  >  USA300TCH1516_ALE/2504899‑2505079

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