Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,808,251 C→A 6.1% intergenic (‑151/+118) pfkA ← / ← accA ATP‑dependent 6‑phosphofructokinase/Acetyl‑coenzyme A carboxylase carboxyl transferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,808,2510CA 6.1% 741.4 / 27.6 215intergenic (‑151/+118)pfkA/accAATP‑dependent 6‑phosphofructokinase/Acetyl‑coenzyme A carboxylase carboxyl transferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (94/107);  new base A (6/7);  total (101/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 227 total aligned reads displayed.

TCATTCCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTAT  >  USA300TCH1516_ALE/1808161‑1808340
                                                                                          |                                                                                         
tCATTCCATCAAAACTAATGTATTGTATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCcac                                                                                           <  1:3611088/91‑1 (MQ=255)
tCATTCCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCcac                                                                                           >  2:1690552/1‑91 (MQ=255)
    tCCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAAtt                                                                                       >  1:3324422/1‑91 (MQ=255)
    tCCATCAAAACTAATGTATTGAATTGATGTAAACGTTATAAATATTGAATTTTATTTTTCCTTGTTTAAAAAAAGCTTAGTAGCTACAAtt                                                                                       >  1:980265/1‑91 (MQ=255)
     ccATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAAttt                                                                                      <  1:1557139/91‑1 (MQ=255)
     ccATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAAttt                                                                                      <  2:905088/91‑1 (MQ=255)
      cATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAAtttt                                                                                     >  1:2652510/1‑91 (MQ=255)
       aTCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTg                                                                                    <  2:1947864/91‑1 (MQ=255)
        tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATTTTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                   >  2:558302/1‑91 (MQ=255)
        tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGTCAAAAttgtgg                                                                                   >  1:390931/1‑87 (MQ=255)
        tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCcac                                                                                           <  2:3089049/83‑1 (MQ=255)
        tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                   >  1:212097/1‑91 (MQ=255)
        tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                   >  1:2265967/1‑91 (MQ=255)
         cAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGac                                                                                  <  2:1536087/91‑1 (MQ=255)
         cAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGac                                                                                  >  2:3589684/1‑91 (MQ=255)
         cAAAACTAATGTATTGAATTGATGTAAACATTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACTATTTTGac                                                                                  >  2:3380004/1‑91 (MQ=255)
           aaaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGacta                                                                                >  2:1341364/1‑91 (MQ=255)
           aaaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGacta                                                                                <  2:1659874/91‑1 (MQ=255)
            aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                               <  1:1756168/91‑1 (MQ=255)
            aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                               <  1:338229/91‑1 (MQ=255)
             aCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGTCAAAATTGTGgctact                                                                              <  2:3107459/91‑1 (MQ=38)
               tAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTACCTTAATTAAAAAAAGATTAGTAGTAAAAATTGTGGCTACTaa                                                                            <  1:593690/91‑1 (MQ=37)
                aaTGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAAt                                                                           >  2:1073634/1‑91 (MQ=255)
                aaTGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAAt                                                                           >  1:1176359/1‑91 (MQ=255)
                   gTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGTCAAAATTGTGGCTACTAATCtt                                                                        <  2:2035531/91‑1 (MQ=37)
                    tATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTa                                                                       <  1:77997/91‑1 (MQ=255)
                    tATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTa                                                                       >  1:3515467/1‑91 (MQ=255)
                     attgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAg                                                                      >  2:3330985/1‑91 (MQ=255)
                     attgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAg                                                                      <  1:1834255/91‑1 (MQ=255)
                      ttgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGt                                                                     <  2:3875898/91‑1 (MQ=255)
                      ttgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGt                                                                     <  2:2554083/91‑1 (MQ=255)
                       tgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGtt                                                                    <  1:1533281/91‑1 (MQ=255)
                           ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                <  1:3802220/91‑1 (MQ=255)
                           ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                <  2:3550837/91‑1 (MQ=255)
                           ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                <  2:1536742/91‑1 (MQ=255)
                           ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                <  2:3550651/91‑1 (MQ=255)
                            tgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAAt                                                               <  2:2589280/91‑1 (MQ=255)
                               tGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCa                                                            >  1:3462380/1‑91 (MQ=255)
                                gTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAg                                                           <  2:1584992/91‑1 (MQ=255)
                                   aaCGTTTTAAATATTGAATTTTATTTTATGTTAATTAAAAAAAGATTAGTAGCCATATTTTTGACTATTAATGTTAGTTAAAATTCAGTAt                                                        >  2:3500573/1‑91 (MQ=255)
                                      tttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTa                                                     <  2:2274591/90‑1 (MQ=255)
                                      tttatAAATATTGAATTTTATTTTAACTAAATTAAAAAAAGATTAGTAGCCACAAATTTTACTACAAATCTTAGTAAAAATTCAGTACCTa                                                     <  2:1837089/87‑1 (MQ=255)
                                       ttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTaa                                                    >  1:563808/1‑91 (MQ=255)
                                       ttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTaa                                                    >  2:3828409/1‑91 (MQ=255)
                                         ttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAat                                                  <  2:600477/91‑1 (MQ=255)
                                           aaaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAatat                                                <  1:586228/91‑1 (MQ=255)
                                           aaaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAatat                                                <  1:3195258/91‑1 (MQ=255)
                                            aaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATAtt                                               >  1:2175838/1‑91 (MQ=255)
                                            aaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATAtt                                               >  2:1016294/1‑91 (MQ=255)
                                             atatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTa                                              <  1:1017106/91‑1 (MQ=255)
                                              tatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTaa                                             >  2:3702100/1‑91 (MQ=255)
                                              tatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTaa                                             <  1:876683/91‑1 (MQ=255)
                                              tatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTaa                                             <  2:2495180/91‑1 (MQ=255)
                                               atTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTCATCTTAGTTAAAATTCAGTATCTAAATATTaca                                            >  2:2744364/1‑89 (MQ=255)
                                                ttGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAat                                           <  1:807398/91‑1 (MQ=255)
                                                ttGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAat                                           >  1:335690/1‑91 (MQ=255)
                                                 tGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAatt                                          <  2:920813/91‑1 (MQ=255)
                                                 tGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAatt                                          <  1:1108209/91‑1 (MQ=255)
                                                  gAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAatta                                         <  1:2599547/91‑1 (MQ=255)
                                                   aaTTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattat                                        <  2:683869/91‑1 (MQ=255)
                                                    attttattttAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattatt                                       >  1:3426570/1‑91 (MQ=255)
                                                    attttattttAAAAAAAATAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattatt                                       <  2:3768071/91‑1 (MQ=255)
                                                     ttttattttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTg                                      >  2:504045/1‑91 (MQ=255)
                                                      tttattttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGt                                     <  1:3169278/91‑1 (MQ=255)
                                                           tttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTaa                                >  2:3565181/1‑91 (MQ=255)
                                                             taACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATg                              <  1:1919948/91‑1 (MQ=255)
                                                              aaCTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGt                             <  2:880566/91‑1 (MQ=255)
                                                              aaCTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTAAGTTAAAATTCAGTATCTATATATTAAATTATTGTTTTAATGt                             >  2:3613011/1‑91 (MQ=255)
                                                                cTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTAc                           >  1:2095965/1‑91 (MQ=255)
                                                                   aattaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATc                                                       >  2:1604781/1‑60 (MQ=255)
                                                                    attaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCa                       <  2:1714429/91‑1 (MQ=255)
                                                                    attaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCa                       >  2:709820/1‑91 (MQ=255)
                                                                     ttaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCaa                      >  1:1955533/1‑91 (MQ=255)
                                                                      taaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAAt                     >  2:301588/1‑91 (MQ=255)
                                                                       aaaaaaaGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAAtt                    <  1:2215072/91‑1 (MQ=255)
                                                                          aaaaGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTAt                 <  1:169234/91‑1 (MQ=255)
                                                                             aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt              <  1:2542918/91‑1 (MQ=255)
                                                                             aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt              >  1:1917992/1‑91 (MQ=255)
                                                                              gATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTa             >  1:2961456/1‑91 (MQ=255)
                                                                              gATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTa             >  2:1425146/1‑91 (MQ=255)
                                                                              gATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTa             >  2:435439/1‑91 (MQ=255)
                                                                               aTTAGTAGTCAAAATTGTGGCTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTCCTTCAATTTATGGTaa            >  1:306818/1‑91 (MQ=37)
                                                                               aTTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaa            <  1:1352606/91‑1 (MQ=255)
                                                                                ttAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaa           >  1:2933083/1‑91 (MQ=255)
                                                                                ttAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaa           >  1:3766882/1‑91 (MQ=255)
                                                                                ttAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaa           >  2:596144/1‑91 (MQ=255)
                                                                                 tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa          <  2:27841/91‑1 (MQ=255)
                                                                                 tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa          >  2:776696/1‑91 (MQ=255)
                                                                                 tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa          >  2:2283601/1‑91 (MQ=255)
                                                                                  agtagtCAAAATTGTGGCTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAc         >  2:1264353/1‑91 (MQ=38)
                                                                                  agtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAc         >  1:493964/1‑91 (MQ=255)
                                                                                   gtagTCAAAATTGTGGCTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAcc        >  2:1217799/1‑91 (MQ=38)
                                                                                   gtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAcc        >  2:2152845/1‑91 (MQ=255)
                                                                                    tagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAccc       >  2:364911/1‑91 (MQ=255)
                                                                                    tagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAccc       <  1:3589684/91‑1 (MQ=255)
                                                                                       ccACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGt    <  2:2265967/91‑1 (MQ=255)
                                                                                       ccACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGt    <  1:837380/91‑1 (MQ=255)
                                                                                       ccACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGt    <  1:3240498/91‑1 (MQ=255)
                                                                                        cacaATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTa   >  2:810216/1‑91 (MQ=255)
                                                                                         acaATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTAt  <  2:1176359/91‑1 (MQ=255)
                                                                                          |                                                                                         
TCATTCCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTAT  >  USA300TCH1516_ALE/1808161‑1808340

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: