Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,973 C→T 5.8% intergenic (+186/+176) asbF → / ← USA300HOU_RS01125 3‑dehydroshikimate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE253,9730CT 5.8% 802.5 / 23.2 241intergenic (+186/+176)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (120/107);  new base T (4/10);  total (124/117)
Fisher's exact test for biased strand distribution p-value = 9.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 248 total aligned reads displayed.

CACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTT  >  USA300TCH1516_ALE/253883‑254063
                                                                                           |                                                                                          
cacCTTAAATCGAGTCATTGTTTATAAAGAAGGTGTACTTTATATAAGTGTATAGCGATTGTCATACCCATTCACAGTAACAATCCTCAcc                                                                                            <  2:3012550/91‑1 (MQ=255)
 acCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCa                                                                                           >  2:3939881/1‑91 (MQ=255)
  ccTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCAt                                                                                          <  2:624935/91‑1 (MQ=255)
  ccTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCAt                                                                                          <  2:2806199/91‑1 (MQ=255)
   cTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCCCCAtt                                                                                         <  2:2715459/91‑1 (MQ=255)
   cTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCAtt                                                                                         <  1:1094872/91‑1 (MQ=255)
   cTTAAAGCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCAtt                                                                                         <  2:110394/91‑1 (MQ=255)
         tCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAag                                                                                   >  2:2491056/1‑91 (MQ=255)
          cAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAaga                                                                                  >  2:2741016/1‑91 (MQ=255)
            aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                >  2:903740/1‑91 (MQ=255)
            aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                >  1:469445/1‑91 (MQ=255)
            aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                >  1:1241316/1‑91 (MQ=255)
             gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                               <  1:2795328/91‑1 (MQ=255)
             gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                               >  1:2220240/1‑91 (MQ=255)
                aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtttat                                                                            >  2:1079139/1‑91 (MQ=255)
                aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                            >  1:968139/1‑91 (MQ=255)
                aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                            >  1:802175/1‑91 (MQ=255)
                aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                            <  1:2736964/91‑1 (MQ=255)
                 ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatata                                                                           <  1:236489/91‑1 (MQ=255)
                 ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                           >  1:2561127/1‑91 (MQ=255)
                  tgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATaa                                                                          >  2:1678741/1‑91 (MQ=255)
                   gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                         >  1:560900/1‑91 (MQ=255)
                    tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                        <  1:1289295/91‑4 (MQ=255)
                    tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                        <  2:2197397/91‑1 (MQ=255)
                    tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                        <  2:3836560/91‑1 (MQ=255)
                      tgtAAGGAAGGTGTACTTTATATAAGTAGATAGCGACGGTCATACCCATTCACAGTATCGATCCTCACCATTGAAAGGAGTATATAACCtt                                                                      <  2:536913/91‑1 (MQ=255)
                       gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                     <  1:3835372/91‑1 (MQ=255)
                        tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                    <  1:3495494/91‑1 (MQ=255)
                        tAAAGAAGGTGTACTTTATAGAAGTATATAGCGATGGTCATACCCATTCACAGTAGCAATCCTCACTATTGAAAAGGTTATATACTCtttt                                                                    <  2:1547909/91‑1 (MQ=25)
                         aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                   >  1:407/1‑91 (MQ=255)
                         aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                   >  2:1783545/1‑91 (MQ=255)
                          aagaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                  <  2:551134/91‑1 (MQ=255)
                           agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                 <  2:3903558/91‑1 (MQ=25)
                            gacggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAGTCCTCACTATTGAAAAGGTTATATACTCTTTTCAAt                                                                <  2:3145141/88‑1 (MQ=21)
                             aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtg                                                               <  1:1464236/91‑2 (MQ=25)
                             aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                               >  2:3132607/1‑91 (MQ=255)
                              agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                              <  1:3815048/91‑1 (MQ=255)
                                gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                            >  1:3938669/1‑91 (MQ=255)
                                   tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                         >  1:753780/1‑91 (MQ=255)
                                     acttATATAAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                        >  1:1247476/3‑91 (MQ=255)
                                        tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                    <  2:3309184/91‑1 (MQ=255)
                                        tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                    <  2:1328875/91‑1 (MQ=255)
                                         atctaAGTATATAGCGATGGTCATACCCATCCACAGTACCAATCCTCACCATTGCAAACCGCCTATAACCTTTTCAACCGTGCGGTATATg                                                   >  2:442158/4‑91 (MQ=255)
                                           ataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                 <  1:2226045/91‑1 (MQ=255)
                                            taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                <  1:1340952/91‑1 (MQ=255)
                                               aGTCTATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGgtaat                                              <  1:3122341/91‑1 (MQ=255)
                                               aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                              >  1:3845178/1‑91 (MQ=255)
                                               aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                              >  1:1307334/1‑91 (MQ=255)
                                               aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                              >  2:3067573/1‑91 (MQ=255)
                                               aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGataat                                              >  2:3696684/1‑91 (MQ=255)
                                                 tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                            <  1:1521936/91‑1 (MQ=255)
                                                  atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                           <  1:93274/91‑1 (MQ=255)
                                                     taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                        <  1:3636465/91‑1 (MQ=255)
                                                     taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                        <  1:815851/91‑1 (MQ=255)
                                                       gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                      >  2:1839859/1‑91 (MQ=255)
                                                         gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGc                                    <  2:3821041/91‑1 (MQ=255)
                                                         gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATAATAATAAAAAAAAGc                                    <  2:767540/91‑1 (MQ=255)
                                                           tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                  >  1:1001048/1‑91 (MQ=255)
                                                             gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                <  1:1481924/91‑1 (MQ=255)
                                                             gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                >  1:3212239/1‑91 (MQ=255)
                                                             gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                >  1:3058752/1‑91 (MQ=255)
                                                             gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                <  1:1803400/91‑1 (MQ=255)
                                                              tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                               <  1:1982180/91‑1 (MQ=255)
                                                               cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                              >  2:1591554/1‑91 (MQ=255)
                                                                aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                             <  1:1398631/91‑1 (MQ=255)
                                                                 tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                            >  1:2974777/1‑91 (MQ=255)
                                                                   cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                          <  1:2862222/91‑1 (MQ=255)
                                                                    ccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                         >  2:634605/1‑91 (MQ=255)
                                                                      aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGATAATAAAAAAAAGCCTGTTGTCACAAt                       <  2:2883067/91‑1 (MQ=255)
                                                                       ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                      >  1:714582/1‑91 (MQ=255)
                                                                         cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                    >  2:1030736/1‑91 (MQ=255)
                                                                          acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAACAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                   >  2:3228872/1‑91 (MQ=255)
                                                                           caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                  <  1:1556345/91‑1 (MQ=255)
                                                                            aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                 >  1:1188798/1‑91 (MQ=255)
                                                                            aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                 >  1:2479701/1‑91 (MQ=255)
                                                                            aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                 >  1:3162596/1‑91 (MQ=255)
                                                                             gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGGGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                >  2:3541906/1‑91 (MQ=255)
                                                                             gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                >  2:3542832/1‑91 (MQ=255)
                                                                             gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                <  1:1952242/91‑1 (MQ=255)
                                                                             gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                >  1:3265126/1‑91 (MQ=255)
                                                                             gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                >  1:884439/1‑91 (MQ=255)
                                                                              tAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAg               >  1:1630185/1‑91 (MQ=255)
                                                                               aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa              >  2:1223765/1‑91 (MQ=255)
                                                                               aaCAATCCACACCATAGAAAAGAGTATATAACCTATTAAATAGTGAGGAAAAAGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa              <  2:2219415/91‑1 (MQ=255)
                                                                                aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac             <  1:1626344/91‑1 (MQ=255)
                                                                                aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac             >  2:430866/1‑91 (MQ=255)
                                                                                 tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca            >  2:1931054/5‑91 (MQ=17)
                                                                                  aacgcTCAACATTTAAAAGAAAAAATAACTTAATCAATACTAAAATAAATAAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac           <  2:2921694/87‑1 (MQ=255)
                                                                                  aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac           <  2:2077378/91‑1 (MQ=255)
                                                                                     ccTCACCATTGAAGAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAGAAAAAGCCTGTTGTCACGATGGTCATAGACACGAc        <  1:203268/91‑1 (MQ=255)
                                                                                     ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc        >  1:581222/1‑91 (MQ=255)
                                                                                      cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa       >  1:2099861/1‑91 (MQ=255)
                                                                                       tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACAt      >  2:1056376/1‑91 (MQ=255)
                                                                                       tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCACAGACACGACAt      >  1:3575164/1‑91 (MQ=255)
                                                                                        cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATa     >  1:187114/1‑91 (MQ=255)
                                                                                         aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc    >  1:2575861/1‑91 (MQ=255)
                                                                                         aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc    >  2:2312320/1‑91 (MQ=255)
                                                                                          ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt   >  2:1657931/1‑91 (MQ=255)
                                                                                          ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt   >  1:3464019/1‑91 (MQ=255)
                                                                                           cATTGAAGAGAGTATATAACCTTTTCAATAGTGTGGTATACGATAATCAAAAAAAGCCTGTTGTCATCATGGTCATCGACACGACATACtt  >  2:3675063/1‑91 (MQ=255)
                                                                                           |                                                                                          
CACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTT  >  USA300TCH1516_ALE/253883‑254063

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: