Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 895,932 G→A 9.0% 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.0% 746.1 / 40.6 243intergenic (+131/+12)metQ_2/Int‑Tn_1Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916
Reads supporting (aligned to +/- strand):  ref base G (90/131);  new base A (8/14);  total (98/145)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 251 total aligned reads displayed.

GGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGT  >  USA300TCH1516_ALE/895842‑896022
                                                                                          |                                                                                          
ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAg                                                                                            >  1:1182125/1‑91 (MQ=255)
ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAg                                                                                            >  1:3791019/1‑91 (MQ=255)
ggggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAg                                                                                            >  1:1953145/1‑91 (MQ=255)
 gggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGa                                                                                           <  1:3938437/91‑1 (MQ=255)
  ggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGAt                                                                                          <  1:98038/91‑1 (MQ=255)
  ggTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGAt                                                                                          <  1:1365271/91‑1 (MQ=255)
    ttatttatGTAGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTATAAACAAAAAAATACGCTCAGAttt                                                                                        <  2:1701977/91‑1 (MQ=255)
       tttttGTAGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAg                                                                                     <  2:641031/91‑1 (MQ=255)
        ttttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgc                                                                                    >  2:1336481/1‑91 (MQ=255)
        ttttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgc                                                                                    >  2:2668920/1‑91 (MQ=255)
        ttttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgc                                                                                    <  2:1869699/91‑1 (MQ=255)
         tttGTGTTATTTAAAATGAAAAATGGCATAAAAGGGGAAGTATATGAAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                   <  2:1211597/91‑1 (MQ=255)
         tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAGATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                   <  2:807110/91‑1 (MQ=255)
         tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcg                                                                                   <  2:91933/91‑1 (MQ=255)
          ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcgc                                                                                  <  1:261421/91‑1 (MQ=255)
          ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcgt                                                                                  <  2:3309290/91‑2 (MQ=255)
           tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATCTGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                 >  2:3837500/1‑91 (MQ=255)
           tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                 <  2:3231869/91‑1 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAGGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAAGCTGAgcgtat                                                                                <  2:486609/91‑4 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                <  1:3285252/91‑1 (MQ=255)
            gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAGCGTAt                                                                                <  2:1059860/91‑1 (MQ=255)
            gtgGTATTGAAAATAAAAAAAGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                <  2:2911325/91‑1 (MQ=255)
                 aaagAAAATAACACAGGGCACAAAAGGGACAGTATATGCAAATCATCGAATATTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                           <  2:1039980/88‑1 (MQ=255)
                 aTTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                           <  1:1441899/91‑1 (MQ=255)
                     aaaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGt                                                                       <  2:2132698/91‑1 (MQ=255)
                      aaaaaaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtt                                                                      <  2:3343105/91‑1 (MQ=255)
                       aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGttt                                                                     <  1:2465940/91‑1 (MQ=255)
                         tAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                   <  2:2345026/91‑1 (MQ=255)
                         tAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                   >  2:3112737/1‑91 (MQ=255)
                           aaaaaGGGCACAAAAGGGCCAGTATATGCAAATTAACGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                 >  2:182229/1‑91 (MQ=255)
                             aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                               <  1:1657678/91‑1 (MQ=255)
                             aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                               <  1:2999167/91‑1 (MQ=255)
                              aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                              >  1:3631421/1‑91 (MQ=255)
                               aGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGtt                                                             <  2:1440996/91‑1 (MQ=255)
                                 ggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCa                                                           <  1:126308/91‑1 (MQ=255)
                                  gCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCaa                                                          >  1:2047606/1‑91 (MQ=255)
                                   cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                         <  1:2794079/91‑1 (MQ=255)
                                   cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                         <  1:2298177/91‑1 (MQ=255)
                                    acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                        >  1:258467/1‑91 (MQ=255)
                                    acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                        <  2:2773758/91‑1 (MQ=255)
                                     caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                       >  1:708292/1‑91 (MQ=255)
                                       aaaGGGGCAGTACATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                     >  1:582125/1‑91 (MQ=255)
                                         aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                   <  2:1805523/91‑1 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  >  2:2775465/1‑91 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  <  1:1052684/91‑1 (MQ=255)
                                          ggggCAGTATATGCAAATTATCGAAGTATAAACAAAACAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                  <  2:1177496/91‑1 (MQ=255)
                                            ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                >  2:1891331/1‑91 (MQ=255)
                                            ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                >  2:1951525/1‑91 (MQ=255)
                                            ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                <  2:2388329/91‑1 (MQ=255)
                                             gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                               >  2:3849749/1‑91 (MQ=255)
                                             gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                               <  1:809007/91‑1 (MQ=255)
                                              cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                              <  1:715438/91‑1 (MQ=255)
                                               aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                             <  2:2237024/91‑1 (MQ=255)
                                                gTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCata                                            >  1:2694226/1‑91 (MQ=255)
                                                  atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                          <  1:420769/91‑1 (MQ=255)
                                                  atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                          <  1:2543189/91‑1 (MQ=255)
                                                    atGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCt                                        >  1:2868879/1‑91 (MQ=255)
                                                    atGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCt                                        <  2:2542004/91‑1 (MQ=255)
                                                      gCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTa                                      <  1:684666/91‑1 (MQ=255)
                                                       cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     <  1:2838118/91‑1 (MQ=255)
                                                       cAAATTATCGAATGATAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                     <  2:571064/91‑1 (MQ=255)
                                                         aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                   <  1:3528094/91‑1 (MQ=255)
                                                          aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                  <  1:1881921/91‑1 (MQ=255)
                                                          aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                  <  2:1294775/91‑1 (MQ=255)
                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              <  2:3732424/91‑1 (MQ=255)
                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              >  1:2316055/1‑91 (MQ=255)
                                                              tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                              <  2:960734/91‑1 (MQ=255)
                                                               cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                             >  2:3553199/1‑91 (MQ=255)
                                                               cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                             >  1:2806021/1‑91 (MQ=255)
                                                               cGAATTTTAAAAAAAAACATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCAGTTCAAGTTTGTTAATAATATCCTTATCCATCTGt                             >  1:2130135/1‑91 (MQ=255)
                                                                gAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGtt                            <  2:26253/91‑1 (MQ=255)
                                                                  aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                          <  2:3213557/91‑1 (MQ=255)
                                                                  aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                          <  1:3187078/91‑1 (MQ=255)
                                                                  aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                          >  1:3903414/1‑91 (MQ=255)
                                                                  aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                          <  2:2484598/91‑1 (MQ=255)
                                                                   ttttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAg                         >  1:2357442/1‑91 (MQ=255)
                                                                         aCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                   >  2:2272212/1‑91 (MQ=255)
                                                                          cAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                  <  2:1127838/91‑1 (MQ=255)
                                                                           aaaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTga                 <  2:363525/91‑1 (MQ=255)
                                                                            aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                <  2:1371062/91‑1 (MQ=255)
                                                                            aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                >  1:1902945/1‑91 (MQ=255)
                                                                            aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                >  1:151082/1‑91 (MQ=255)
                                                                              aaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGaa              <  1:3747013/91‑1 (MQ=255)
                                                                               aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt             <  2:444197/91‑1 (MQ=255)
                                                                                  tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa          <  1:752490/91‑1 (MQ=255)
                                                                                   gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt         >  1:1776365/2‑91 (MQ=255)
                                                                                   aCGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt         <  1:460210/91‑1 (MQ=255)
                                                                                    cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg        >  2:1613915/1‑91 (MQ=255)
                                                                                    cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg        <  1:3802784/91‑1 (MQ=255)
                                                                                    cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg        >  1:1731692/1‑91 (MQ=255)
                                                                                     gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       <  1:3323235/91‑1 (MQ=255)
                                                                                     gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       <  2:1159968/91‑1 (MQ=255)
                                                                                     gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAATTTGTTCATCATATCCTTATCCATCGGTTCAGTAACGTGAGAATAAATgg       >  2:3665164/1‑91 (MQ=255)
                                                                                     gCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       >  1:1761222/1‑91 (MQ=255)
                                                                                     gCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg       >  1:2162487/1‑91 (MQ=255)
                                                                                      cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATggc      >  2:202579/1‑90 (MQ=255)
                                                                                      cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa      >  2:2083269/1‑91 (MQ=255)
                                                                                      cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa      <  1:208672/91‑1 (MQ=255)
                                                                                         aGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg   <  1:1873607/91‑1 (MQ=255)
                                                                                          gATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt  <  1:350768/91‑1 (MQ=255)
                                                                                          |                                                                                          
GGGGTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGT  >  USA300TCH1516_ALE/895842‑896022

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