Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,210 C→G 8.8% intergenic (‑181/+68) 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,2100CG 8.8% 748.3 / 45.5 216intergenic (‑181/+68)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base C (98/99);  new base G (7/12);  total (105/111)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 223 total aligned reads displayed.

AAGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGT  >  USA300TCH1516_ALE/942120‑942300
                                                                                          |                                                                                          
aaGAGATATAGGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAGTc                                                                                            <  2:1373978/91‑1 (MQ=255)
  gagaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                          <  1:2100922/91‑1 (MQ=255)
  gagaTATATGTAGATAGGGCGAATCTGAAAAAAAAAAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAt                                                                                          <  2:2904325/91‑1 (MQ=255)
   agaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTTAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAtt                                                                                         >  1:3145427/1‑91 (MQ=255)
   agaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCAtt                                                                                         >  1:2979958/1‑91 (MQ=255)
    gaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTc                                                                                        <  1:844275/91‑1 (MQ=255)
    gaTATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTc                                                                                        <  2:2107578/91‑1 (MQ=255)
     atatatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATTTCAAAAGAAAACAGCAGTAAAATCATTcc                                                                                       <  2:1546332/91‑1 (MQ=255)
      tatatGTAGGTAGGGCGAATCAGTAGTCAAAGAAAAAATCAATGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCt                                                                                      <  2:2087675/91‑1 (MQ=255)
      tatatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCt                                                                                      >  2:3350531/1‑91 (MQ=255)
       atatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTa                                                                                     >  2:115729/1‑91 (MQ=255)
        tatGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTaa                                                                                    <  2:142608/91‑1 (MQ=255)
          tGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAAtt                                                                                  <  1:3213326/91‑1 (MQ=255)
          tGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAAACATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCAAAtt                                                                                  <  2:281299/91‑1 (MQ=255)
           gTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAAttt                                                                                 <  2:939194/91‑1 (MQ=255)
             aGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTgg                                                                               >  2:714224/1‑91 (MQ=255)
                tAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAat                                                                            >  2:2160891/1‑91 (MQ=255)
                tAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAat                                                                            <  1:3515562/91‑1 (MQ=255)
                  gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                          <  2:1035788/91‑1 (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)
                        aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCTTTAAGATGATTCCAAATTAGGAATGATTTTAc                                                                    <  2:2272484/91‑1 (MQ=21)
                        aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                    <  1:776236/91‑1 (MQ=255)
                          cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActg                                                                  <  1:1847696/91‑1 (MQ=18)
                          cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActg                                                                  >  2:3396971/1‑91 (MQ=255)
                          cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTAActg                                                                  >  1:2939607/1‑91 (MQ=255)
                           tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                 <  1:1294841/91‑1 (MQ=18)
                           tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                 <  2:2018719/91‑1 (MQ=18)
                           tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                 <  1:427434/91‑1 (MQ=18)
                           tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                 <  2:2831688/91‑1 (MQ=18)
                            gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                >  1:1830558/1‑91 (MQ=255)
                             tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                               >  1:624472/1‑91 (MQ=255)
                             tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                               >  2:3584328/1‑91 (MQ=255)
                             tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                               >  2:2599297/1‑91 (MQ=255)
                               gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                             >  2:597386/1‑91 (MQ=255)
                                tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTACTGCtgttt                                                            <  2:2214811/91‑1 (MQ=18)
                                tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                            <  2:2309792/91‑1 (MQ=255)
                                  aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                          <  2:3368878/91‑1 (MQ=255)
                                  aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                          <  1:2029837/91‑1 (MQ=255)
                                    aGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttg                                                        <  1:1023401/91‑1 (MQ=255)
                                     gAAAAAATCATTGAAAAAGAGGTCACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                       >  2:77278/1‑91 (MQ=255)
                                     gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                       <  1:3413883/91‑1 (MQ=255)
                                     gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                       >  2:3820928/1‑91 (MQ=255)
                                        aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                    <  1:965232/91‑1 (MQ=255)
                                        aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                    >  2:1245220/1‑91 (MQ=255)
                                        aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                    >  2:3089040/1‑91 (MQ=255)
                                           aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                 >  1:2699011/1‑91 (MQ=255)
                                           aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                 >  1:3278777/1‑91 (MQ=255)
                                            tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                >  1:1563311/1‑91 (MQ=255)
                                            tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                >  1:2443777/1‑91 (MQ=255)
                                             cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                               >  2:1313843/1‑91 (MQ=255)
                                             cATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTtat                                               >  1:2938071/1‑91 (MQ=255)
                                                 gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                           <  1:1647655/91‑1 (MQ=255)
                                                 gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                           <  1:2131774/91‑1 (MQ=255)
                                                     aaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGa                                       >  1:1686623/1‑91 (MQ=255)
                                                       gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                     >  2:219131/1‑91 (MQ=255)
                                                       gagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtt                                     >  1:635358/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                                 >  2:741797/1‑91 (MQ=255)
                                                            aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                >  1:3540336/1‑91 (MQ=255)
                                                             aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                               <  1:3377569/91‑1 (MQ=255)
                                                             aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                               <  1:651873/91‑1 (MQ=255)
                                                              cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                              >  2:2334927/1‑91 (MQ=255)
                                                              cAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtat                              >  1:2359238/1‑91 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  2:1993772/91‑1 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  2:2870156/91‑1 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        >  2:3724020/1‑91 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  2:2694668/91‑1 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  2:1995149/91‑1 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  1:446948/91‑1 (MQ=255)
                                                                    cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                        <  1:601818/91‑1 (MQ=255)
                                                                     aaaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCt                       <  2:1813999/91‑1 (MQ=255)
                                                                      taaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                      <  2:3932765/90‑1 (MQ=255)
                                                                      aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTCTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                      >  2:3467508/1‑91 (MQ=255)
                                                                      aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGATGTTTGTTGTTGATTTATATTCATGATTTTGTTATATATTCTa                      >  2:63907/1‑91 (MQ=255)
                                                                      aaaGAAAACAGCAGTAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                     >  2:3883361/1‑91 (MQ=255)
                                                                       aaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                     <  2:2703566/91‑1 (MQ=255)
                                                                        aGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACa                    >  2:2038606/1‑91 (MQ=255)
                                                                          aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                  >  1:1409266/1‑91 (MQ=255)
                                                                           aaaCAGCAGTAAGATGATTCCAAATTAGGCATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTGATATAAGCTACAAtt                 >  2:3436236/1‑91 (MQ=17)
                                                                           aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                 >  1:2606810/1‑91 (MQ=255)
                                                                            aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                >  1:566119/1‑91 (MQ=255)
                                                                             aCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtttt               <  2:3459678/91‑1 (MQ=255)
                                                                              cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg              >  1:3073034/1‑91 (MQ=255)
                                                                              cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg              >  2:317030/1‑91 (MQ=255)
                                                                                gcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg            >  2:429798/1‑91 (MQ=21)
                                                                                gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg            >  1:1188598/1‑91 (MQ=255)
                                                                                gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg            >  1:334443/1‑91 (MQ=255)
                                                                                 cagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt           >  2:2403148/1‑91 (MQ=255)
                                                                                   gTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCt         >  1:3769195/1‑91 (MQ=255)
                                                                                      aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt      <  1:3718355/91‑1 (MQ=255)
                                                                                       aaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTa     >  2:943398/1‑91 (MQ=255)
                                                                                        aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTaa    <  1:2500399/91‑1 (MQ=255)
                                                                                         tCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAg   >  2:86079/1‑91 (MQ=255)
                                                                                          cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt  >  2:2436504/1‑91 (MQ=255)
                                                                                          cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt  >  1:3457228/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAGAGATATATGTAGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGT  >  USA300TCH1516_ALE/942120‑942300

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