Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 76,969 G→A 53.3% intergenic (+509/‑107) USA300HOU_RS00140 → / → USA300HOU_RS00350 hypothetical protein/Monoacylglycerol lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE76,9690GA53.3% 15.5 / 157.6 110intergenic (+509/‑107)USA300HOU_RS00140/USA300HOU_RS00350hypothetical protein/Monoacylglycerol lipase
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (35/23);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 111 total aligned reads displayed.

AAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCT  >  USA300TCH1516_ALE/76882‑77059
                                                                                       |                                                                                          
aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGt                                                                                         >  1:2044931/1‑91 (MQ=255)
aaGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGt                                                                                         >  1:2200853/1‑91 (MQ=255)
 aGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGtt                                                                                        <  1:3444029/91‑1 (MQ=255)
 aGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGtt                                                                                        <  1:1917655/91‑1 (MQ=255)
  gTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTa                                                                                       <  1:2895994/91‑1 (MQ=255)
  gTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTa                                                                                       <  1:1663511/91‑1 (MQ=255)
   ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTaa                                                                                      <  1:1952300/91‑1 (MQ=255)
   ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTaa                                                                                      <  1:2207958/91‑1 (MQ=255)
   ttAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTaa                                                                                      <  1:2207868/91‑1 (MQ=255)
     aaGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAAtt                                                                                    >  1:2501539/1‑91 (MQ=255)
        cGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtat                                                                                 <  1:1452532/91‑1 (MQ=255)
         gAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtata                                                                                <  1:3452302/91‑1 (MQ=255)
         gAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtata                                                                                >  1:4009328/1‑91 (MQ=255)
         gAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTtata                                                                                >  1:4009921/1‑91 (MQ=255)
           acacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatata                                                                              <  1:2778633/91‑1 (MQ=255)
            cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAg                                                                             <  1:2798988/91‑1 (MQ=255)
            cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAg                                                                             >  1:1668526/1‑91 (MQ=255)
            cacTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAg                                                                             >  1:1668527/1‑91 (MQ=255)
               tGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAt                                                                          <  1:3335273/91‑1 (MQ=255)
                gACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAtt                                                                         >  1:1441583/1‑91 (MQ=255)
                gACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTAtt                                                                         <  1:1930629/91‑1 (MQ=255)
                 aCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTa                                                                        >  1:89944/1‑91 (MQ=255)
                 aCATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTa                                                                        >  1:1111851/1‑91 (MQ=255)
                  cATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTAc                                                                       <  1:3720017/91‑1 (MQ=255)
                   aTGATAAATCAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATACTATTACt                                                                      >  1:2087993/1‑91 (MQ=255)
                    tGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTg                                                                     >  1:133838/1‑91 (MQ=255)
                        aaaTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCaaa                                                                 >  1:3805196/1‑91 (MQ=255)
                          aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACa                                                               <  1:1154901/91‑1 (MQ=255)
                          aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACa                                                               >  1:2902127/1‑91 (MQ=255)
                          aTTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACa                                                               <  1:1154866/91‑1 (MQ=255)
                           ttAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAt                                                              >  1:3361567/1‑91 (MQ=255)
                           ttAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAt                                                              <  1:37249/91‑1 (MQ=255)
                            tAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAtt                                                             >  1:2021592/1‑91 (MQ=255)
                            tAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACAtt                                                             <  1:2644666/91‑1 (MQ=255)
                              gTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGa                                                           >  1:529356/1‑91 (MQ=255)
                              gTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGa                                                           >  1:2644643/1‑91 (MQ=255)
                               tGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGAt                                                          >  1:2780245/1‑91 (MQ=255)
                               tGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATCGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAt                                                          >  1:4133670/1‑91 (MQ=255)
                                ggTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGATTAATTTATATAGTATTACTGCAAACATTGAtt                                                         >  1:3675035/1‑91 (MQ=255)
                                 gTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAttt                                                        >  1:940533/1‑91 (MQ=255)
                                 gTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGAttt                                                        >  1:940515/1‑91 (MQ=255)
                                      cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                   <  1:846081/91‑1 (MQ=255)
                                      cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                   <  1:4122885/91‑1 (MQ=255)
                                      cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                   <  1:2672897/91‑1 (MQ=255)
                                      cTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTc                                                   <  1:1095157/91‑1 (MQ=255)
                                       tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                                                  >  1:1164652/1‑91 (MQ=255)
                                       tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                                                  <  1:2731982/91‑1 (MQ=255)
                                       tataTTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                                                  <  1:2924002/91‑1 (MQ=255)
                                          atTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaaa                                               >  1:787193/1‑91 (MQ=255)
                                           tttttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaaaa                                              >  1:1646869/1‑91 (MQ=255)
                                           tttttttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCaaaaa                                              >  1:1316760/1‑91 (MQ=255)
                                             ttcttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTAATGCAAACATTGATTTGGGTCAAAAAtt                                            >  1:3513746/4‑91 (MQ=255)
                                               tttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAg                                          >  1:2193130/1‑91 (MQ=255)
                                                ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgg                                         <  1:2043678/91‑1 (MQ=255)
                                                ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgg                                         <  1:3155069/91‑1 (MQ=255)
                                                ttACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAggg                                        >  1:3832522/1‑91 (MQ=255)
                                                 tACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAggg                                        <  1:2484859/91‑1 (MQ=255)
                                                  aCTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAgggg                                       >  1:3947898/1‑91 (MQ=255)
                                                    tttGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTa                                     >  1:2171161/1‑91 (MQ=255)
                                                     ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                    <  1:1008334/91‑1 (MQ=255)
                                                     ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                    <  1:1255763/91‑1 (MQ=255)
                                                     ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                    <  1:1255613/91‑1 (MQ=255)
                                                     ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                    <  1:429236/91‑1 (MQ=255)
                                                     ttGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAt                                    <  1:429252/91‑1 (MQ=255)
                                                      tGCAACAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTAtt                                   >  1:3160405/1‑91 (MQ=255)
                                                       gCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa                                  <  1:583708/91‑1 (MQ=255)
                                                       gCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa                                  >  1:2851235/1‑91 (MQ=255)
                                                         aaCAGAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaa                                <  1:723900/91‑1 (MQ=255)
                                                            aGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaa                             >  1:3126848/1‑91 (MQ=255)
                                                             gAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                            >  1:2157531/1‑91 (MQ=255)
                                                             gAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                            >  1:2961403/1‑91 (MQ=255)
                                                             gAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                            <  1:2086/91‑1 (MQ=255)
                                                             gAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                            <  1:2223/91‑1 (MQ=255)
                                                             gAACCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTaaaaaaa                            <  1:3883/91‑1 (MQ=255)
                                                              aaCCTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaat                           >  1:2972409/1‑91 (MQ=255)
                                                                ccTATAAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatga                         >  1:1448564/1‑91 (MQ=255)
                                                                  tataAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaat                       >  1:4074095/1‑91 (MQ=255)
                                                                  tataAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaat                       >  1:1977010/1‑91 (MQ=255)
                                                                   ataAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaatg                      >  1:1650135/1‑91 (MQ=255)
                                                                   ataAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAaatgaatg                      <  1:3336393/91‑1 (MQ=255)
                                                                    taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAatgaatga                     >  1:288180/1‑91 (MQ=255)
                                                                    taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAatgaatga                     >  1:288279/1‑91 (MQ=255)
                                                                    taAAAAACATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAatgaatga                     <  1:3132555/91‑1 (MQ=255)
                                                                      aaaaaCATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATgaatgaat                   <  1:3624247/91‑1 (MQ=255)
                                                                       aaaaCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATc                  <  1:2530727/91‑1 (MQ=255)
                                                                       aaaaCATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATc                  <  1:2337780/91‑1 (MQ=255)
                                                                       aaaaCATGGGCATAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATc                  <  1:3758339/91‑1 (MQ=255)
                                                                           cATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATaa              >  1:3001858/1‑91 (MQ=255)
                                                                           cATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATaa              >  1:3001423/1‑91 (MQ=255)
                                                                           cATGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATaa              >  1:3064518/1‑91 (MQ=255)
                                                                            aTGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAt             >  1:1933754/1‑91 (MQ=255)
                                                                             tGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtt            <  1:56639/91‑1 (MQ=255)
                                                                             tGGGCTTAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAAtt            <  1:901729/91‑1 (MQ=255)
                                                                                   tAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCa      >  1:2329647/1‑91 (MQ=255)
                                                                                   tAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCa      <  1:1894235/91‑1 (MQ=255)
                                                                                   tAAAAGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCa      >  1:816662/1‑91 (MQ=255)
                                                                                      aGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGc   <  1:448467/91‑1 (MQ=255)
                                                                                      aGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGc   <  1:1751018/91‑1 (MQ=255)
                                                                                      aGGGTTAATTTATATAGTATTACTGCAAACATAGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGc   <  1:3532247/91‑1 (MQ=255)
                                                                                       gggTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCt  <  1:941565/91‑1 (MQ=255)
                                                                                       |                                                                                          
AAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTAAAAAAATGAATGAATCATAATTTTATCAAGCT  >  USA300TCH1516_ALE/76882‑77059

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