Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 299,460 A→C 5.6% intergenic (+2/‑274) tarF_1 → / → tarI Teichoic acid glycerol‑phosphate transferase/Ribitol‑5‑phosphate cytidylyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE299,4600AC 5.6% 779.7 / 33.6 235intergenic (+2/‑274)tarF_1/tarITeichoic acid glycerol‑phosphate transferase/Ribitol‑5‑phosphate cytidylyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (95/126);  new base C (3/10);  total (99/136)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 244 total aligned reads displayed.

TTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTATG  >  USA300TCH1516_ALE/299370‑299550
                                                                                          |                                                                                          
 tGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATaa                                                                                           <  2:3365382/91‑1 (MQ=255)
  gAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAAc                                                                                          >  1:182255/1‑91 (MQ=255)
   aaGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACt                                                                                         <  2:1284393/91‑1 (MQ=255)
    aGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACtt                                                                                        <  1:73322/91‑1 (MQ=255)
    aGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACtt                                                                                        >  1:2372024/1‑91 (MQ=255)
      cACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTATGTTATCACTTaa                                                                                      <  1:2979335/91‑1 (MQ=255)
      cACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTaa                                                                                      <  2:631506/91‑1 (MQ=255)
      cACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTTATCACTTaa                                                                                      <  2:2935695/91‑1 (MQ=255)
        cGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTaaaa                                                                                    >  2:778140/1‑91 (MQ=255)
          aaCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAAc                                                                                  <  2:862996/91‑1 (MQ=255)
          aaCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAAc                                                                                  <  1:2512327/91‑1 (MQ=255)
           aCGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACa                                                                                 >  1:483802/1‑91 (MQ=255)
            cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                                >  1:3805426/1‑91 (MQ=255)
            cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                                >  2:2600723/1‑91 (MQ=255)
            cGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACaa                                                                                >  1:928943/1‑91 (MQ=255)
             ggTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAATTTAAAAAAAAt                                                                               <  2:1118463/91‑1 (MQ=255)
              gTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATa                                                                              <  1:2961010/91‑1 (MQ=255)
               tATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATaa                                                                             >  1:1765100/1‑91 (MQ=255)
                 ttttATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATaaaa                                                                           >  1:1703459/1‑91 (MQ=255)
                  tttATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATaaaaa                                                                          <  2:3798566/91‑1 (MQ=255)
                  tttATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTGTTAAGTGATAACTTAAAAACAATaaaaa                                                                          >  2:1177792/1‑91 (MQ=255)
                    tATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAAt                                                                        >  1:2802658/1‑91 (MQ=255)
                      tCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATta                                                                      >  1:437114/1‑91 (MQ=255)
                       cAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATtat                                                                     <  1:628119/91‑1 (MQ=255)
                        aaaaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATtata                                                                    <  2:1839735/91‑1 (MQ=255)
                        aaaaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATtata                                                                    >  1:3575935/1‑91 (MQ=255)
                          aaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATaaa                                                                  <  2:2968740/91‑1 (MQ=255)
                          aaaTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATaaa                                                                  <  1:2111143/91‑1 (MQ=255)
                            aTTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatt                                                                >  1:2385625/1‑91 (MQ=255)
                             ttCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatta                                                               <  2:2706123/91‑1 (MQ=255)
                             ttCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatta                                                               >  2:3648908/1‑91 (MQ=255)
                             ttCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaatta                                                               <  1:2553608/91‑1 (MQ=255)
                              tCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaa                                                              <  2:1661996/91‑1 (MQ=255)
                               cTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaat                                                             <  2:3329631/91‑1 (MQ=255)
                               cTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaat                                                             <  1:1275255/91‑1 (MQ=255)
                               cTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAaattaat                                                             <  1:1277840/91‑1 (MQ=255)
                                 gTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAattaatta                                                           >  2:2443118/1‑91 (MQ=255)
                                 gTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAattaatta                                                           <  1:699415/91‑1 (MQ=255)
                                  ttCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAg                                                          >  1:3501789/1‑91 (MQ=255)
                                  ttCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAg                                                          >  2:1168579/1‑91 (MQ=255)
                                    cATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGtt                                                        <  1:1128177/91‑1 (MQ=255)
                                     aTTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTa                                                       <  1:3246685/91‑1 (MQ=255)
                                        agaagaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGt                                                    <  1:1918912/91‑1 (MQ=255)
                                        agaagaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGt                                                    >  1:1678640/1‑91 (MQ=255)
                                          aagaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGa                                                  <  1:1003987/91‑1 (MQ=255)
                                          aagaTGGACAAGCGTCACAACAAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAAAAAAAAAAAATTATAAATTAATTAGTTAAGTGa                                                  <  2:2652216/91‑1 (MQ=255)
                                            gaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGata                                                >  1:3734337/1‑91 (MQ=255)
                                            gaTGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGata                                                <  1:3348943/91‑1 (MQ=255)
                                                gACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataaa                                            <  2:171912/91‑1 (MQ=255)
                                                gACAAGCGTCACAACGAAATTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataaa                                            <  2:1779368/91‑1 (MQ=255)
                                                 aCAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAACATTATAAATTAATTAGTTAAGTGATataaat                                           >  2:2425521/1‑91 (MQ=255)
                                                  cAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataaata                                          >  1:203541/1‑91 (MQ=255)
                                                   aaGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataaataa                                         >  1:873222/1‑91 (MQ=255)
                                                   aaGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATataaataa                                         <  1:201878/91‑1 (MQ=255)
                                                       gTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATaaa                                     <  2:2048610/91‑1 (MQ=255)
                                                        tCACAACGAATTTGCCAAACGATTTTTAAGTGCTAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAAc                                    >  2:3813481/1‑91 (MQ=255)
                                                        acacaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAAc                                    <  2:2171063/90‑1 (MQ=255)
                                                          acaACGAATTTGCCAAACGATTTTTAAGTTATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                                  <  2:321296/91‑1 (MQ=255)
                                                          acaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                                  <  2:1457013/91‑1 (MQ=255)
                                                          acaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                                  >  1:895021/1‑91 (MQ=255)
                                                          acaACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGa                                  <  1:271405/91‑1 (MQ=255)
                                                           caACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaa                                 <  2:44978/91‑1 (MQ=255)
                                                            aaCGAATTTGCCAAACGATTTTTAAGTGATAACTTAGAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaaa                                <  2:3062344/91‑1 (MQ=255)
                                                            aaCGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATCAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaaa                                <  2:3475861/91‑1 (MQ=255)
                                                            aaCGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGaaa                                >  2:2838362/1‑91 (MQ=255)
                                                             aCGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAAt                               >  2:596222/1‑91 (MQ=255)
                                                             aCGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAAt                               >  1:3178328/1‑91 (MQ=255)
                                                             aCGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAAt                               <  1:1038734/91‑1 (MQ=255)
                                                              cGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATg                              >  2:326820/1‑91 (MQ=255)
                                                              cGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATg                              <  2:1326535/91‑1 (MQ=255)
                                                              cGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATg                              <  2:2906350/91‑1 (MQ=255)
                                                                 aTTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGttt                           >  2:900758/1‑91 (MQ=255)
                                                                   ttGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAACCGAAATGTTTGc                         >  2:3227053/1‑91 (MQ=255)
                                                                   ttGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGc                         <  1:1513652/91‑1 (MQ=255)
                                                                    tGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCt                        <  1:747428/91‑1 (MQ=255)
                                                                     gCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCtt                       >  1:53103/1‑91 (MQ=255)
                                                                      ccAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTg                      >  2:2481525/1‑91 (MQ=255)
                                                                      ccAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTCTGCTTg                      <  2:3496688/91‑1 (MQ=255)
                                                                        aaaCGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTa                    <  2:182255/91‑1 (MQ=255)
                                                                        aaaCGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTa                    <  2:2333740/91‑1 (MQ=255)
                                                                          aCGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATg                  <  2:2683343/91‑1 (MQ=255)
                                                                          aCGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATg                  <  2:736647/91‑1 (MQ=255)
                                                                          aCGATTTTTAAGTGATAACTTAAAAAAAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATg                  <  2:3741090/91‑1 (MQ=255)
                                                                             aTTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGtta               <  2:2385625/91‑1 (MQ=255)
                                                                              tttttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttat              >  2:2079367/1‑91 (MQ=255)
                                                                              tttttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttat              >  1:3425147/1‑91 (MQ=255)
                                                                              tttttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttat              >  2:3501805/1‑91 (MQ=255)
                                                                                tttAAGTTATCACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta            >  2:1933018/1‑91 (MQ=255)
                                                                                tttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta            <  1:371854/91‑1 (MQ=255)
                                                                                tttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGttatta            >  1:3548499/1‑91 (MQ=255)
                                                                                 ttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTtattat           <  1:827958/91‑1 (MQ=255)
                                                                                 ttAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTtattat           >  2:1237748/1‑91 (MQ=255)
                                                                                   aaGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTAttt         <  2:1119338/91‑1 (MQ=255)
                                                                                      tGATAACATAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTtgtg      <  2:3736004/91‑1 (MQ=255)
                                                                                        aTAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTa    <  1:841066/91‑1 (MQ=255)
                                                                                        aTAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTa    >  1:381494/1‑91 (MQ=255)
                                                                                        aTAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTa    <  2:1403586/91‑1 (MQ=255)
                                                                                          aaCTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTATg  <  2:2316715/91‑1 (MQ=255)
                                                                                          aaCTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTATg  >  1:3050888/1‑91 (MQ=255)
                                                                                          aaCTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTATg  <  2:314543/91‑1 (MQ=255)
                                                                                          |                                                                                          
TTGAAGCACGAACGGTATTTTATCAAAAATTCTGTTCATTAGAAGATGGACAAGCGTCACAACGAATTTGCCAAACGATTTTTAAGTGATAACTTAAAAACAATAAAAAATTATAAATTAATTAGTTAAGTGATATAAATAATAAACGAAATGTTTGCTTGTATGTTATTATTTGTGTATG  >  USA300TCH1516_ALE/299370‑299550

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: