Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,808,256 T→G 5.3% intergenic (‑156/+113) 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,2560TG 5.3% 704.4 / 20.9 208intergenic (‑156/+113)pfkA/accAATP‑dependent 6‑phosphofructokinase/Acetyl‑coenzyme A carboxylase carboxyl transferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (91/106);  new base G (4/7);  total (95/113)
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.98e-01

Only 100 of 218 total aligned reads displayed.

CCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTATCGATT  >  USA300TCH1516_ALE/1808166‑1808345
                                                                                          |                                                                                         
ccATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAAttt                                                                                           <  2:905088/91‑1 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATTTTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  2:558302/1‑91 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  2:1209748/1‑91 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  1:1454384/1‑91 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  1:2265967/1‑91 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  1:212097/1‑91 (MQ=255)
   tCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGa                                                                                        >  1:3277109/1‑91 (MQ=255)
    cAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGac                                                                                       >  2:3589684/1‑91 (MQ=255)
    cAAAACTAATGTATTGAATTGATGTAAACATTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACTATTTTGac                                                                                       >  2:3380004/1‑91 (MQ=255)
      aaaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGacta                                                                                     <  2:1658757/91‑1 (MQ=255)
      aaaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGacta                                                                                     <  2:1500873/91‑1 (MQ=255)
      aaaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGacta                                                                                     <  1:1929463/91‑1 (MQ=255)
       aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTtactac                                                                                    >  2:3622863/1‑91 (MQ=255)
       aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                                    <  1:338229/91‑1 (MQ=255)
       aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                                    <  1:1756152/91‑1 (MQ=255)
       aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                                    >  1:880566/1‑91 (MQ=255)
       aaCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactac                                                                                    <  2:3775999/91‑1 (MQ=255)
        aCTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGactact                                                                                   <  2:2063397/91‑1 (MQ=255)
          tAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTACCTTAATTAAAAAAAGATTAGTAGTAAAAATTGTGGCTACTaa                                                                                 <  1:593690/91‑1 (MQ=37)
           aaTGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAAt                                                                                >  2:1073634/1‑91 (MQ=255)
           aaTGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAAt                                                                                >  2:3105354/1‑91 (MQ=255)
               tATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGTCAAAATTGTGGCTACTAATCttt                                                                            <  1:1540881/91‑2 (MQ=25)
               tATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTa                                                                            <  1:77997/91‑1 (MQ=255)
                attgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGTCAAAATTGTGGCTACTAATCtttt                                                                           <  1:180964/91‑3 (MQ=25)
                attgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAg                                                                           <  1:1834255/91‑1 (MQ=255)
                attgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAg                                                                           >  2:3330985/1‑91 (MQ=255)
                 ttgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGt                                                                          <  2:3875898/91‑1 (MQ=255)
                 ttgaattgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGt                                                                          <  2:251930/91‑1 (MQ=255)
                      ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                     <  2:1536655/91‑1 (MQ=255)
                      ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                     <  2:3550837/91‑1 (MQ=255)
                      ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                     <  2:3550651/91‑1 (MQ=255)
                      ttgaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTaaaa                                                                     <  1:3802220/91‑1 (MQ=255)
                        gaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTCGTTAAAAtt                                                                   >  1:3717310/1‑91 (MQ=255)
                        gaTGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAAtt                                                                   >  2:474974/1‑91 (MQ=255)
                          tGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCa                                                                 >  2:3195258/1‑91 (MQ=255)
                          tGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCa                                                                 >  1:3462794/1‑91 (MQ=255)
                           gTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAg                                                                <  2:1584992/91‑1 (MQ=255)
                              aaCGTTTTAAATATTGAATTTTATTTTATGTTAATTAAAAAAAGATTAGTAGCCATATTTTTGACTATTAATGTTAGTTAAAATTCAGTAt                                                             >  2:3500573/1‑91 (MQ=255)
                              aaCGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTAt                                                             <  1:844727/91‑1 (MQ=255)
                              aaCGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTAt                                                             >  2:1232322/1‑91 (MQ=255)
                                 tttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTa                                                          <  2:2274591/90‑1 (MQ=255)
                                  ttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTaa                                                         >  1:563808/1‑91 (MQ=255)
                                  ttttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTaa                                                         >  2:3828409/1‑91 (MQ=255)
                                    ttAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAat                                                       <  2:600477/91‑1 (MQ=255)
                                     tAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAata                                                      <  1:1329795/91‑1 (MQ=255)
                                      aaaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAatat                                                     <  2:2554987/91‑1 (MQ=255)
                                      aaaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAatat                                                     <  1:586228/91‑1 (MQ=255)
                                      aaaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAatat                                                     <  2:813466/91‑1 (MQ=255)
                                       aaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATAtt                                                    >  2:1736105/1‑91 (MQ=255)
                                       aaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATAtt                                                    >  1:2175838/1‑91 (MQ=255)
                                       aaTATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATAtt                                                    >  2:1016294/1‑91 (MQ=255)
                                        atatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTa                                                   <  1:1017106/91‑1 (MQ=255)
                                        atatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTa                                                   <  1:1163770/91‑1 (MQ=255)
                                        atatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATATTAGTTAAAATTCAGTATCTAAATATTa                                                   >  2:1667883/1‑91 (MQ=255)
                                         tatTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTaa                                                  >  2:3702100/1‑91 (MQ=255)
                                           ttGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAat                                                <  1:807398/91‑1 (MQ=255)
                                           ttGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAat                                                >  1:335690/1‑91 (MQ=255)
                                            tGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAatt                                               <  2:920813/91‑1 (MQ=255)
                                             gAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAatta                                              <  1:2599547/91‑1 (MQ=255)
                                               attttattttAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattatt                                            <  2:129172/91‑1 (MQ=255)
                                               attttattttAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattatt                                            >  1:3426570/1‑91 (MQ=255)
                                               attttattttAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTAAAAATTCAGTATCTAAATATTAAattatt                                            <  2:275176/91‑1 (MQ=255)
                                               attttattttAAAAAAAATAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAattatt                                            <  2:3768071/91‑1 (MQ=255)
                                                ttttattttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTg                                           >  2:504045/1‑91 (MQ=255)
                                                 tttattttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGt                                          <  1:3169404/91‑1 (MQ=255)
                                                     ttttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTa                                      >  2:1171388/1‑91 (MQ=255)
                                                      tttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTaa                                     >  2:3565181/1‑91 (MQ=255)
                                                      tttaACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTaa                                     >  2:1863403/1‑91 (MQ=255)
                                                         aaCTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGt                                  <  2:880566/91‑1 (MQ=255)
                                                           cTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTAc                                <  1:3220753/91‑1 (MQ=255)
                                                           cTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTAc                                >  1:2095965/1‑91 (MQ=255)
                                                              aattaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATc                                                            <  1:1604781/60‑1 (MQ=255)
                                                               attaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCa                            <  2:1714429/91‑1 (MQ=255)
                                                               attaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCa                            >  2:708476/1‑91 (MQ=255)
                                                                ttaaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCaa                           >  1:1955534/1‑91 (MQ=255)
                                                                 taaAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAAt                          >  2:301588/1‑91 (MQ=255)
                                                                  aaaaaaaGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAAtt                         <  1:2215072/91‑1 (MQ=255)
                                                                        aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt                   <  1:2542918/91‑1 (MQ=255)
                                                                        aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt                   >  1:1917992/1‑91 (MQ=255)
                                                                        aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt                   <  1:2616823/91‑1 (MQ=255)
                                                                        aGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGt                   <  2:3594981/91‑1 (MQ=255)
                                                                         gATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTa                  >  1:2961456/1‑91 (MQ=255)
                                                                          aTTAGTAGTCAAAATTGTGGCTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTCCTTCAATTTATGGTaa                 >  1:306818/1‑91 (MQ=37)
                                                                           ttAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaa                >  1:3766882/1‑91 (MQ=255)
                                                                           ttAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaa                >  2:596144/1‑91 (MQ=255)
                                                                            tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa               <  2:27841/91‑1 (MQ=255)
                                                                            tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa               <  2:2628690/91‑1 (MQ=255)
                                                                            tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa               >  2:2283601/1‑91 (MQ=255)
                                                                            tagtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTaaaa               >  2:776696/1‑91 (MQ=255)
                                                                             agtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAc              >  1:2659261/1‑91 (MQ=255)
                                                                             agtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAc              >  1:493964/1‑91 (MQ=255)
                                                                              gtagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAcc             >  2:2152845/1‑91 (MQ=255)
                                                                               tagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAccc            <  1:3436337/91‑1 (MQ=255)
                                                                               tagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAccc            <  1:3589684/91‑1 (MQ=255)
                                                                               tagCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAAccc            <  1:762371/91‑1 (MQ=255)
                                                                                  ccACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGt         <  2:2265967/91‑1 (MQ=255)
                                                                                   cacaATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTa        <  2:1417463/91‑1 (MQ=255)
                                                                                       aTTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTATCGa    <  1:108967/91‑1 (MQ=255)
                                                                                        ttttGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTATCGAt   <  1:3105354/91‑1 (MQ=255)
                                                                                        ttttGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTATCGAt   >  1:1400468/1‑91 (MQ=255)
                                                                                          |                                                                                         
CCATCAAAACTAATGTATTGAATTGATGTAAACGTTTTAAATATTGAATTTTATTTTAACTTAATTAAAAAAAGATTAGTAGCCACAATTTTGACTACTAATCTTAGTTAAAATTCAGTATCTAAATATTAAATTATTGTTTTAATGTACTTCAATTTATGGTAAAACCCAGTATCGATT  >  USA300TCH1516_ALE/1808166‑1808345

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