Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,230 C→T 9.1% intergenic (‑201/+48) 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,2300CT 9.1% 687.0 / 38.6 218intergenic (‑201/+48)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base C (106/92);  new base T (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 = 9.96e-01

Only 100 of 234 total aligned reads displayed.

GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACT  >  USA300TCH1516_ALE/942140‑942320
                                                                                          |                                                                                          
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           <  2:918691/91‑1 (MQ=255)
  gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                                          <  1:2859545/91‑1 (MQ=255)
    aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        <  1:1755380/91‑1 (MQ=255)
      cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActg                                                                                      <  1:1847696/91‑1 (MQ=18)
      cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                                                      >  2:3396971/1‑91 (MQ=255)
      cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                                                      <  1:3722239/91‑1 (MQ=255)
      cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                                                      >  1:2939744/1‑91 (MQ=255)
       tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                     <  2:2831688/91‑1 (MQ=18)
       tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                     <  2:2018719/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:2599297/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)
            tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTACTGCtgttt                                                                                <  2:2214746/91‑1 (MQ=17)
            tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                                <  2:2309792/91‑1 (MQ=255)
              aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                                              <  2:3368878/91‑1 (MQ=255)
                aGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttg                                                                            <  1:1023401/91‑1 (MQ=255)
                       aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                     >  1:2699011/1‑91 (MQ=255)
                       aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                     >  1:2105255/1‑91 (MQ=255)
                        tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                    >  1:2443777/1‑91 (MQ=255)
                        tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                    >  1:772604/1‑91 (MQ=255)
                        tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                    >  1:1563311/1‑91 (MQ=255)
                         cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                   >  2:3600538/1‑91 (MQ=255)
                         cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                   >  2:1313843/1‑91 (MQ=255)
                         cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                                                   >  1:2938944/1‑91 (MQ=255)
                          aTTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtata                                                                  >  1:3323985/1‑91 (MQ=255)
                             gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                               <  1:2131774/91‑1 (MQ=255)
                             gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                               >  1:951196/1‑91 (MQ=255)
                             gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                               <  1:1647655/91‑1 (MQ=255)
                             gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                               >  1:3195858/1‑91 (MQ=255)
                                aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                                            >  2:2376218/1‑91 (MQ=255)
                                 aaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAGTTATATTCATGa                                                           <  1:3037372/91‑1 (MQ=255)
                                   gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                                         >  1:3497015/1‑91 (MQ=255)
                                   gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                                         >  1:635358/1‑91 (MQ=255)
                                    agGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAttt                                                        <  1:3274824/91‑1 (MQ=255)
                                     ggTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtttt                                                       >  1:2015914/1‑91 (MQ=255)
                                      gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCCTGATTCTg                                                      >  2:1703/1‑91 (MQ=255)
                                      gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                                                      >  1:3615597/1‑91 (MQ=255)
                                       tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                     >  2:741797/1‑91 (MQ=255)
                                         aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                                   <  1:651873/91‑1 (MQ=255)
                                         aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                                   <  1:3377780/91‑1 (MQ=255)
                                         aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                                   >  1:3538738/1‑91 (MQ=255)
                                          cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                                                  >  2:2334927/1‑91 (MQ=255)
                                          cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                                                  >  1:2359238/1‑91 (MQ=255)
                                               tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                             >  2:685618/1‑91 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:601818/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  2:1993772/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            >  2:3724020/1‑91 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:844346/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  2:2870156/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:446948/91‑1 (MQ=255)
                                                  taaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          <  2:3932765/90‑1 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTCTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          >  2:3467508/1‑91 (MQ=255)
                                                  aaaGAAAACAGCAGTAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                                         >  2:3883361/1‑91 (MQ=255)
                                                   aaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                                         <  2:2703650/91‑1 (MQ=255)
                                                    aGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACa                                        >  2:2038606/1‑91 (MQ=255)
                                                      aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                                      >  1:1409266/1‑91 (MQ=255)
                                                       aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                     >  2:2190118/1‑91 (MQ=255)
                                                       aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                     >  1:2606810/1‑91 (MQ=255)
                                                        aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACACttt                                    >  2:3029634/1‑91 (MQ=255)
                                                         aCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtttt                                   <  2:3459678/91‑1 (MQ=255)
                                                          cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATGTTCATGATTTTGTTATATAATCTACAATTTtg                                  >  1:3835247/1‑91 (MQ=255)
                                                           agcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                                 >  1:1213524/1‑91 (MQ=18)
                                                           agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                                 >  1:1675862/1‑91 (MQ=255)
                                                           agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                                 <  1:2760071/91‑1 (MQ=255)
                                                            gcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                                >  2:429798/1‑91 (MQ=21)
                                                            gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                                <  2:2253999/91‑1 (MQ=255)
                                                            gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                                <  1:2033988/91‑1 (MQ=255)
                                                            gcagTAAAATCATTCATAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATCCATGATTTTGTTATATAATCTACAATTTtgtg                                >  1:3374304/1‑91 (MQ=255)
                                                              agTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                              >  2:1852440/1‑91 (MQ=21)
                                                              agTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                              <  2:636321/91‑1 (MQ=255)
                                                              agTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                              >  2:2694996/1‑91 (MQ=255)
                                                                 aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt                           <  2:599959/91‑1 (MQ=255)
                                                                      gaTTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                      >  1:132177/2‑91 (MQ=37)
                                                                      cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                      >  2:2436504/1‑91 (MQ=255)
                                                                      cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                      >  1:3457228/1‑91 (MQ=255)
                                                                       aTTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTc                     >  2:2669357/1‑91 (MQ=255)
                                                                        ttCCTAATTTGGAGTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTATTATATAATCTACAATTTTGTGTCTTTTAAGTCt                    <  2:1161492/91‑1 (MQ=255)
                                                                        ttCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCt                    >  2:1986850/1‑91 (MQ=255)
                                                                        ttCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCt                    >  2:1356603/1‑91 (MQ=255)
                                                                         tCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTTAGTCtt                   >  2:10666/1‑91 (MQ=255)
                                                                           cTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTcc                 >  1:578478/1‑91 (MQ=255)
                                                                            tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg                >  2:1475864/1‑91 (MQ=255)
                                                                                ttGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAt            >  2:1404092/1‑91 (MQ=255)
                                                                                 tGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAtt           <  1:223753/91‑1 (MQ=255)
                                                                                  ggAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAttt          >  1:2660057/1‑91 (MQ=255)
                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        <  1:1901769/91‑1 (MQ=255)
                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        <  2:3227077/91‑1 (MQ=255)
                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        >  1:3405831/1‑91 (MQ=255)
                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        <  2:839674/91‑1 (MQ=255)
                                                                                     atcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCAt       >  2:3446139/1‑91 (MQ=255)
                                                                                     atcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCAt       >  2:567125/1‑91 (MQ=255)
                                                                                      tcatcaTACTGCTGTTTGTTGTTGAGTTATATTCATGATATTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc      <  1:3670962/91‑1 (MQ=255)
                                                                                      tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc      <  1:2306390/91‑1 (MQ=255)
                                                                                          cTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACt  >  1:441125/1‑91 (MQ=255)
                                                                                          |                                                                                          
GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACT  >  USA300TCH1516_ALE/942140‑942320

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: