Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 895,935 T→C 9.3% intergenic (+134/+9) 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,9350TC 9.3% 687.7 / 48.3 236intergenic (+134/+9)metQ_2/Int‑Tn_1Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916
Reads supporting (aligned to +/- strand):  ref base T (88/125);  new base C (8/14);  total (96/140)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 244 total aligned reads displayed.

GTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTT  >  USA300TCH1516_ALE/895845‑896025
                                                                                          |                                                                                          
gTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGAtt                                                                                            >  2:550067/1‑91 (MQ=255)
 ttatttatGTAGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTATAAACAAAAAAATACGCTCAGAttt                                                                                           <  2:1701977/91‑1 (MQ=255)
     ttttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgc                                                                                       >  2:1336481/1‑91 (MQ=255)
     ttttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgc                                                                                       <  2:1868514/91‑1 (MQ=255)
      tttGTGTTATTTAAAATGAAAAATGGCATAAAAGGGGAAGTATATGAAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                      <  2:1211597/91‑1 (MQ=255)
      tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAGATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                      <  2:807110/91‑1 (MQ=255)
      tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                      <  1:1950508/91‑1 (MQ=255)
       ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcgc                                                                                     <  1:261421/91‑1 (MQ=255)
        tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcgta                                                                                    <  1:2162220/91‑3 (MQ=255)
         gtgGTATTGAATAGAAAAAAGGGTACAAAAGGGGCAGTATATGCAAATTATCGAATTTTTAACAAAAAAATGCGCTCAAATCTGAgcgtat                                                                                   <  2:3017621/91‑4 (MQ=255)
         gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                   >  2:1607932/1‑91 (MQ=255)
         gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                   <  1:3285252/91‑1 (MQ=255)
         gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAGCGTAt                                                                                   <  2:1059860/91‑1 (MQ=255)
            gTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                                <  2:1370936/91‑1 (MQ=255)
            gTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                                >  2:3322987/1‑91 (MQ=255)
             tATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATc                                                                               <  2:1924442/91‑1 (MQ=255)
              aaagAAAATAACACAGGGCACAAAAGGGACAGTATATGCAAATCATCGAATATTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                              <  2:1039980/88‑1 (MQ=255)
                tGAAAATAAAAAATGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCaa                                                                            <  1:3962079/91‑1 (MQ=255)
                  aaaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGt                                                                          <  2:2132698/91‑1 (MQ=255)
                    aaaaaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGttt                                                                        <  2:2472562/88‑1 (MQ=255)
                    aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGttt                                                                        <  1:2465940/91‑1 (MQ=255)
                       aaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAc                                                                     <  2:3842905/91‑1 (MQ=255)
                        aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  1:1766995/1‑91 (MQ=255)
                        aaaaaGGGCACAAAAGGGCCAGTATATGCAAATTAACGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  2:182229/1‑91 (MQ=255)
                         aaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCt                                                                   >  2:3219551/1‑91 (MQ=255)
                          aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  <  1:1657678/91‑1 (MQ=255)
                          aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  >  2:2764644/1‑91 (MQ=255)
                          aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  <  1:2999167/91‑1 (MQ=255)
                           aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 >  1:3631421/1‑91 (MQ=255)
                           aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 >  1:985869/1‑91 (MQ=255)
                           aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 >  1:3839125/1‑91 (MQ=255)
                           aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 >  1:1413452/1‑91 (MQ=255)
                            aGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGtt                                                                <  2:1440996/91‑1 (MQ=255)
                              ggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCa                                                              <  1:126308/91‑1 (MQ=255)
                               gCACATAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCaa                                                             <  2:3202427/91‑1 (MQ=255)
                               gCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCaa                                                             >  1:2047606/1‑91 (MQ=255)
                               gCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCaa                                                             >  2:925335/1‑91 (MQ=255)
                                cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                            <  1:2298177/91‑1 (MQ=255)
                                cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                            >  2:2644515/1‑91 (MQ=255)
                                 acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                           <  2:2773884/91‑1 (MQ=255)
                                 acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                           <  1:1227185/91‑1 (MQ=255)
                                 acaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGt                                                           >  1:258467/1‑91 (MQ=255)
                                  caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                          <  1:380043/91‑1 (MQ=255)
                                  caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                          >  1:708292/1‑91 (MQ=255)
                                    aaaGGGGCAGTACATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                        >  1:582125/1‑91 (MQ=255)
                                       ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                     <  1:1052684/91‑1 (MQ=255)
                                       ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                     <  2:3802511/91‑1 (MQ=255)
                                       ggggCAGTATATGCAAATTATCGAAGTATAAACAAAACAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                     <  2:1177496/91‑1 (MQ=255)
                                         ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                   <  2:2388329/91‑1 (MQ=255)
                                         ggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcat                                                   >  2:1891331/1‑91 (MQ=255)
                                          gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                                  <  1:809007/91‑1 (MQ=255)
                                           cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                 <  1:715438/91‑1 (MQ=255)
                                           cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                 <  2:3723319/91‑1 (MQ=255)
                                            aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                                <  2:2237024/91‑1 (MQ=255)
                                             gTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCata                                               >  1:2694226/1‑91 (MQ=255)
                                               atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                             <  1:2543189/91‑1 (MQ=255)
                                               atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                             <  1:420769/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:2837928/91‑1 (MQ=255)
                                                    cAAATTATCGAATGATAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                        <  2:571064/91‑1 (MQ=255)
                                                      aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                      <  1:3528094/91‑1 (MQ=255)
                                                      aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                      <  1:1869952/91‑1 (MQ=255)
                                                       aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                     <  2:1294775/91‑1 (MQ=255)
                                                       aTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCa                                     <  1:1881906/91‑1 (MQ=255)
                                                          aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCt                                  <  2:3260877/91‑1 (MQ=255)
                                                           tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                 <  2:356114/91‑1 (MQ=255)
                                                           tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                 >  1:2316055/1‑91 (MQ=255)
                                                           tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                 <  2:960734/91‑1 (MQ=255)
                                                            cGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGt                                >  1:2806021/1‑91 (MQ=255)
                                                            cGAATTTTAAAAAAAAACATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCAGTTCAAGTTTGTTAATAATATCCTTATCCATCTGt                                >  1:2130135/1‑91 (MQ=255)
                                                               aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAGGTTTGTTCATCATATCCTTATCCATCTGTTCa                             <  2:3885448/91‑1 (MQ=255)
                                                                ttttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAg                            >  1:2357442/1‑91 (MQ=255)
                                                                  ttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTa                          <  2:3375980/91‑1 (MQ=255)
                                                                     aaCAATAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACg                       <  2:3653651/91‑1 (MQ=255)
                                                                      aCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                      >  2:2272212/1‑91 (MQ=255)
                                                                       cAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTg                     <  2:1127838/91‑1 (MQ=255)
                                                                        aaaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTga                    <  2:2814188/91‑1 (MQ=255)
                                                                         aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                   >  1:151082/1‑91 (MQ=255)
                                                                         aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                   >  1:1902606/1‑91 (MQ=255)
                                                                           aaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGaa                 <  1:3747013/91‑1 (MQ=255)
                                                                            aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                <  2:444197/91‑1 (MQ=255)
                                                                             aaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATa               >  1:3200006/1‑91 (MQ=255)
                                                                               tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             <  1:752490/91‑1 (MQ=255)
                                                                                aCGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt            <  1:460210/91‑1 (MQ=255)
                                                                                 cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           <  2:677588/91‑1 (MQ=255)
                                                                                 cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           >  1:2788900/1‑91 (MQ=255)
                                                                                 cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           >  1:1731692/1‑91 (MQ=255)
                                                                                  gcgcAGATTTGAGCGCATTATCCAAGTTTCACATGTTCAAGTTTGTTCATGATATCCTCGTCCACCTGTTCAGTAACGTGCGAAGAAACgg          >  2:2585960/4‑91 (MQ=255)
                                                                                  gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg          <  2:1159968/91‑1 (MQ=255)
                                                                                  gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg          <  1:3323235/91‑1 (MQ=255)
                                                                                   cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATggc         >  2:202579/1‑90 (MQ=255)
                                                                                   cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa         >  2:2083269/1‑91 (MQ=255)
                                                                                   cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa         <  1:208672/91‑1 (MQ=255)
                                                                                      aGATTTGAGCGCATTATCCAAGTTTTACCTGATCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg      <  2:1721969/91‑1 (MQ=255)
                                                                                       gATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt     <  1:350768/91‑1 (MQ=255)
                                                                                        aTTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtg    <  1:1748439/91‑1 (MQ=255)
                                                                                        aTTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtg    <  1:2242627/91‑1 (MQ=255)
                                                                                         tttGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgtgt   <  1:1791956/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTTATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTT  >  USA300TCH1516_ALE/895845‑896025

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: