Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 942,230 C→T 8.0% 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 8.0% 722.9 / 33.5 226intergenic (‑201/+48)ptlE/mnhG1Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand):  ref base C (110/97);  new base T (9/9);  total (120/106)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 238 total aligned reads displayed.

GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACT  >  USA300TCH1516_ALE/942140‑942320
                                                                                          |                                                                                          
gCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatcatc                                                                                            >  2:3790911/1‑91 (MQ=255)
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           <  1:2562029/91‑1 (MQ=255)
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           >  1:2340185/1‑91 (MQ=255)
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           >  1:1579280/1‑91 (MQ=255)
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           <  2:1899993/91‑1 (MQ=255)
 cGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCt                                                                                           >  1:3871184/1‑91 (MQ=255)
   aaTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTa                                                                                         <  2:731823/91‑1 (MQ=255)
    gtCTGTAGTCAAAGAAGAAATCATTGAAAAAGAGGTAGCAATGTCAAAAGAAGACGGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        <  2:3847855/90‑1 (MQ=255)
    aTCTGTAGTCAAAGTAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACCGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        <  1:3624709/91‑1 (MQ=255)
    aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        <  1:3576815/91‑1 (MQ=255)
    aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        >  1:720919/1‑91 (MQ=255)
    aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAc                                                                                        >  2:349952/1‑91 (MQ=255)
      cTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActg                                                                                      <  2:1494860/91‑1 (MQ=18)
       tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgc                                                                                     >  2:292071/1‑91 (MQ=255)
       tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgc                                                                                     <  1:1645221/91‑1 (MQ=255)
        gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgct                                                                                    <  1:1114950/91‑1 (MQ=18)
        gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgct                                                                                    <  1:2623985/91‑1 (MQ=18)
        gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                    >  1:2156241/1‑91 (MQ=255)
        gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                                    >  1:3862837/1‑91 (MQ=255)
         tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                   <  1:2427497/91‑1 (MQ=255)
         tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                                   <  1:2483460/91‑1 (MQ=255)
          aGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgt                                                                                  >  1:2867346/1‑91 (MQ=255)
           gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                                 <  2:1179966/91‑1 (MQ=255)
           gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                                 >  2:3220414/1‑91 (MQ=255)
           gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                                 <  1:371054/91‑1 (MQ=255)
           gTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtt                                                                                 <  1:2908869/91‑1 (MQ=255)
            tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                                <  1:1379149/91‑1 (MQ=255)
              aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                                              <  1:2197485/91‑1 (MQ=255)
              aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                                              <  1:3527840/91‑1 (MQ=255)
                aGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttttcac                                                                            >  1:2129535/1‑87 (MQ=255)
                 gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                                           >  1:881559/1‑91 (MQ=255)
                 gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                                           <  2:3697506/91‑1 (MQ=255)
                    aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAGCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                        <  1:711648/91‑1 (MQ=255)
                    aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                                        >  1:172040/1‑91 (MQ=255)
                        tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                                    >  2:2822513/1‑91 (MQ=255)
                          aTTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgttg                                                                         >  1:3654715/1‑84 (MQ=255)
                           ttGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTatat                                                                 <  1:3363024/91‑1 (MQ=255)
                            tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                                <  1:3589078/91‑1 (MQ=255)
                            tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                                <  2:2182861/91‑1 (MQ=255)
                            tGAAAAAGAGGGAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                                <  1:3512437/91‑1 (MQ=255)
                              aaaaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCa                                                              >  1:2834081/1‑91 (MQ=255)
                              aaaaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCa                                                              >  2:1335388/1‑91 (MQ=255)
                                aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                                            <  1:2269803/91‑1 (MQ=255)
                                 aaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGa                                                           <  2:2605317/91‑1 (MQ=255)
                                     ggTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAtttt                                                       <  2:2232835/91‑1 (MQ=255)
                                      gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                                                      >  1:3781681/1‑91 (MQ=255)
                                      gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                                                      >  1:1144356/1‑91 (MQ=255)
                                      gTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTg                                                      >  2:1055118/1‑91 (MQ=255)
                                       tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                     >  2:115813/1‑91 (MQ=255)
                                       tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                                     >  2:1127521/1‑91 (MQ=255)
                                             tGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTtatata                                               >  2:1351140/1‑91 (MQ=255)
                                               tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                             >  2:2108516/1‑91 (MQ=255)
                                               tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                             >  1:98548/1‑91 (MQ=255)
                                               tCAAAAGAAAACAGCAGTAAAATCATTCCTAACTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                             >  2:3813426/1‑91 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:2934750/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:2241698/91‑1 (MQ=255)
                                                cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                            <  1:1692758/91‑1 (MQ=255)
                                                 aaaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCt                                           <  2:144589/91‑1 (MQ=255)
                                                 aaaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCt                                           <  1:798510/91‑1 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          <  2:1570681/91‑1 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          <  2:29745/91‑1 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          <  2:3153164/91‑1 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          >  2:796099/1‑91 (MQ=255)
                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                          <  2:888086/91‑1 (MQ=255)
                                                    aGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACa                                        >  2:862959/1‑91 (MQ=255)
                                                     gAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACaa                                       >  1:2767407/1‑91 (MQ=255)
                                                       aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                                     <  1:2975964/91‑1 (MQ=255)
                                                         aCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtttt                                   <  1:1067813/91‑1 (MQ=255)
                                                          cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACCATTTtg                                  >  1:3027806/1‑91 (MQ=255)
                                                          cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                                  >  1:1456654/1‑91 (MQ=255)
                                                           agcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                                 <  1:3162342/91‑1 (MQ=255)
                                                            gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                                <  2:3272655/91‑1 (MQ=255)
                                                            gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                                <  2:3444113/91‑1 (MQ=255)
                                                             cagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                               >  1:2039867/1‑91 (MQ=18)
                                                             cagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                               >  2:2069746/1‑91 (MQ=255)
                                                             cagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                               <  1:1732988/91‑1 (MQ=255)
                                                             cagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                               >  2:2538286/1‑91 (MQ=255)
                                                             cagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTgtgt                               >  1:3846455/1‑91 (MQ=255)
                                                              agTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                              >  2:1797794/1‑91 (MQ=21)
                                                                 aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt                           >  2:196545/1‑91 (MQ=255)
                                                                  aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                          >  1:1542368/1‑91 (MQ=255)
                                                                  aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                          >  2:2860943/1‑91 (MQ=255)
                                                                  aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                          <  1:2108516/91‑1 (MQ=255)
                                                                    aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTaa                        <  2:453926/91‑1 (MQ=255)
                                                                    aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCGTTTaa                        >  2:3577159/1‑91 (MQ=255)
                                                                     tCATTCCTAATTTCGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATGTTGTGTCTTTTAAg                       <  2:1718416/91‑1 (MQ=255)
                                                                      cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTTAGt                      >  2:1072331/1‑91 (MQ=255)
                                                                      cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt                      >  1:1476335/1‑91 (MQ=255)
                                                                          ccTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTc                  <  1:1431722/91‑1 (MQ=255)
                                                                             aaTTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGa               >  1:3798221/1‑91 (MQ=255)
                                                                                ttGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAt            >  2:1084228/1‑91 (MQ=255)
                                                                                ttGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAt            >  2:1692174/1‑91 (MQ=255)
                                                                                 tGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAtt           >  1:310884/1‑91 (MQ=255)
                                                                                  ggAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAttt          <  2:3105966/91‑1 (MQ=255)
                                                                                    gatcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        <  2:2355164/90‑1 (MQ=255)
                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa        >  2:1828597/1‑91 (MQ=255)
                                                                                     atcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCAt       >  1:1059230/1‑91 (MQ=255)
                                                                                     atcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCAt       >  1:1870374/1‑91 (MQ=255)
                                                                                       catcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCg     <  2:526787/91‑1 (MQ=255)
                                                                                          cTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACt  >  1:2722836/1‑91 (MQ=255)
                                                                                          |                                                                                          
GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGACT  >  USA300TCH1516_ALE/942140‑942320

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: