Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,227 C→G 9.2% intergenic (‑198/+51) ptlE ← / ← mnhG1 Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE942,2270CG 9.2% 737.3 / 51.0 218intergenic (‑198/+51)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base C (106/92);  new base G (9/11);  total (115/103)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

Only 100 of 236 total aligned reads displayed.

GGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCG  >  USA300TCH1516_ALE/942138‑942317
                                                                                         |                                                                                          
gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                           <  2:2825898/91‑1 (MQ=255)
gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                           >  2:2083041/1‑91 (MQ=255)
gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                           >  2:1518998/1‑91 (MQ=255)
   cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                        <  2:918691/91‑1 (MQ=255)
      aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                     <  1:1755380/91‑1 (MQ=255)
        cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActg                                                                                   <  1:1847696/91‑1 (MQ=18)
        cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATCTGGAATCATCTTActg                                                                                   >  2:3396895/1‑91 (MQ=255)
         tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                  <  2:3408712/91‑1 (MQ=18)
         tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                  <  2:2018719/91‑1 (MQ=18)
         tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                  <  1:427434/91‑1 (MQ=18)
          gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                 >  2:2735366/1‑91 (MQ=255)
          gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                 >  2:1028095/1‑91 (MQ=255)
           tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                >  1:1103463/1‑91 (MQ=255)
           tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                >  2:200451/1‑91 (MQ=255)
           tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                >  1:624472/1‑91 (MQ=255)
           aagacAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                <  2:945642/87‑1 (MQ=255)
             gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                              <  2:1568508/91‑1 (MQ=255)
             gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                              <  2:3799528/91‑1 (MQ=255)
                  aGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttg                                                                         <  1:1023401/91‑1 (MQ=255)
                   gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                                        <  1:3413883/91‑1 (MQ=255)
                     aaaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgttg                                                                      <  1:2147698/91‑1 (MQ=255)
                      aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                     >  2:1245220/1‑91 (MQ=255)
                      aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                     >  2:3089040/1‑91 (MQ=255)
                      aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                     >  2:2123709/1‑91 (MQ=255)
                         aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                  >  1:2105255/1‑91 (MQ=255)
                         aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                  >  1:2699011/1‑91 (MQ=255)
                         aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                  >  1:3278777/1‑91 (MQ=255)
                          tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                 >  1:772604/1‑91 (MQ=255)
                          tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                 >  1:2443777/1‑91 (MQ=255)
                          tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                 >  1:1563311/1‑91 (MQ=255)
                           cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                >  1:2938944/1‑91 (MQ=255)
                           cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                <  1:3520433/91‑1 (MQ=255)
                           cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                >  2:14003/1‑91 (MQ=255)
                               gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                            >  1:3195858/1‑91 (MQ=255)
                               gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                            >  2:24635/1‑91 (MQ=255)
                               gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                            <  1:1647655/91‑1 (MQ=255)
                               gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                            >  1:951196/1‑91 (MQ=255)
                               gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                            <  1:785567/91‑1 (MQ=255)
                                  aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                                         >  2:2376218/1‑91 (MQ=255)
                                     gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                                      >  2:219131/1‑91 (MQ=255)
                                       ggTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtttt                                                    >  1:2015914/1‑91 (MQ=255)
                                        gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCCTGATTCTg                                                   >  2:1703/1‑91 (MQ=255)
                                        gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                                                   >  1:3615597/1‑91 (MQ=255)
                                         tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                  <  1:2399873/91‑1 (MQ=255)
                                         tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                  >  2:741797/1‑91 (MQ=255)
                                          aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                                 >  1:3540336/1‑91 (MQ=255)
                                           aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                                <  1:651873/91‑1 (MQ=255)
                                            cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                                               >  2:3463001/1‑91 (MQ=255)
                                                 tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                          >  2:685618/1‑91 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         <  2:1993772/91‑1 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         <  1:601818/91‑1 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         <  2:2694668/91‑1 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         >  2:3724020/1‑91 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         <  1:630534/91‑1 (MQ=255)
                                                  cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                         <  1:446948/91‑1 (MQ=255)
                                                   aaaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCt                                        <  2:2331328/91‑1 (MQ=255)
                                                    taaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                       <  2:3932765/90‑1 (MQ=255)
                                                    aaaGAAAACAGCAGTAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                                      >  2:3883361/1‑91 (MQ=255)
                                                     aaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                                      <  2:2703650/91‑1 (MQ=255)
                                                        aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                                   <  2:242772/91‑1 (MQ=255)
                                                        aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                                   >  1:1409266/1‑91 (MQ=255)
                                                        aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                                   <  1:3005678/91‑1 (MQ=255)
                                                         aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                  >  1:2606810/1‑91 (MQ=255)
                                                          aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                 >  1:566119/1‑91 (MQ=255)
                                                           aCAGCAGTAAGATCATTCCTTATTTGGATTCTTCATACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtttt                                <  2:3342672/91‑1 (MQ=255)
                                                           aCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtttt                                <  2:3458941/91‑1 (MQ=255)
                                                            cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                               >  1:3073034/1‑91 (MQ=255)
                                                             agcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                              >  1:1213524/1‑91 (MQ=18)
                                                             agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                              >  1:1675862/1‑91 (MQ=255)
                                                             agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                              <  2:3580040/91‑1 (MQ=255)
                                                             agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                              <  1:2760071/91‑1 (MQ=255)
                                                              gcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                             >  2:429798/1‑91 (MQ=21)
                                                              gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                             >  1:334443/1‑91 (MQ=255)
                                                              gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                             <  2:2253999/91‑1 (MQ=255)
                                                              gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                             <  1:2033988/91‑1 (MQ=255)
                                                              gcagTAAAATCATTCATAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATCCATGATTTTGTTATATAATCTACAATTTtgtg                             >  1:3374304/1‑91 (MQ=255)
                                                               cagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                            >  2:2279246/1‑91 (MQ=21)
                                                                agTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                           >  2:1852440/1‑91 (MQ=21)
                                                                agTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                           <  2:636321/91‑1 (MQ=255)
                                                                 gTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCt                          >  1:3769195/1‑91 (MQ=255)
                                                                  tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                         <  1:2927425/91‑1 (MQ=255)
                                                                   aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt                        <  2:599959/91‑1 (MQ=255)
                                                                    aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                       <  1:3718355/91‑1 (MQ=255)
                                                                        gaTTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                   >  1:132177/2‑91 (MQ=37)
                                                                        cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                   >  1:3457228/1‑91 (MQ=255)
                                                                         aTTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTc                  >  2:2669357/1‑91 (MQ=255)
                                                                          ttCCTAATTTGGAGTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTATTATATAATCTACAATTTTGTGTCTTTTAAGTCt                 <  2:1161492/91‑1 (MQ=255)
                                                                          ttCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCt                 >  2:1354313/1‑91 (MQ=255)
                                                                          ttCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCt                 >  2:1986850/1‑91 (MQ=255)
                                                                           tCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTTAGTCtt                >  2:10666/1‑91 (MQ=255)
                                                                              tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg             <  1:3381834/91‑1 (MQ=255)
                                                                              tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg             <  2:2854167/91‑1 (MQ=255)
                                                                               aaTTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGa            >  2:141005/1‑91 (MQ=255)
                                                                               aaTTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGa            >  2:2983459/1‑91 (MQ=255)
                                                                                   tGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAtt        <  1:223753/91‑1 (MQ=255)
                                                                                    ggAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAttt       >  1:2660096/1‑91 (MQ=255)
                                                                                      aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa     <  2:839674/91‑1 (MQ=255)
                                                                                      aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa     <  1:1901769/91‑1 (MQ=255)
                                                                                       atcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCAt    >  2:567125/1‑91 (MQ=255)
                                                                                        tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc   <  1:2306390/91‑1 (MQ=255)
                                                                                         |                                                                                          
GGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCG  >  USA300TCH1516_ALE/942138‑942317

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