Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 895,932 G→A 9.7% intergenic (+131/+12) metQ_2 → / ← Int‑Tn_1 Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE895,9320GA 9.7% 706.5 / 45.6 236intergenic (+131/+12)metQ_2/Int‑Tn_1Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916
Reads supporting (aligned to +/- strand):  ref base G (102/111);  new base A (8/15);  total (110/126)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 240 total aligned reads displayed.

GGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGT  >  USA300TCH1516_ALE/895842‑896022
                                                                                           |                                                                                          
ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAg                                                                                            >  1:3434412/1‑91 (MQ=255)
ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAg                                                                                            >  1:1621096/1‑91 (MQ=255)
  gggtATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATAGAATTTTAAACAAAAAAATACGCTGAGAt                                                                                          >  2:2028927/4‑91 (MQ=255)
  ggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACGAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAAt                                                                                          <  2:3829970/91‑1 (MQ=255)
         tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                   >  2:3178598/1‑91 (MQ=255)
         tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcg                                                                                   <  2:1170005/91‑1 (MQ=255)
          ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcgc                                                                                  <  1:2179962/91‑1 (MQ=255)
          ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcgt                                                                                  <  1:1049143/91‑2 (MQ=255)
           tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                 <  2:2865614/91‑1 (MQ=255)
           tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                 >  1:1785175/1‑91 (MQ=255)
           tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcgta                                                                                 <  2:192603/91‑3 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                >  1:431205/1‑91 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                >  1:948098/1‑91 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                >  2:2732428/1‑91 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                >  1:1924045/1‑91 (MQ=255)
             tgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAtt                                                                               >  1:2133446/1‑91 (MQ=255)
              ggTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAAACAAAAAAATACGCTCAGATTTGAGCGCAtt                                                                               <  2:587601/91‑1 (MQ=255)
                 aTTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                           <  2:3439236/91‑1 (MQ=255)
                   tGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCaa                                                                         >  2:1705715/1‑91 (MQ=255)
                   tGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCaa                                                                         <  2:1491844/91‑1 (MQ=255)
                     aaaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGt                                                                       <  1:2248222/91‑1 (MQ=255)
                      aaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtt                                                                      >  1:612378/1‑91 (MQ=255)
                        aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                    >  1:1279729/1‑91 (MQ=255)
                         tAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                   <  2:3790755/91‑1 (MQ=255)
                         aaaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                   <  2:147098/90‑1 (MQ=255)
                           aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                 >  2:672989/1‑91 (MQ=255)
                           aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                 <  2:3696273/91‑1 (MQ=255)
                           aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                 >  1:1685076/1‑91 (MQ=255)
                           aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                 >  2:2261184/1‑91 (MQ=255)
                             aaaTGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                               <  1:3026327/91‑1 (MQ=255)
                             aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                               <  1:1395452/91‑1 (MQ=255)
                              aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                              >  1:225290/1‑91 (MQ=255)
                              aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                              >  2:3350695/1‑91 (MQ=255)
                                gggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTc                                                            >  1:1509799/1‑91 (MQ=255)
                                 ggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCa                                                           <  1:2623031/91‑1 (MQ=255)
                                    acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                        >  1:963136/1‑91 (MQ=255)
                                     caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                       >  1:813876/1‑91 (MQ=255)
                                     caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                       >  2:304237/1‑91 (MQ=255)
                                       aaaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACTCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                     >  1:3089258/1‑91 (MQ=255)
                                       aaaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                     <  1:1717245/91‑1 (MQ=255)
                                       aaaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                     >  1:2848903/1‑91 (MQ=255)
                                         aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                   <  2:428848/91‑1 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  <  2:3527559/91‑1 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  <  1:3475163/91‑1 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  <  2:456275/91‑1 (MQ=255)
                                           gggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtca                                                 >  1:3686928/1‑91 (MQ=255)
                                             gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                               >  1:3537494/1‑91 (MQ=255)
                                             gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                               >  1:1688018/1‑91 (MQ=255)
                                              cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                              >  1:1000500/1‑91 (MQ=255)
                                               aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGTGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                             >  2:3265263/1‑91 (MQ=255)
                                               aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                             <  1:797367/59‑1 (MQ=255)
                                               aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                             >  2:797367/1‑59 (MQ=255)
                                               aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                             >  1:1082078/1‑91 (MQ=255)
                                                  atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                          <  2:3490928/91‑1 (MQ=255)
                                                     tGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCtt                                       <  2:1452315/91‑1 (MQ=255)
                                                     tGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCtt                                       >  2:1591239/1‑91 (MQ=255)
                                                      gCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTa                                      <  1:3699891/91‑1 (MQ=255)
                                                      gCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTa                                      >  2:890717/1‑91 (MQ=255)
                                                       cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     >  2:2227912/1‑91 (MQ=255)
                                                       cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     <  2:563521/91‑1 (MQ=255)
                                                       cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     <  2:742873/91‑1 (MQ=255)
                                                       cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     <  2:3667402/91‑1 (MQ=255)
                                                         aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                   <  1:3173423/91‑1 (MQ=255)
                                                          aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                  >  1:199179/1‑91 (MQ=255)
                                                             aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCt                               <  1:672989/91‑1 (MQ=255)
                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              <  1:667607/91‑1 (MQ=255)
                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              <  2:2133446/91‑1 (MQ=255)
                                                              tCGAAGTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              <  2:3148538/91‑1 (MQ=255)
                                                               cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                             >  1:2035730/1‑91 (MQ=255)
                                                               cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                             <  2:372277/91‑1 (MQ=255)
                                                                gAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGtt                            >  2:2060299/1‑91 (MQ=255)
                                                                   ttttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAg                         <  1:1886966/91‑1 (MQ=255)
                                                                    tttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGt                        <  2:3879733/91‑1 (MQ=255)
                                                                    tttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGt                        <  1:487194/91‑1 (MQ=255)
                                                                      tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                      <  2:1722281/91‑1 (MQ=255)
                                                                      tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                      <  2:1255465/91‑1 (MQ=255)
                                                                      tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                      <  2:1779734/91‑1 (MQ=255)
                                                                           cAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATGCTTATCCATCTGTTCAGTAACGTg                  >  2:3616405/1‑91 (MQ=255)
                                                                           cAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                  >  1:144183/1‑91 (MQ=255)
                                                                            aaaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTga                 <  1:1213671/91‑1 (MQ=255)
                                                                            aaaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTga                 >  1:696639/1‑91 (MQ=255)
                                                                             aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGCTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgtg                <  1:3735947/91‑3 (MQ=255)
                                                                             aaaaaaTACGCCCAGATTTGCGCGCATTATCCACGTTTTACCTGTTCAAGTTTGCACATCCTATCCTCATCCCTCTGTTCCGTGACGTgcg                >  2:1803008/1‑89 (MQ=255)
                                                                              taaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgaga               >  1:2492046/2‑91 (MQ=255)
                                                                                aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt             >  2:1598997/1‑91 (MQ=255)
                                                                                   tgcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa          >  2:2312475/3‑91 (MQ=255)
                                                                                   tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa          <  1:1144589/91‑1 (MQ=255)
                                                                                   tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa          <  2:714681/91‑1 (MQ=255)
                                                                                    gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt         >  1:1743596/2‑91 (MQ=255)
                                                                                    gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt         >  1:2393864/2‑91 (MQ=255)
                                                                                    aCGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGACAATAAAt         >  1:1930529/1‑91 (MQ=255)
                                                                                     cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg        >  1:16933/1‑91 (MQ=255)
                                                                                     cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg        <  2:1839692/91‑1 (MQ=255)
                                                                                      gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       <  2:2991523/91‑1 (MQ=255)
                                                                                      gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       <  2:2522553/91‑1 (MQ=255)
                                                                                       cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa      >  1:2748686/1‑91 (MQ=255)
                                                                                       cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa      <  1:2275032/91‑1 (MQ=255)
                                                                                        tCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTATTCAGTAACGTGAGAATAAATGGaa     >  1:3600613/1‑91 (MQ=255)
                                                                                          aGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg   <  1:2787237/91‑1 (MQ=255)
                                                                                          aGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg   <  1:2677189/91‑1 (MQ=255)
                                                                                           |                                                                                          
GGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGT  >  USA300TCH1516_ALE/895842‑896022

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