Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,407 G→A 50.8% E18K (GAA→AAA)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,4070GA50.8% ‑1.6 / 84.9 65E18K (GAA→AAA) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base G (11/21);  new base A (18/15);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAA  >  USA300TCH1516_ALE/1231271‑1231545
                                                                                                                                        |                                                                                                                                           
aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGa                                                                                                                                         <  1:528025/141‑1 (MQ=255)
  gtcgtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGTTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATg                                                                                                                                       <  1:525858/141‑1 (MQ=255)
       gAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGggt                                                                                                                                  <  1:72650/141‑1 (MQ=255)
         aCAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGggtgg                                                                                                                                <  1:517935/141‑1 (MQ=255)
         aCAAGGTTCGGATAGTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGggtgg                                                                                                                                <  2:639788/141‑1 (MQ=255)
            aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAt                                                                                                                             >  1:37037/1‑141 (MQ=255)
             ggTTCGGATATTGACGCAGCTTGTGGTCAATTAAGCGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGGATGAAGATGCGGGTGGTAtt                                                                                                                            <  1:212958/141‑1 (MQ=255)
                  gggatgtGACGCAGCTTGGGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                                       <  2:60553/135‑1 (MQ=255)
                  gggagttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                                       <  2:229655/136‑1 (MQ=255)
                  ggctatTGACGCAGCTTGTGGTCAATTAAGAGGAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                                       <  2:509099/138‑1 (MQ=255)
                  ggATATTGACGCTGCTTGTGGTCGATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACTATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTta                                                                                                                       <  2:229927/141‑1 (MQ=255)
                   gATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTtat                                                                                                                      <  1:643734/141‑1 (MQ=255)
                      atTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAAt                                                                                                                   >  1:236510/1‑141 (MQ=255)
                          aCGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCaaa                                                                                                               >  1:726982/1‑141 (MQ=255)
                           cGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAAc                                                                                                              >  1:378390/1‑141 (MQ=255)
                             cAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTa                                                                                                            <  1:677789/141‑1 (MQ=255)
                              aGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTaa                                                                                                           >  1:77991/1‑141 (MQ=255)
                               ccTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAAt                                                                                                          <  1:390791/140‑1 (MQ=255)
                                  tgtgGTCAATTAAGAGCAAAGGGACGACAAGTAGACACGAGGTCAAGATAAATGCTAGCGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATCATCAAAGTAATcaa                                                                                                       >  2:188302/1‑141 (MQ=255)
                                   gtgGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaac                                                                                                      >  1:404951/1‑141 (MQ=255)
                                   ctgGTCAATGAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAAGGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaac                                                                                                      <  2:656793/140‑1 (MQ=255)
                                     ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                    >  2:706768/1‑141 (MQ=255)
                                     ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATACTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                    >  1:521562/1‑141 (MQ=255)
                                               gagCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTc                                                                                          >  1:147868/1‑141 (MQ=255)
                                                agCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCt                                                                                         >  1:164699/1‑141 (MQ=255)
                                                        aaCGACAAGTAGAAACACGCAAGAGGTTAATGTTTGCGGCACAATATTTCGTTGAGTCTGTCCACCATAGAGTCATGTATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATg                                                                                 <  1:372310/141‑1 (MQ=255)
                                                         aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTTTGAttg                                                                                <  2:682528/141‑3 (MQ=255)
                                                         aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                                >  1:198958/1‑141 (MQ=255)
                                                         aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                                >  1:449959/1‑141 (MQ=255)
                                                          cGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGt                                                                               >  1:348248/1‑141 (MQ=255)
                                                          cGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGt                                                                               >  1:709379/1‑141 (MQ=255)
                                                           caCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTa                                                                              <  1:733744/140‑1 (MQ=255)
                                                                gTAGAAACGAGGTAAAAACCAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGGATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgggg                                                                         >  1:441947/1‑140 (MQ=255)
                                                                  aGAAATGAGGTAAAGACAAATGCTAGAGGCCAAAATTTTATGGGGTTCCGGACATCATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTgg                                                                       <  2:754922/141‑1 (MQ=255)
                                                                     aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAAGTTTTAGTTCTGTGTGATGGTATGGGCGG‑GCa                                                                   >  1:292675/1‑141 (MQ=255)
                                                                     aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCa                                                                   >  1:39283/1‑141 (MQ=255)
                                                                      aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCAt                                                                  <  2:688366/141‑1 (MQ=255)
                                                                      aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTGAGTTCTGTGTGATGGTATGGGTGG‑CCAt                                                                  <  1:567333/141‑1 (MQ=255)
                                                                              aaGACAAATGCAGTAGGCACAATTTTTGACCGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATCTGTTATAATCAAACTAATCAACAACTTTTAGATCTGTGTGATGGTATGGGTGG‑CCATAAAGCacg                                                          <  2:217281/141‑3 (MQ=255)
                                                                                 aCAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGgaga                                                       >  1:135133/1‑141 (MQ=255)
                                                                                      tGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAGCGAGTCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTg                                                  <  1:738975/141‑1 (MQ=255)
                                                                                       gCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAACCAGGATAAGTTGc                                                 >  2:287437/1‑141 (MQ=255)
                                                                                        cTAGAGGCACAATTTTTTTCTTTTACTGGACAACAAATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCa                                                >  2:503824/1‑141 (MQ=255)
                                                                                          agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCTTAAATCAGGAGAAGTTGCAAg                                              <  1:405285/141‑1 (MQ=255)
                                                                                           gcggCACAATTTTTTCCTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGt                                             <  1:92872/139‑1 (MQ=255)
                                                                                            agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTa                                            >  2:139033/1‑141 (MQ=255)
                                                                                             ggCACAATTTGTTACTGATACGGAACAACATATAAATAAGAATGAAGATGCGGGTAGTATTTTTTATAAGCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTaa                                           <  2:772516/141‑1 (MQ=255)
                                                                                                    ttttttACTAATACTGGACAACATAGAGATAAGAATAAAGAAGCGGGTGGTATTTTTTCTAATCAAACTAATCAACAACTTTTAGTTCGGTGTGATGGTATGAGTGG‑CCATAAAGAGGAAGAAGTTGCAGATTAATTTGtt                                    >  2:339339/1‑141 (MQ=255)
                                                                                                       tttACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                                 >  1:191975/1‑141 (MQ=255)
                                                                                                            tGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATgg                            >  1:320551/1‑140 (MQ=255)
                                                                                                             gATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAg                           >  1:183831/1‑141 (MQ=255)
                                                                                                                   ggACAACAAAGGGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTTTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGaaa                     <  2:275265/141‑1 (MQ=255)
                                                                                                                      caacaaAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATcc                  <  2:649914/141‑1 (MQ=255)
                                                                                                                         aatatAGATTAGGATGAAGGTGCTGGTTGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCCTAAAGCGGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt               <  2:539923/140‑1 (MQ=255)
                                                                                                                           tttagaTATGACTGAAGTTGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt             <  1:4887/138‑1 (MQ=255)
                                                                                                                           tatagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt             >  2:160162/4‑141 (MQ=255)
                                                                                                                           taaagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCAGTCAAGTGGTGGGGGTGGAGCAT‑CGCCAGTATAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt             <  1:185996/138‑1 (MQ=255)
                                                                                                                                aTAAGAATGAAGATGCGGGTGGTATTTTTTATAAACAACCTAATCAACAACATTGATTTCTGTGAAATGTTAGGGGTGG‑CCATAAAGCAGGAGAAGTTTCATGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg        <  2:458793/141‑1 (MQ=255)
                                                                                                                                aTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg        <  2:726982/141‑1 (MQ=255)
                                                                                                                                  aaGAATGTAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg      <  1:456109/141‑1 (MQ=255)
                                                                                                                                  aaGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg      >  2:239656/1‑141 (MQ=255)
                                                                                                                                   aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg     <  2:263436/141‑1 (MQ=255)
                                                                                                                                   aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg     <  1:659909/141‑1 (MQ=255)
                                                                                                                                   aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg     <  1:56459/141‑1 (MQ=255)
                                                                                                                                      aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaa  >  1:415140/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                           
ACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGG‑CCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAA  >  USA300TCH1516_ALE/1231271‑1231545

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: ‑