Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 973,213 C→A E175* (GAG→TAG)  yitU_1 ← 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE973,2130CA97.5% 705.1 / 6.3 198E175* (GAG→TAG) yitU_15‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (85/108);  total (87/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 198 total aligned reads displayed.

GTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGAT  >  USA300TCH1516_ALE/973126‑973301
                                                                                       |                                                                                        
gTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAat                                                                                       >  2:1782918/1‑91 (MQ=255)
gTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAat                                                                                       <  1:2844745/91‑1 (MQ=255)
gTCTAACTTGCTCAATGACTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAat                                                                                       >  2:855339/1‑91 (MQ=255)
 tCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAata                                                                                      >  2:2333863/1‑91 (MQ=255)
   tAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATg                                                                                    >  2:445796/1‑91 (MQ=255)
       ttGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCg                                                                                <  2:2988542/91‑1 (MQ=255)
        tGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCgg                                                                               <  2:2549728/91‑1 (MQ=255)
        tGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCgg                                                                               <  2:127541/91‑1 (MQ=255)
         gCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCTAGCGTCGATGCTAAATATGATCGGc                                                                              >  1:2894076/1‑91 (MQ=255)
         gCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGc                                                                              <  1:1408372/91‑1 (MQ=255)
          cTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCa                                                                             <  2:2067096/91‑1 (MQ=255)
          cTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCa                                                                             >  2:2640931/1‑91 (MQ=255)
          cTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCa                                                                             <  1:2830350/91‑1 (MQ=255)
              aTGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaa                                                                         >  2:197381/1‑91 (MQ=255)
              aTGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaa                                                                         <  1:345915/91‑1 (MQ=255)
              aTGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaa                                                                         >  1:1875201/1‑91 (MQ=255)
               tGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCATCGTCGATGCTAAATATGATCGGCATaaaa                                                                        >  1:908166/1‑91 (MQ=255)
                gCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaaaa                                                                       <  2:2825427/91‑1 (MQ=255)
                  ctctTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATg                                                                     >  1:1984425/1‑91 (MQ=255)
                  ctctTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATg                                                                     >  1:1200958/1‑91 (MQ=255)
                    ctTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAg                                                                   >  2:1297072/1‑91 (MQ=255)
                    ctTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAg                                                                   <  1:564750/91‑1 (MQ=255)
                     ttGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGt                                                                  >  1:2009943/1‑91 (MQ=255)
                      tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCAACATAAAAATGAGTa                                                                 <  2:2976154/91‑1 (MQ=255)
                       gCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTaa                                                                <  1:252095/91‑1 (MQ=255)
                       gCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTaa                                                                <  1:2737205/91‑1 (MQ=255)
                        cTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAg                                                               <  2:1805613/91‑1 (MQ=255)
                          ttattaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCa                                                             >  2:1543/1‑91 (MQ=255)
                           tattaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCat                                                            <  1:2296576/91‑1 (MQ=255)
                             ttaATACCAAGTTTTATAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCatat                                                          <  1:1093573/91‑1 (MQ=255)
                             ttaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCatat                                                          <  2:1124784/91‑1 (MQ=255)
                              taATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtt                                                         <  1:2180160/91‑1 (MQ=255)
                              taATACCAAGTTTTACAATTTCAATGACAGGGAATGATGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtt                                                         <  2:760134/91‑1 (MQ=255)
                              taATACCAAGTGTTACAATTTCACTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAATTATGATCGGCATAAAAATGAGTAAGCATAtt                                                         <  2:2361802/91‑1 (MQ=255)
                               aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGTCCCAGCGTCGATGCTAAATATGATCGGCATAAAAGTGAGTAAGCATAttt                                                        <  1:79958/91‑1 (MQ=255)
                               aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttt                                                        <  2:2345995/91‑1 (MQ=255)
                               aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttt                                                        <  2:2760869/91‑1 (MQ=255)
                               aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttt                                                        <  2:2386940/91‑1 (MQ=255)
                               aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttt                                                        <  2:1250570/91‑1 (MQ=255)
                                aTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtttt                                                       <  1:412061/91‑1 (MQ=255)
                                aTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtttt                                                       <  1:164442/91‑1 (MQ=255)
                                aTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtttt                                                       <  2:634245/91‑1 (MQ=255)
                                 tACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttttt                                                      >  2:554730/1‑91 (MQ=255)
                                 tACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAttttt                                                      >  1:1506487/1‑91 (MQ=255)
                                  aCCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTg                                                     >  1:564412/1‑91 (MQ=255)
                                      aGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAttt                                                 >  2:482498/1‑91 (MQ=255)
                                       gTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTc                                                <  2:2465368/91‑1 (MQ=255)
                                       gTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTc                                                <  1:2922771/91‑1 (MQ=255)
                                       gTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTc                                                <  2:1916182/91‑1 (MQ=255)
                                        ttttACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCa                                               >  1:639616/1‑91 (MQ=255)
                                            aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTa                                           >  1:954919/1‑91 (MQ=255)
                                            aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTa                                           >  2:2944943/1‑91 (MQ=255)
                                             cAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAt                                          >  1:50122/1‑91 (MQ=255)
                                               aTTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAttt                                        <  2:1031091/91‑1 (MQ=255)
                                                tttCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAGTGAGTAAGCATATTTTTGATTTCAGGTAtttt                                       <  2:3040726/91‑1 (MQ=255)
                                                tttCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGGGTAAGCATATTTTTGATTTCAGGTAtttt                                       <  2:3056616/91‑1 (MQ=255)
                                                 ttCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTa                                      <  2:349282/91‑1 (MQ=255)
                                                  tCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTAc                                     <  1:878479/91‑1 (MQ=255)
                                                   cAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACt                                    <  2:917301/91‑1 (MQ=255)
                                                   cAATAACAGAAAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCTGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACt                                    <  1:1163243/91‑1 (MQ=255)
                                                    aaTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                   >  1:504523/1‑91 (MQ=255)
                                                     aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACttt                                  >  2:2184411/1‑91 (MQ=255)
                                                      tGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttc                                 <  2:1584891/91‑1 (MQ=255)
                                                       gACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttct                                <  2:1598727/91‑1 (MQ=255)
                                                        aCAGGGAATGGTGCGCCCCCGCGTCTATGCCAAATATGATCGGCATAACAATGAGTCAGCACATTTTTGGTGTCAGGTATCTTACTttcct                               >  2:1371383/1‑89 (MQ=255)
                                                        aCAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttctt                               <  2:1236019/91‑1 (MQ=255)
                                                        aCAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttctt                               <  1:1187562/91‑1 (MQ=255)
                                                         cAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttcttc                              >  2:2100247/1‑91 (MQ=255)
                                                         cAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttcttc                              >  2:1618531/1‑91 (MQ=255)
                                                             gAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                          >  1:481507/1‑91 (MQ=255)
                                                              aaTGGTGCGCCCCAGCGTGGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCtt                         <  2:1225457/91‑1 (MQ=255)
                                                              aaTGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCtt                         <  1:336908/91‑1 (MQ=255)
                                                               aTGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTc                        >  1:236719/1‑91 (MQ=255)
                                                                tGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCa                       <  1:62305/91‑1 (MQ=255)
                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                     <  1:1112143/91‑1 (MQ=255)
                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                     <  1:1260557/91‑1 (MQ=255)
                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                     <  1:1286874/91‑1 (MQ=255)
                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                     <  2:2974644/91‑1 (MQ=255)
                                                                   tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                    >  2:2925859/1‑91 (MQ=255)
                                                                   tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                    >  2:1136949/1‑91 (MQ=255)
                                                                   tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                    <  1:2335132/91‑1 (MQ=255)
                                                                   tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                    >  1:2095160/1‑91 (MQ=255)
                                                                   tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                    <  1:1921162/91‑1 (MQ=255)
                                                                     cgcCCCAGCGTCGATGTTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaa                  <  2:2374935/91‑1 (MQ=255)
                                                                     cgcCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaa                  >  1:986045/1‑91 (MQ=255)
                                                                     cgcCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaa                  >  2:726206/1‑91 (MQ=255)
                                                                     cgcCCCAGCGTCCATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaa                  >  1:2385460/1‑91 (MQ=255)
                                                                       ccccAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaaa                >  2:1895664/1‑91 (MQ=255)
                                                                        cccAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAt               <  1:1291923/91‑1 (MQ=255)
                                                                         ccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAtt              <  1:380576/91‑1 (MQ=255)
                                                                            gCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAg           >  2:1889778/1‑91 (MQ=255)
                                                                             cGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAgg          <  2:394808/91‑1 (MQ=255)
                                                                               tCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTa        <  1:257622/91‑1 (MQ=255)
                                                                               tCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTa        >  2:804223/1‑91 (MQ=255)
                                                                                cGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAg       <  2:2385710/91‑1 (MQ=255)
                                                                                cGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAg       >  2:1939773/1‑91 (MQ=255)
                                                                                 gATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg      >  2:590380/1‑91 (MQ=255)
                                                                                 gATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg      <  2:2520343/91‑1 (MQ=255)
                                                                                  aTGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg     <  1:2931420/91‑1 (MQ=255)
                                                                                     cTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGAt  <  1:1928261/91‑1 (MQ=255)
                                                                                       |                                                                                        
GTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGAT  >  USA300TCH1516_ALE/973126‑973301

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: