Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,468,887 A→T 5.6% intergenic (‑76/‑132) lipY ← / → USA300HOU_RS12715 Triacylglycerol lipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,468,8870AT 5.6% 203.0 / 2.5 66intergenic (‑76/‑132)lipY/USA300HOU_RS12715Triacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/27);  new base T (2/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ATAATCGAGCGATTGTGAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAGGTACCTTGCATGAC  >  USA300TCH1516_ALE/2468747‑2469023
                                                                                                                                            |                                                                                                                                        
aTAATCGAGCGATTGTGAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGtata                                                                                                                                          >  1:443045/1‑141 (MQ=255)
            ttGTGAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATTTagaaga                                                                                                                              >  1:240331/1‑141 (MQ=255)
            ttGTGAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTagaaga                                                                                                                              <  1:373706/141‑1 (MQ=255)
              gtgAAGGAAATATTTATTGACAACGCTATTAATAACACGGGTTCTAATAAATGGCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTAGCCCCTACATTCTTTGCCAAACAAAGTATATAATGTTGGAGatt                                                                                                                            >  2:194398/1‑140 (MQ=255)
                gAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATTATGTagacattatt                                                                                                                          >  2:493751/1‑134 (MQ=255)
                      aaTATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCttt                                                                                                                    <  1:379644/141‑1 (MQ=255)
                         atttattGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTc                                                                                                                 >  2:231979/1‑141 (MQ=255)
                            tattGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACt                                                                                                              <  2:146326/141‑1 (MQ=255)
                            tattGAAAACGCCATCCATAACACGGTTCCTAATCAATGTCCTAACCTACCTTTATATACTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACt                                                                                                              <  1:228840/141‑1 (MQ=255)
                               tGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCTTGCATTTTTTTCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTc                                                                                                           >  1:707057/1‑141 (MQ=255)
                                     cGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTCCATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTAACTTTCAAAttt                                                                                                     >  2:420761/1‑141 (MQ=255)
                                        cactCATACCACGGTGTCTAATCCATGTCTAAACCAACCTTAATAAATGTTATTTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAg                                                                                                  <  1:119890/138‑1 (MQ=255)
                                          ttCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGAc                                                                                                >  1:748195/1‑141 (MQ=255)
                                              tAACACGGTTTCTAATCAATGACTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACa                                                                                            >  1:473536/1‑141 (MQ=255)
                                                acacGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAAt                                                                                          >  2:247504/1‑141 (MQ=255)
                                                acacGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAAt                                                                                          >  1:281087/1‑141 (MQ=255)
                                                  acGGTTTCTAATCTATTTCATAACCTACCTTTATATATTTGATGTATACAATGGTTGTCTCTACCCTACATTCTTGGCCAAAAAAATTATATAATGTAGAAGATATTTTCTTTTTCACCTTCATATTTAAGACTACAATTg                                                                                        <  2:36943/141‑1 (MQ=255)
                                                      tttCTAATCAATGTCTCAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACa                                                                                    <  1:11009/141‑1 (MQ=255)
                                                        tCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTAAATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGt                                                                                  >  2:97911/1‑141 (MQ=255)
                                                         cTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTg                                                                                 <  1:505691/141‑1 (MQ=255)
                                                          tAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGa                                                                                <  2:620047/141‑1 (MQ=255)
                                                          tAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGa                                                                                >  2:362005/1‑141 (MQ=255)
                                                          tAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGACAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGa                                                                                <  2:465567/141‑1 (MQ=255)
                                                                 tGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTtca                                                                         <  1:820110/141‑1 (MQ=255)
                                                                  gTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTtcat                                                                        >  1:174925/1‑141 (MQ=255)
                                                                   tCTTAACCTACCTTTCTATATTTTATCTATCCAACGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGCAGATATTTTCTTTCTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTtcatc                                                                       <  1:551859/141‑1 (MQ=255)
                                                                   tCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATTTAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTCAGACTACAATTGAACAGTGCTTTTtcatt                                                                       >  2:824334/1‑140 (MQ=255)
                                                                   tCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATTTAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTtcatc                                                                       >  2:156143/1‑141 (MQ=255)
                                                                     ttAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTATAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTCAACAGTGATTTTTcttcat                                                                     >  1:355791/1‑141 (MQ=255)
                                                                     ttAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTtcatcct                                                                     >  2:183361/1‑139 (MQ=255)
                                                                      tAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCAtt                                                                    <  1:578704/141‑1 (MQ=255)
                                                                        acctacctTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATCTATACTTTTTCCCTTTCAAATTTAAGACTCCAATTGAACAGTGATTTTTCATCATtat                                                                  >  1:803590/1‑141 (MQ=255)
                                                                        acctacctTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATtat                                                                  >  1:219065/1‑141 (MQ=255)
                                                                           tacctTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATaaa                                                               >  1:748411/1‑140 (MQ=255)
                                                                                  ctatatTTTATGTATCCAATGATTGTCTATCACCTACATTCTTTGCCAAAAAAAGTATATAACGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTGTTCATCATTATAACAGACAAc                                                        >  2:367628/2‑141 (MQ=255)
                                                                                    atatTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTAAAAAAATTAGAGGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTa                                                      <  1:897931/141‑1 (MQ=255)
                                                                                      gttttATGTACACAATGTTTGCCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGa                                                    <  2:428489/140‑1 (MQ=255)
                                                                                       ttttATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGAc                                                   >  2:688575/1‑141 (MQ=255)
                                                                                        tttATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACa                                                  <  1:688575/141‑1 (MQ=255)
                                                                                          ctgggaTCCAATGTGTTTTTATCGCTTACATGTTTTGCCAAATAAAGTATATAATGTAGAAGATATTTTTTTTTTCTCTTTCAAATTTAAGACTATAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACata                                                <  2:707057/136‑1 (MQ=255)
                                                                                           aTGTATCCAATGATTGTCTATCCCCTATCTTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTTCAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACatat                                               <  2:590646/141‑1 (MQ=255)
                                                                                            tGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATAtt                                              <  1:299511/141‑1 (MQ=255)
                                                                                               aTCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGAt                                           >  1:291678/1‑141 (MQ=255)
                                                                                               aTCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTCTCAAATTTAAGACTACAATTGAACAGTGATTTTTCTTCATTTTAACAGACAACTAGACTTATTgtt                                           >  2:678572/1‑139 (MQ=255)
                                                                                                  cAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAg                                        <  1:817465/141‑1 (MQ=255)
                                                                                                   aaTGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGt                                       >  2:614373/1‑141 (MQ=255)
                                                                                                     tGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTaa                                     <  1:335737/141‑1 (MQ=255)
                                                                                                      gATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTaaa                                    >  1:458398/1‑141 (MQ=255)
                                                                                                      gATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTCCAAATTGATGGCTAAAAAGGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTaaa                                    >  1:172071/1‑141 (MQ=255)
                                                                                                      gATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGAAATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATAATTATAACTAACAACTAGACATTTTGGTAAGTgaa                                    >  1:750301/1‑141 (MQ=255)
                                                                                                          gTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGCCAACTAGACATATTGATAAGTaaagaaa                                >  2:375512/1‑141 (MQ=255)
                                                                                                          gTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTaaagaaa                                >  1:331315/1‑141 (MQ=255)
                                                                                                                ccccTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACt                          >  2:850705/1‑141 (MQ=255)
                                                                                                                 cccTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACtt                         >  2:425670/1‑141 (MQ=255)
                                                                                                                 cccTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACtt                         <  1:361236/141‑1 (MQ=255)
                                                                                                                 cccTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACtt                         >  1:696555/1‑141 (MQ=255)
                                                                                                                 cccTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCAGTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACtt                         <  1:176350/141‑1 (MQ=255)
                                                                                                                 cccTACATTCTTTGCCAAAAAAAGTATATAATGAAGAAGAACTTTCCTTGTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACtt                         <  1:852698/141‑1 (MQ=255)
                                                                                                                   cTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACACTTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACGTATTGATAAGTAAAGAAAAGAACTTta                       >  1:655992/1‑141 (MQ=255)
                                                                                                                    tACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTtat                      >  1:338350/1‑141 (MQ=255)
                                                                                                                    tACATTCTTTCCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTtat                      >  1:660507/1‑141 (MQ=255)
                                                                                                                      aaTTCTTTGACAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACATTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATAc                    <  2:494255/140‑1 (MQ=255)
                                                                                                                          cTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCGAATTTAAGACTACAATTGAACAGTAATTTTTCATCATTATAACAGACAACTAGGCATATTGATAAGTAAAGAAAAGAACTTTATACGGAg                >  1:602309/1‑141 (MQ=255)
                                                                                                                           tttGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAgg               <  2:473536/141‑1 (MQ=255)
                                                                                                                              gCCAAAAAAAGTATTTAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAGGTAc            <  1:156143/141‑1 (MQ=255)
                                                                                                                              gCCAAAAAAAGTATTTAATGTAGAAAATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAGGTAc            <  1:824334/141‑1 (MQ=255)
                                                                                                                                        gTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAGGTACCTTGCATGAc  <  1:231979/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
ATAATCGAGCGATTGTGAAGCAAATATTTATTGACAACGCTATTCATAACACGGTTTCTAATCAATGTCTTAACCTACCTTTATATATTTTATGTATCCAATGATTGTCTATCCCCTACATTCTTTGCCAAAAAAAGTATATAATGTAGAAGATATTTTCTTTTTCACTTTCAAATTTAAGACTACAATTGAACAGTGATTTTTCATCATTATAACAGACAACTAGACATATTGATAAGTAAAGAAAAGAACTTTATACGGAGGTACCTTGCATGAC  >  USA300TCH1516_ALE/2468747‑2469023

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