Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 895,935 T→C 9.4% 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.4% 675.7 / 48.3 234intergenic (+134/+9)metQ_2/Int‑Tn_1Methionine‑binding lipoprotein MetQ/Transposase from transposon Tn916
Reads supporting (aligned to +/- strand):  ref base T (100/112);  new base C (8/14);  total (108/126)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 235 total aligned reads displayed.

ATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTT  >  USA300TCH1516_ALE/895848‑896025
                                                                                        |                                                                                          
   tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcg                                                                                      >  2:3473444/1‑91 (MQ=255)
   tttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcg                                                                                      <  2:1170005/91‑1 (MQ=255)
    ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcgc                                                                                     <  1:2180199/91‑1 (MQ=255)
    ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAgcgc                                                                                     >  2:3566182/1‑91 (MQ=255)
    ttGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAgcgt                                                                                     <  1:1049143/91‑2 (MQ=255)
     tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                    >  2:3296505/1‑91 (MQ=255)
     tgtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCa                                                                                    >  1:1785175/1‑91 (MQ=255)
      gtgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCAt                                                                                   >  1:948098/1‑91 (MQ=255)
       tgGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAAATCTGAGCGTAtt                                                                                  <  2:491481/91‑1 (MQ=255)
        ggTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAAACAAAAAAATACGCTCAGATTTGAGCGCAtt                                                                                  <  2:587601/91‑1 (MQ=255)
        ggTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATGCGCTCAAATCTGAGCGTAttt                                                                                 <  1:2253193/91‑2 (MQ=38)
           aTTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATcc                                                                              >  1:1466577/1‑91 (MQ=255)
             tGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCaa                                                                            <  2:1491844/91‑1 (MQ=255)
             tGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCaa                                                                            >  2:1705560/1‑91 (MQ=255)
               aaaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGt                                                                          <  1:2248712/91‑1 (MQ=255)
                aaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtt                                                                         >  1:612378/1‑91 (MQ=255)
                aaaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtt                                                                         >  1:2016109/1‑91 (MQ=255)
                 aaTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATAAAAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGttt                                                                        <  2:2888561/91‑1 (MQ=255)
                  aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                       >  1:1279350/1‑91 (MQ=255)
                  aTAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGtttt                                                                       <  2:2019759/91‑1 (MQ=255)
                   tAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                      <  2:3790755/91‑1 (MQ=255)
                   aaaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                      <  2:147098/90‑1 (MQ=255)
                   aaaaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTa                                                                      <  2:2095311/90‑1 (MQ=255)
                     aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  1:1685076/1‑91 (MQ=255)
                     aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  1:3502632/1‑91 (MQ=255)
                     aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  1:3696273/1‑91 (MQ=255)
                     aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  2:672989/1‑91 (MQ=255)
                     aaaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTAcc                                                                    >  2:2261184/1‑91 (MQ=255)
                      aaaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCt                                                                   <  2:2699358/91‑1 (MQ=255)
                       aaaTGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  <  1:3026327/91‑1 (MQ=255)
                       aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  <  2:3514582/91‑1 (MQ=255)
                       aaaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTg                                                                  <  1:1395452/91‑1 (MQ=255)
                        aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 <  2:341074/91‑1 (MQ=255)
                        aaGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGt                                                                 <  1:2175585/91‑1 (MQ=255)
                           ggCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCa                                                              <  1:2623031/91‑1 (MQ=255)
                           ggCACAAAAGGGGCAGTATATGCAAATTATCGAATAATAAAAAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCa                                                              <  2:2751564/91‑1 (MQ=255)
                             cacaAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAg                                                            >  2:1015664/1‑91 (MQ=255)
                               caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                          >  1:813876/1‑91 (MQ=255)
                               caAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGtt                                                          <  1:2053732/91‑1 (MQ=255)
                                 aaaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACTCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                        >  1:3089258/1‑91 (MQ=255)
                                 aaaGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTg                                                        <  1:1717245/91‑1 (MQ=255)
                                   aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                      <  2:429009/91‑1 (MQ=255)
                                   aGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGtt                                                      <  2:428673/91‑1 (MQ=255)
                                    ggggCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtc                                                     <  1:3728841/91‑1 (MQ=255)
                                       gCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatc                                                  >  1:1688018/1‑91 (MQ=255)
                                        cAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTtcatca                                                 >  1:1000500/1‑91 (MQ=255)
                                         aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                                >  2:797367/1‑59 (MQ=255)
                                         aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTAt                                                                                <  1:799249/59‑1 (MQ=255)
                                         aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                                >  1:1082078/1‑91 (MQ=255)
                                         aGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTAGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTcatcat                                                >  1:514476/1‑91 (MQ=255)
                                            atatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATc                                             <  2:3703005/91‑1 (MQ=255)
                                             tatGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATcc                                            <  1:1975581/91‑1 (MQ=255)
                                               tGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCtt                                          <  2:1452315/91‑1 (MQ=255)
                                                gCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTa                                         >  2:890717/1‑91 (MQ=255)
                                                 cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                        <  2:742873/91‑1 (MQ=255)
                                                 cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                        <  2:563521/91‑1 (MQ=255)
                                                 cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                        <  2:3667402/91‑1 (MQ=255)
                                                 cAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTAt                                        >  2:2227912/1‑91 (MQ=255)
                                                   aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                      <  1:482609/91‑1 (MQ=255)
                                                   aaTTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATcc                                      <  1:3173423/91‑1 (MQ=255)
                                                       aTCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCt                                  <  2:306637/91‑1 (MQ=255)
                                                        tCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTg                                 <  2:2133446/91‑1 (MQ=255)
                                                          gAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGtt                               >  2:3860144/1‑91 (MQ=255)
                                                          gAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGtt                               >  2:249310/1‑91 (MQ=255)
                                                           aaTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTc                              <  1:2678073/91‑1 (MQ=255)
                                                            aTTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCa                             <  2:366051/91‑1 (MQ=255)
                                                             ttttAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAg                            <  1:1886966/91‑1 (MQ=255)
                                                                tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                         <  2:1779734/91‑1 (MQ=255)
                                                                tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                         <  2:1255465/91‑1 (MQ=255)
                                                                tAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTaa                         <  2:1722281/91‑1 (MQ=255)
                                                                  aaaCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAAc                        >  2:3769039/1‑91 (MQ=255)
                                                                    aCAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGt                      >  2:3159207/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)
                                                                       aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                   >  2:1802951/1‑91 (MQ=255)
                                                                       aaaaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgag                   <  1:1629629/91‑1 (MQ=255)
                                                                        taaaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTgaga                  >  1:2492046/2‑91 (MQ=255)
                                                                          aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                >  2:1910498/1‑91 (MQ=255)
                                                                          aaaTACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAAt                >  2:1598997/1‑91 (MQ=255)
                                                                            atgcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaa              >  1:669719/4‑91 (MQ=255)
                                                                             tgcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             >  2:2312475/3‑91 (MQ=255)
                                                                             tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             <  2:714681/91‑1 (MQ=255)
                                                                             tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             <  1:3435180/91‑1 (MQ=255)
                                                                             tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             <  1:170477/91‑1 (MQ=255)
                                                                             tACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATaaa             <  1:1144589/91‑1 (MQ=255)
                                                                              gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt            >  1:1743596/2‑91 (MQ=255)
                                                                              gcGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAAt            >  1:2393864/2‑91 (MQ=255)
                                                                               cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           <  2:1839692/91‑1 (MQ=255)
                                                                               cGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           >  2:324427/1‑91 (MQ=255)
                                                                               cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           >  1:3796966/1‑91 (MQ=255)
                                                                               cGCTCAAATCTGAGCGTATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATg           >  1:1925908/1‑91 (MQ=255)
                                                                                gCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATgg          <  2:2991562/91‑1 (MQ=255)
                                                                                 cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa         <  1:2275032/91‑1 (MQ=255)
                                                                                 cTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGa         >  1:2748686/1‑91 (MQ=255)
                                                                                  tCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGaa        >  1:362673/1‑91 (MQ=255)
                                                                                  tCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGaa        <  1:67949/91‑1 (MQ=255)
                                                                                    aGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg      <  1:2787237/91‑1 (MQ=255)
                                                                                    aGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAg      <  1:2858962/91‑1 (MQ=255)
                                                                                     gATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAgt     <  1:3545211/91‑1 (MQ=255)
                                                                                        |                                                                                          
ATTTTTGTGGTATTGAAAATAAAAAAGGGCACAAAAGGGGCAGTATATGCAAATTATCGAATTTTAAACAAAAAAATACGCTCAGATTTGAGCGCATTATCCAAGTTTTACCTGTTCAAGTTTGTTCATCATATCCTTATCCATCTGTTCAGTAACGTGAGAATAAATGGAAAGTGTT  >  USA300TCH1516_ALE/895848‑896025

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: