Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,207 A→G 8.9% intergenic (‑178/+71) 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,2070AG 8.9% 654.2 / 37.6 216intergenic (‑178/+71)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base A (98/97);  new base G (7/13);  total (106/110)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

TTAAAGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTA  >  USA300TCH1516_ALE/942117‑942297
                                                                                          |                                                                                          
ttaaaGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTaaa                                                                                            <  2:2558621/91‑1 (MQ=255)
  aaaGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAAt                                                                                          >  1:471430/1‑91 (MQ=255)
   aaGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATc                                                                                         <  1:430074/91‑1 (MQ=255)
   aaGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATc                                                                                         >  1:252201/1‑91 (MQ=255)
   aaGAGATATAGGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAGTc                                                                                         <  2:1373978/91‑1 (MQ=255)
    agagATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCa                                                                                        >  2:2841354/1‑91 (MQ=255)
     gagaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                       <  2:231934/91‑1 (MQ=255)
     gagaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                       <  2:2091508/91‑1 (MQ=255)
     gagaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                       <  1:2100922/91‑1 (MQ=255)
     gagaTATATGTAGATAGGGCGAATCTGAAAAAAAAAAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                       <  2:2904076/91‑1 (MQ=255)
      agaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTTAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAtt                                                                                      >  1:3145427/1‑91 (MQ=255)
        atatatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATTTCAAAAGAAAACAGCAGTAAAATCATTcc                                                                                    <  2:1546332/91‑1 (MQ=255)
        atatatGTAGAGAGGGCGAATCAAAAGAAAAAAAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTcc                                                                                    <  2:83945/91‑1 (MQ=255)
         tatatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCt                                                                                   >  2:3351806/1‑91 (MQ=255)
          atatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTa                                                                                  >  2:115729/1‑91 (MQ=255)
           tatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTaa                                                                                 <  2:142608/91‑1 (MQ=255)
             tGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTCAtt                                                                               >  2:762636/1‑91 (MQ=255)
             tGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAAACATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCAAAtt                                                                               <  2:281299/91‑1 (MQ=255)
              gTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAAttt                                                                              <  2:939194/91‑1 (MQ=255)
                   tAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAat                                                                         <  1:3515562/91‑1 (MQ=255)
                     gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                       <  2:1035788/91‑1 (MQ=255)
                     gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                       >  2:2083041/1‑91 (MQ=255)
                       gCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatcatc                                                                     >  1:935335/1‑91 (MQ=255)
                        ggAATCTGTAGTCAAAGAAAAAAACATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                    <  2:2711448/90‑1 (MQ=255)
                        cGAATCTGTAGTCAAAGAAAAAATCAAAGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                    <  2:918376/91‑1 (MQ=255)
                         gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                   <  1:2859545/91‑1 (MQ=255)
                           aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                 <  1:776236/91‑1 (MQ=255)
                             cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActg                                                               <  1:1847696/91‑1 (MQ=18)
                             cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                               <  1:3722239/91‑1 (MQ=255)
                             cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                               >  1:2939744/1‑91 (MQ=255)
                             cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                               >  2:3396971/1‑91 (MQ=255)
                              tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                              <  2:2831599/91‑1 (MQ=18)
                              tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                              <  1:1294841/91‑1 (MQ=18)
                               gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                             >  1:1830558/1‑91 (MQ=255)
                                tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgctg                                                            <  2:200649/91‑1 (MQ=18)
                                  gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                          <  2:3799528/91‑1 (MQ=255)
                                  gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                          >  2:597386/1‑91 (MQ=255)
                                   tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTACTGCtgttt                                                         <  2:2214746/91‑1 (MQ=17)
                                     aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                       <  1:2029837/91‑1 (MQ=255)
                                       aGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttg                                                     <  1:1023401/91‑1 (MQ=255)
                                        gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                    <  1:3413883/91‑1 (MQ=255)
                                        gAAAAAATCATAAAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                    <  2:2021903/91‑1 (MQ=255)
                                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                 >  2:3089040/1‑91 (MQ=255)
                                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                 <  1:965232/91‑1 (MQ=255)
                                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                 >  2:1245220/1‑91 (MQ=255)
                                               tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                             >  1:1563311/1‑91 (MQ=255)
                                               tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                             >  1:3656123/1‑91 (MQ=255)
                                                cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                            >  1:2938071/1‑91 (MQ=255)
                                                cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                            >  2:1313843/1‑91 (MQ=255)
                                                 aTTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtata                                           >  1:3323985/1‑91 (MQ=255)
                                                  ttGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTatat                                          >  1:3558320/1‑91 (MQ=255)
                                                    gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                        <  1:2131774/91‑1 (MQ=255)
                                                    gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                        <  1:1647655/91‑1 (MQ=255)
                                                      aaaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCAt                                      <  1:2709779/91‑1 (MQ=255)
                                                       aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                     >  2:2376218/1‑91 (MQ=255)
                                                        aaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGa                                    >  1:1686623/1‑91 (MQ=255)
                                                          gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                  >  1:635358/1‑91 (MQ=255)
                                                           agGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAttt                                 <  1:3274824/91‑1 (MQ=255)
                                                             gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCCTGATTCTg                               >  2:1703/1‑91 (MQ=255)
                                                             gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                               >  2:2114607/1‑91 (MQ=255)
                                                                aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                            <  1:651873/91‑1 (MQ=255)
                                                                aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                            >  1:3538738/1‑91 (MQ=255)
                                                                aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                            <  1:3377569/91‑1 (MQ=255)
                                                                 cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                           >  2:3463001/1‑91 (MQ=255)
                                                                 cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                           >  1:2359238/1‑91 (MQ=255)
                                                                   aTGTCAAAAGAAAGCCGCAGTAACATCATTCATCATTTGGAATCATCCTACTGCTGTTTGATGTTGATTTATATTCATGATCTTGTtatat                         >  2:405978/1‑91 (MQ=255)
                                                                   aTGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtatat                         <  2:3845815/91‑1 (MQ=255)
                                                                      tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                      >  2:685618/1‑91 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  1:601818/91‑1 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     >  2:3724020/1‑91 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  1:630337/91‑1 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  1:844346/91‑1 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  2:1995149/91‑1 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  2:2694291/91‑1 (MQ=255)
                                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                     <  2:1993772/91‑1 (MQ=255)
                                                                        aaaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCt                    <  2:1813999/91‑1 (MQ=255)
                                                                         aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATcaa                   <  2:733744/91‑3 (MQ=255)
                                                                         aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                   >  2:3886480/1‑91 (MQ=255)
                                                                             aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt               >  1:1409266/1‑91 (MQ=255)
                                                                             aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt               <  1:3005678/91‑1 (MQ=255)
                                                                             aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt               <  2:242772/91‑1 (MQ=255)
                                                                              aaaCAGCAGTAAGATGATTCCAAATTAGGCATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTGATATAAGCTACAAtt              >  2:3436236/1‑91 (MQ=17)
                                                                              aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt              >  1:2606973/1‑91 (MQ=255)
                                                                              aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt              >  1:1940634/1‑91 (MQ=255)
                                                                              aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt              >  2:2190118/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)
                                                                                 cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATGTTCATGATTTTGTTATATAATCTACAATTTtg           >  1:3835247/1‑91 (MQ=255)
                                                                                 cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg           >  2:317030/1‑91 (MQ=255)
                                                                                 cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg           >  1:3073034/1‑91 (MQ=255)
                                                                                  agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt          <  2:3580040/91‑1 (MQ=255)
                                                                                   gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg         >  1:1188598/1‑91 (MQ=255)
                                                                                   gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg         <  2:2253999/91‑1 (MQ=255)
                                                                                   gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg         >  1:334473/1‑91 (MQ=255)
                                                                                    cagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt        >  2:2279246/1‑91 (MQ=21)
                                                                                    cagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt        >  1:243952/1‑91 (MQ=18)
                                                                                     agTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc       <  2:636321/91‑1 (MQ=255)
                                                                                     agTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc       >  2:2694996/1‑91 (MQ=255)
                                                                                          aaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTa  >  2:943398/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTAAAGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTA  >  USA300TCH1516_ALE/942117‑942297

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: