Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 597,163 A→C 5.6% intergenic (+113/‑104) fusA → / → tuf Elongation factor G/Elongation factor Tu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE597,1630AC 5.6% 648.1 / 24.2 194intergenic (+113/‑104)fusA/tufElongation factor G/Elongation factor Tu
Reads supporting (aligned to +/- strand):  ref base A (83/100);  new base C (4/7);  total (87/107)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 210 total aligned reads displayed.

AGCCTAGGTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTAGAATAG  >  USA300TCH1516_ALE/597072‑597253
                                                                                           |                                                                                          
 gcCTAGGTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTaa                                                                                            <  2:1717133/91‑1 (MQ=255)
  ccTAGGTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTaaa                                                                                           >  2:1909574/1‑91 (MQ=255)
   cTAGGTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTaaaa                                                                                          <  1:2200086/91‑1 (MQ=255)
     aGGCTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAAtt                                                                                        <  1:3128625/91‑1 (MQ=255)
      ggTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTTTACAATTc                                                                                       <  1:2727469/91‑1 (MQ=255)
       gTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTTTACAATTCt                                                                                      <  1:885847/91‑1 (MQ=255)
       gTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCt                                                                                      >  1:1520108/1‑91 (MQ=255)
         tAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTTTACAATTCttt                                                                                    <  1:578198/91‑1 (MQ=255)
         tAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCttt                                                                                    >  1:2687317/1‑91 (MQ=255)
         tAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCttt                                                                                    <  2:1541967/91‑1 (MQ=255)
          aaaaTACAAGGTGAGCGTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTTTACAATTCtttt                                                                                   <  2:3706940/91‑1 (MQ=255)
              tACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCa                                                                               <  1:1444204/91‑1 (MQ=255)
              tACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCa                                                                               >  1:1761635/1‑91 (MQ=255)
                cAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATc                                                                             <  2:3071909/91‑1 (MQ=255)
                cAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATc                                                                             <  2:2471726/91‑1 (MQ=255)
                cAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATc                                                                             >  2:3729253/1‑91 (MQ=255)
                 aaGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCg                                                                            <  2:1097012/91‑1 (MQ=255)
                   ggTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCt                                                                          <  2:1656351/91‑1 (MQ=255)
                   ggTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCt                                                                          <  1:2302323/91‑1 (MQ=255)
                   ggTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCt                                                                          <  1:79937/91‑1 (MQ=255)
                    gTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCta                                                                         >  1:57049/1‑91 (MQ=255)
                    gTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCta                                                                         >  2:3585406/1‑91 (MQ=255)
                     tGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCtat                                                                        >  2:3652224/1‑91 (MQ=255)
                     tGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCtat                                                                        <  2:2064229/91‑1 (MQ=255)
                      gAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTTGTTTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCtata                                                                       >  2:1485621/1‑91 (MQ=255)
                        gCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaa                                                                     <  1:2952934/91‑1 (MQ=255)
                         cTTACATTTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    >  2:2708186/1‑91 (MQ=255)
                         cTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTTGTTTGTATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    >  2:2760996/1‑91 (MQ=255)
                         cTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    <  1:1572515/91‑1 (MQ=255)
                         cTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    <  2:2578133/91‑1 (MQ=255)
                         cTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    >  2:656004/1‑91 (MQ=255)
                         cTTAAATGTAAGATATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTataaat                                                                    >  1:2975877/1‑91 (MQ=255)
                            aaaTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTaaa                                                                                           >  1:1121705/1‑65 (MQ=255)
                            aaaTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTaaa                                                                                           <  2:1121705/65‑1 (MQ=255)
                            aaaTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAaataat                                                                 >  1:400917/1‑91 (MQ=255)
                             aaTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAAtt                                                                >  2:1846483/1‑91 (MQ=255)
                              aTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGTATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAAttt                                                               <  1:2842988/91‑1 (MQ=255)
                              aTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAAttt                                                               >  2:504584/1‑91 (MQ=255)
                              aTGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAAttt                                                               <  1:1107920/91‑1 (MQ=255)
                               tGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTtc                                                              >  2:2863624/1‑91 (MQ=255)
                                gTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTtct                                                             <  2:44603/91‑1 (MQ=255)
                                gTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTtct                                                             <  2:833999/91‑1 (MQ=255)
                                 tAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTtctc                                                            <  2:590122/91‑1 (MQ=255)
                                  aaGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCa                                                           >  2:3514411/1‑91 (MQ=255)
                                   aGCTACCATCTCTATAGTTTGATTTTTCGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCCCat                                                          >  2:1779847/1‑91 (MQ=255)
                                     cTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCatga                                                        >  1:169836/1‑91 (MQ=255)
                                      tATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCatgat                                                       >  1:1344637/1‑91 (MQ=255)
                                         catcTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGt                                                    >  1:1200778/1‑91 (MQ=255)
                                         catcTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGt                                                    <  2:1607535/91‑1 (MQ=255)
                                             tttATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTgaga                                                >  2:2014065/1‑91 (MQ=255)
                                             tttATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTgaga                                                >  1:537786/1‑91 (MQ=255)
                                             tttATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTgaga                                                <  1:786223/91‑1 (MQ=255)
                                             tttATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTgaga                                                <  2:2951460/91‑1 (MQ=255)
                                              ttATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGaa                                               <  1:245489/91‑1 (MQ=255)
                                               tataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGaaa                                              <  2:897615/91‑1 (MQ=255)
                                                ataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAAc                                             >  2:1301249/1‑91 (MQ=255)
                                                ataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAAc                                             <  2:725505/91‑1 (MQ=255)
                                                ataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAAc                                             <  1:3283644/91‑1 (MQ=255)
                                                ataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAAc                                             <  1:3514944/91‑1 (MQ=255)
                                                ataGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTGTTTGCATCGCTATAAGTAGTTTCTCATGATGGTGAGAAAc                                             <  2:1367717/91‑1 (MQ=255)
                                                 tcgTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTAATGATGGTGAGAAACt                                            >  2:3796727/3‑91 (MQ=255)
                                                 taGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACt                                            >  1:2529724/1‑91 (MQ=255)
                                                   gTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTAt                                          >  1:3540414/1‑91 (MQ=255)
                                                      tGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCAt                                       >  2:2063065/1‑91 (MQ=255)
                                                        aTTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATga                                     >  1:1161382/1‑91 (MQ=255)
                                                        aTTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATga                                     <  1:3638053/91‑1 (MQ=255)
                                                         ttttttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATgag                                    >  1:692436/1‑91 (MQ=255)
                                                           ttttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATgagag                                  <  2:402436/91‑1 (MQ=255)
                                                            tttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATgagaga                                 <  1:1024781/91‑1 (MQ=255)
                                                             ttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGAt                                <  2:1200778/91‑1 (MQ=255)
                                                             ttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGAt                                <  1:1805482/91‑1 (MQ=255)
                                                             ttGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGAt                                >  2:3281193/1‑91 (MQ=255)
                                                               ggggTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATaa                              >  1:989519/1‑91 (MQ=255)
                                                               ggggTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATaa                              >  1:3483492/1‑91 (MQ=255)
                                                                  gTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAAtt                           >  2:3404988/1‑91 (MQ=255)
                                                                  gTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAAtt                           >  1:3857418/1‑91 (MQ=255)
                                                                    gAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTa                         >  1:1863997/1‑91 (MQ=255)
                                                                    gAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTa                         <  2:2806242/91‑1 (MQ=255)
                                                                    gAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTa                         <  1:3047139/91‑1 (MQ=255)
                                                                      aTGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTaaa                       <  1:2137878/91‑1 (MQ=255)
                                                                      aTGCATTATAAAAGAAGTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTaaa                       <  2:2751502/91‑1 (MQ=255)
                                                                        gCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAata                     <  2:3525830/91‑1 (MQ=255)
                                                                         cATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAatat                    <  1:2120153/91‑1 (MQ=255)
                                                                           ttATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAAtatta                  >  2:3700911/1‑91 (MQ=255)
                                                                           ttATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAAtatta                  >  1:2089973/1‑91 (MQ=255)
                                                                           ttATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAAtatta                  >  1:1975350/1‑91 (MQ=255)
                                                                           ttATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAAtatta                  <  2:710352/91‑1 (MQ=255)
                                                                           ttATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAAtatta                  <  1:1435593/91‑1 (MQ=255)
                                                                               aaaaGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTAtttt              >  2:15543/1‑91 (MQ=255)
                                                                               aaaaGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTAtttt              <  2:2847372/91‑1 (MQ=255)
                                                                                aaaGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTAttttt             <  1:280766/91‑1 (MQ=255)
                                                                                  aGAATTTTACAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTaa           >  1:642687/1‑91 (MQ=255)
                                                                                  aGAATTTTACAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTaa           >  2:2734988/1‑91 (MQ=255)
                                                                                  aGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTaa           <  1:433876/91‑1 (MQ=255)
                                                                                   gAATTTTACAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAAt          >  1:3837086/1‑91 (MQ=255)
                                                                                    aaTTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAAtt         <  1:720637/91‑1 (MQ=255)
                                                                                     aTTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTa        >  1:3376078/1‑91 (MQ=255)
                                                                                      ttGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTAg       <  1:3787325/91‑1 (MQ=255)
                                                                                          aaaaTTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTAGAATa   >  2:1911220/1‑91 (MQ=255)
                                                                                          aaaaTTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTAGAATa   >  2:2158396/1‑91 (MQ=255)
                                                                                           |                                                                                          
AGCCTAGGTTAAAATACAAGGTGAGCTTAAATGTAAGCTATCATCTTTATAGTTTGATTTTTTGGGGTGAATGCATTATAAAAGAATTGTAAAATTCTTTTTGCATCGCTATAAATAATTTCTCATGATGGTGAGAAACTATCATGAGAGATAAATTTAAATATTATTTTTAATTAGAATAG  >  USA300TCH1516_ALE/597072‑597253

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: