Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,227 C→G 8.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 8.2% 335.1 / 17.1 97intergenic (‑198/+51)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base C (39/50);  new base G (3/5);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 104 total aligned reads displayed.

AGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATC  >  USA300TCH1516_ALE/942137‑942316
                                                                                          |                                                                                         
aGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatc                                                                                           >  1:3601013/1‑91 (MQ=255)
 gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                          >  1:2796726/1‑91 (MQ=255)
 gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                          >  1:3021717/1‑91 (MQ=255)
   gCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatcatc                                                                                        >  1:104701/1‑91 (MQ=255)
     gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                                      >  1:533918/1‑91 (MQ=255)
     gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                                      <  1:1489511/91‑1 (MQ=255)
       aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAATAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTAc                                                                                    <  1:3598812/91‑1 (MQ=18)
        tCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAct                                                                                   >  1:1225609/1‑91 (MQ=255)
          tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                                 <  1:1700773/91‑1 (MQ=18)
          tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgc                                                                                 >  1:1985952/1‑91 (MQ=255)
           gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                <  1:4827/91‑1 (MQ=255)
           gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                >  1:2490894/1‑91 (MQ=255)
            tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgctg                                                                               <  1:753994/91‑1 (MQ=18)
            tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                               >  1:1266410/1‑91 (MQ=255)
               tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                            <  1:3919181/91‑1 (MQ=255)
               tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                            <  1:3919152/91‑1 (MQ=255)
               tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                            >  1:875880/1‑91 (MQ=255)
               tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                            <  1:2032707/91‑1 (MQ=255)
                cAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttg                                                                           <  1:350009/91‑1 (MQ=255)
                 aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                                          <  1:535250/91‑1 (MQ=255)
                    gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                                       <  1:3657406/91‑1 (MQ=255)
                    gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                                       <  1:1607426/91‑1 (MQ=255)
                       aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTTCTAATTTGGAATCATCTTACTGCTGTTTGATGTTGa                                                                    >  1:458101/1‑91 (MQ=255)
                       aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                    <  1:1060791/91‑1 (MQ=255)
                       aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                    <  1:253345/91‑1 (MQ=255)
                        aaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAt                                                                   >  1:1684026/1‑91 (MQ=255)
                        aaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAt                                                                   >  1:295199/1‑91 (MQ=255)
                         aaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAtt                                                                  >  1:16992/1‑91 (MQ=255)
                         aaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAtt                                                                  >  1:16554/1‑91 (MQ=255)
                          aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                 >  1:2822081/1‑91 (MQ=255)
                          aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                 >  1:3301526/1‑91 (MQ=255)
                           tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                <  1:1870451/91‑1 (MQ=255)
                               tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                            <  1:2916775/91‑1 (MQ=255)
                               tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                            <  1:3521870/91‑1 (MQ=255)
                                gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                           <  1:2779837/91‑1 (MQ=255)
                                gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                           <  1:2017584/91‑1 (MQ=255)
                                  aaaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCAt                                                         <  1:2588816/91‑1 (MQ=255)
                                   aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                                        <  1:2199202/91‑1 (MQ=255)
                                     agagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAt                                                      <  1:2337100/91‑1 (MQ=255)
                                       agGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAttt                                                    >  1:3705692/1‑91 (MQ=255)
                                          tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                 >  1:2498002/1‑91 (MQ=255)
                                           aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                                <  1:2573813/91‑1 (MQ=255)
                                           aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                                >  1:3721588/1‑91 (MQ=255)
                                            aCAATGTCAAAATAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                               >  1:3189429/1‑91 (MQ=255)
                                            aCAATGTCAAAATAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                               >  1:3189444/1‑91 (MQ=255)
                                            aCAATGTCAAAAGAAATCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                               <  1:4120684/91‑1 (MQ=255)
                                            aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                               >  1:506879/1‑91 (MQ=255)
                                            aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                               >  1:507051/1‑91 (MQ=255)
                                                  tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                         >  1:2553355/1‑91 (MQ=255)
                                                  tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                         >  1:3156022/1‑91 (MQ=255)
                                                   cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                        <  1:1307216/91‑1 (MQ=255)
                                                   cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                        <  1:3058748/91‑1 (MQ=255)
                                                   cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                        >  1:1488563/1‑91 (MQ=255)
                                                   cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                        <  1:2702648/91‑1 (MQ=255)
                                                   cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                        <  1:1181091/91‑1 (MQ=255)
                                                     aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                      <  1:3445736/91‑1 (MQ=255)
                                                      aaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                                     >  1:2800365/1‑91 (MQ=255)
                                                       aGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACa                                    >  1:3711113/1‑91 (MQ=255)
                                                        gAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACaa                                   <  1:2708133/91‑1 (MQ=255)
                                                         aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                                  <  1:563914/91‑1 (MQ=255)
                                                          aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                 >  1:2444220/1‑91 (MQ=255)
                                                          aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGGTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                 <  1:2395792/91‑1 (MQ=255)
                                                           aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                <  1:635570/91‑1 (MQ=255)
                                                           aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                >  1:1968583/1‑91 (MQ=255)
                                                           aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                <  1:636692/91‑1 (MQ=255)
                                                           aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTGGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                >  1:3297950/1‑91 (MQ=255)
                                                             cagcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                              >  1:517423/1‑91 (MQ=18)
                                                             cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                              <  1:1905627/91‑1 (MQ=255)
                                                             cagcagTAAAATAATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                              >  1:3178549/1‑91 (MQ=255)
                                                              agcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                             >  1:2833859/1‑91 (MQ=18)
                                                               gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                            >  1:443696/1‑91 (MQ=255)
                                                                 agTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                          >  1:321213/1‑91 (MQ=18)
                                                                  gTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCt                         >  1:3069051/1‑91 (MQ=255)
                                                                   tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                        <  1:4009488/91‑1 (MQ=255)
                                                                   tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                        >  1:103929/1‑91 (MQ=255)
                                                                   tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                        <  1:4009350/91‑1 (MQ=255)
                                                                   tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                        <  1:3966882/91‑1 (MQ=255)
                                                                   caaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                        <  1:4009677/90‑1 (MQ=255)
                                                                    aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt                       <  1:304073/91‑1 (MQ=255)
                                                                    aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt                       <  1:359488/91‑1 (MQ=255)
                                                                     aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                      <  1:1246288/91‑1 (MQ=255)
                                                                     aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                      >  1:3641303/1‑91 (MQ=255)
                                                                     aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                      <  1:1246102/91‑1 (MQ=255)
                                                                       aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTat                    >  1:3607708/1‑90 (MQ=255)
                                                                       aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTaa                    <  1:143984/91‑1 (MQ=255)
                                                                        tCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAg                   >  1:2025626/1‑91 (MQ=255)
                                                                        tCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAg                   >  1:2025657/1‑91 (MQ=255)
                                                                         cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                  <  1:2198293/91‑1 (MQ=255)
                                                                          aTTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTc                 <  1:2378101/91‑1 (MQ=255)
                                                                            tCCTAATTTGGAATCATCTTACTGCTGTTTGCTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTTTCTTTTAAGTCtt               >  1:3434241/1‑91 (MQ=255)
                                                                              cTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTcc             >  1:2582358/1‑91 (MQ=255)
                                                                              cTAATTTGGAATCAGCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCGACAATTTTGTGTCTTTTAAGTCTTcc             <  1:3742117/91‑1 (MQ=255)
                                                                               tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg            <  1:1462518/91‑1 (MQ=255)
                                                                               tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg            <  1:1990169/91‑1 (MQ=255)
                                                                                 aTTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGaa          >  1:1877810/1‑91 (MQ=255)
                                                                                     ggAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAttt      <  1:1545232/91‑1 (MQ=255)
                                                                                       aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa    <  1:2222437/91‑1 (MQ=255)
                                                                                       aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa    <  1:382065/91‑1 (MQ=255)
                                                                                         tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc  <  1:2319523/91‑1 (MQ=255)
                                                                                         tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc  <  1:1260196/91‑1 (MQ=255)
                                                                                          |                                                                                         
AGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATC  >  USA300TCH1516_ALE/942137‑942316

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