Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,619 Δ1 bp 7.2% intergenic (+55/+85) pepF1_2 → / ← phoU Oligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,421,6190T. 7.2% 855.4 / 54.3 262intergenic (+55/+85)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base T (124/116);  new base . (15/4);  total (141/121)
Fisher's exact test for biased strand distribution p-value = 2.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

AAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTG  >  USA300TCH1516_ALE/1421529‑1421709
                                                                                           |                                                                                          
 aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                                           <  1:459555/91‑1 (MQ=255)
 aaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAtt                                                                                           >  2:2469107/1‑91 (MQ=255)
  aTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtg                                                                                          >  1:537807/1‑91 (MQ=255)
   ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                         <  1:1505755/91‑1 (MQ=255)
   ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                         <  2:2874761/91‑1 (MQ=255)
     aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                       >  1:1740393/1‑91 (MQ=255)
      tcAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      <  2:3927884/90‑1 (MQ=255)
      aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      >  2:1339813/1‑91 (MQ=255)
      aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      >  2:2016187/1‑91 (MQ=255)
      aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      >  1:3640930/1‑91 (MQ=255)
      aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      >  1:3499566/1‑91 (MQ=255)
        aGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCCCATTACTAATTGTGAGaa                                                                                    >  2:2783442/1‑91 (MQ=255)
         gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                                   <  2:657978/91‑1 (MQ=255)
         gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                                   >  1:1704375/1‑91 (MQ=255)
          aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                                  <  2:83355/91‑1 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              <  2:1757071/91‑1 (MQ=255)
               aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                             <  2:1401605/91‑1 (MQ=255)
                 tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                           >  2:3463679/1‑91 (MQ=255)
                    gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATACACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                        <  1:2131477/91‑1 (MQ=255)
                    gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                        <  1:3432738/91‑1 (MQ=255)
                     caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                                       >  2:3564908/1‑91 (MQ=255)
                     caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                                       >  2:1532282/1‑91 (MQ=255)
                     caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                                       <  2:488310/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  2:1857342/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  1:1098936/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  2:2146947/91‑1 (MQ=255)
                         caaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAgg                                                                   <  1:2281496/91‑1 (MQ=255)
                            aTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGtt                                                                <  2:2086903/91‑1 (MQ=255)
                            aTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGtt                                                                >  1:675605/1‑91 (MQ=255)
                             tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                               <  1:3926344/91‑1 (MQ=255)
                               aCTAAGATGCGAGACAATTTATAGTAGCACAATATAAACCAAGATTTTCACATTACTAAGTGTGAGAATCTTGGGTTTGATGAGGGGTTAt                                                             >  1:3585246/1‑91 (MQ=255)
                                cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                            >  1:1481480/1‑91 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                         >  2:2492960/1‑91 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                         <  1:1580552/91‑1 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                         <  2:1276230/91‑1 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                         <  1:2882078/91‑1 (MQ=255)
                                     aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                                       <  1:673818/91‑1 (MQ=255)
                                     aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                                       <  1:3774040/91‑1 (MQ=255)
                                      tGCGAGACAATTTATAGTAGAACAATATAAACTAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataat                                                      <  1:3370748/91‑1 (MQ=255)
                                          agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                                  >  1:2154306/1‑91 (MQ=255)
                                          agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                                  >  1:2576916/1‑91 (MQ=255)
                                           gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                                 <  1:3204035/91‑1 (MQ=255)
                                           gaCAATGTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                                 <  2:1076708/91‑1 (MQ=255)
                                                  tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                                          >  1:1201671/1‑91 (MQ=255)
                                                   ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                                         >  2:589462/1‑91 (MQ=255)
                                                   ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                                         >  2:2673720/1‑91 (MQ=255)
                                                       tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCCTGGTTTTTTTGTGGGTTTTTGAATAATAGAAATAGATATTTATc                                     >  2:2543210/1‑91 (MQ=255)
                                                        agAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGCTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCg                                    >  2:2327535/1‑91 (MQ=255)
                                                           aCAATATAAACCAATATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                                 <  2:2128450/91‑1 (MQ=255)
                                                           aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                                 <  1:2695533/91‑1 (MQ=255)
                                                            cAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAg                                <  1:1353270/91‑1 (MQ=255)
                                                             aaTATAAACCAAGGTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                               <  2:2381812/91‑1 (MQ=255)
                                                             aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                               <  2:1907179/91‑1 (MQ=255)
                                                             aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                               <  2:310216/91‑1 (MQ=255)
                                                             aaTATAAACCAAGATGTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                               <  1:1283752/91‑1 (MQ=255)
                                                               tataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGa                             >  1:756261/1‑91 (MQ=255)
                                                                ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATCCTTATCGTAAGGGAt                            >  2:646929/1‑91 (MQ=255)
                                                                ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                            >  2:901977/1‑91 (MQ=255)
                                                                ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                            <  1:3729485/91‑1 (MQ=255)
                                                                  aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                          >  1:1892842/1‑91 (MQ=255)
                                                                  aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                          <  2:3738257/91‑1 (MQ=255)
                                                                  aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                          >  1:653130/1‑91 (MQ=255)
                                                                   aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                         >  2:1781863/1‑91 (MQ=255)
                                                                   aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGACATATATATTTATCGTAAGGGATTTta                        >  1:2871837/1‑91 (MQ=18)
                                                                   aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                        >  1:1002238/1‑91 (MQ=18)
                                                                    aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                        >  2:3136506/1‑91 (MQ=255)
                                                                    aCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttga                        >  1:819584/1‑89 (MQ=21)
                                                                    aCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                        >  1:3258251/1‑91 (MQ=21)
                                                                        aGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                    >  2:1363444/1‑91 (MQ=255)
                                                                        aGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                    <  2:3300074/91‑1 (MQ=255)
                                                                        aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGCGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                    >  2:3003447/1‑91 (MQ=21)
                                                                        aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTCAGGGATTTTATAGt                    >  2:2627261/1‑91 (MQ=21)
                                                                          aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa                  <  1:2390302/91‑1 (MQ=255)
                                                                           ttttCACGTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat                 <  2:3510106/91‑1 (MQ=255)
                                                                           ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGACTAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat                 <  2:3001485/91‑1 (MQ=255)
                                                                           ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat                 >  1:3400636/1‑91 (MQ=255)
                                                                           ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat                 >  1:3000722/1‑91 (MQ=255)
                                                                            tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata                >  2:375319/1‑91 (MQ=255)
                                                                            tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata                >  2:3402215/1‑91 (MQ=255)
                                                                             ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa               >  2:3395743/1‑91 (MQ=255)
                                                                             ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa               >  1:3165704/1‑91 (MQ=255)
                                                                               cacaTTACTAATTGTGAGCATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt             >  1:2171763/1‑91 (MQ=255)
                                                                               cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt             >  2:2822258/1‑91 (MQ=255)
                                                                               cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt             >  2:727675/1‑91 (MQ=255)
                                                                               cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt             >  1:1051044/1‑91 (MQ=255)
                                                                                acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            >  1:1251581/1‑91 (MQ=255)
                                                                                acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            >  1:1530486/1‑91 (MQ=255)
                                                                                acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            >  1:2344572/1‑91 (MQ=255)
                                                                                 caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACttt           <  1:2928046/91‑1 (MQ=255)
                                                                                     tACTAATTGTGAGAAGCTTGGTTTTGGTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAg        <  2:3422998/91‑1 (MQ=255)
                                                                                      aCTAATTGTGAGAATCTTGGTTTTTTTTTTTTTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt       >  2:2207721/1‑91 (MQ=255)
                                                                                      aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt       >  1:632754/1‑91 (MQ=255)
                                                                                       cgaaTTGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGtt      <  2:1170823/89‑1 (MQ=255)
                                                                                         aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAt    <  1:1944738/91‑1 (MQ=255)
                                                                                         aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAt    <  1:885341/91‑1 (MQ=255)
                                                                                         aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAt    <  2:1739184/91‑1 (MQ=255)
                                                                                          aTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAtt   >  1:1834181/1‑91 (MQ=255)
                                                                                          aTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAtt   <  2:3173499/91‑1 (MQ=255)
                                                                                          aTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAtt   >  1:2515089/1‑91 (MQ=255)
                                                                                           ttGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTg  >  2:3600620/1‑91 (MQ=255)
                                                                                           |                                                                                          
AAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTG  >  USA300TCH1516_ALE/1421529‑1421709

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