Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,119 C→A 5.2% G355C (GGT→TGT)  patA_2 ← Putative N‑acetyl‑LL‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,055,1190CA 5.2% 951.8 / 26.1 270G355C (GGT→TGT) patA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (143/113);  new base A (2/12);  total (145/125)
Fisher's exact test for biased strand distribution p-value = 4.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 276 total aligned reads displayed.

ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATAT  >  USA300TCH1516_ALE/1055029‑1055208
                                                                                          |                                                                                         
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAACACATCCATTCCCTCAGCCTAGGAAATCCTTACAAATCCTTTGcc                                                                                           >  2:1697089/1‑91 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                           <  2:3575418/91‑1 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATACATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                           >  1:447605/1‑91 (MQ=255)
 tttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCg                                                                                          >  2:1600122/1‑91 (MQ=255)
   taGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGaa                                                                                        >  1:3624170/1‑91 (MQ=255)
   taGTAATTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGaa                                                                                        <  2:3651871/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTc                                                                                     <  2:2298820/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTc                                                                                     <  1:1017211/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTc                                                                                     >  1:2511469/1‑91 (MQ=255)
        tttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTg                                                                                   >  2:1974776/1‑91 (MQ=255)
         ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTcggt                                                                                  <  1:646510/91‑4 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTa                                                                                 <  1:1877923/91‑1 (MQ=255)
           aaaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaa                                                                                >  1:403986/1‑91 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                              <  2:1957679/91‑1 (MQ=255)
               atTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGa                                                                            <  2:2706706/91‑1 (MQ=255)
                tttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGAt                                                                           >  2:1913532/1‑91 (MQ=255)
                tttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGAt                                                                           <  1:243373/91‑1 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                          >  2:644224/1‑91 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                          >  1:792687/1‑91 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAAGCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                          <  2:2230883/91‑1 (MQ=255)
                  tAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGa                                                                         <  1:3347097/91‑1 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                      >  2:1579576/1‑91 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                      >  1:1066431/1‑91 (MQ=255)
                       aGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTg                                                                    >  2:139134/1‑91 (MQ=255)
                       aGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTg                                                                    <  1:1668999/91‑1 (MQ=255)
                        gTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTgg                                                                   >  1:678304/1‑91 (MQ=255)
                         tctcTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGa                                                                  <  1:2949032/91‑1 (MQ=255)
                           tctTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAAc                                                                <  2:2564081/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAGTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                              <  2:484928/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                              >  1:2177101/1‑91 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                              >  1:861276/1‑91 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                              >  1:2017657/1‑91 (MQ=255)
                              tCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAt                                                             >  2:1005474/1‑91 (MQ=255)
                               cATACCTGCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAtt                                                            <  1:1965966/91‑1 (MQ=255)
                                aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACACATCCTATACCGAATTCTGTAAAGGATGAGCCTGGAACGAttt                                                           >  2:3915384/1‑90 (MQ=255)
                                aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTg                                                           >  1:758965/1‑91 (MQ=255)
                                 tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATGCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                          <  2:31013/91‑1 (MQ=255)
                                 tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                          <  2:2389802/91‑1 (MQ=255)
                                 tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                          <  1:1289644/91‑1 (MQ=255)
                                   ccTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTa                                                        >  2:2804016/1‑91 (MQ=255)
                                    cttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACACATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTaa                                                       >  2:755573/1‑91 (MQ=255)
                                    cttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTaa                                                       <  1:1367647/91‑1 (MQ=255)
                                     ttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTaaa                                                      >  1:2249241/1‑91 (MQ=255)
                                      tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                     >  1:2775685/1‑91 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                    >  1:568428/1‑91 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                    <  2:1845139/91‑1 (MQ=255)
                                         tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                  <  2:3747706/91‑1 (MQ=255)
                                         tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCAAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                  <  2:125022/91‑1 (MQ=255)
                                          ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                 <  2:3805457/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                >  1:306495/1‑91 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                <  2:1657749/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAACTCTGTAAAGGATGAGCCTGGAACGATCGCTAAATGTGtt                                                >  2:3805732/1‑91 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCATTTACCGAATTCTGTAAAGGATGAGTCTGGAACGATTACTAAATGTGtt                                                >  2:3011237/1‑91 (MQ=255)
                                            aaaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTg                                               >  2:2617463/1‑91 (MQ=255)
                                             aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                              >  1:3519242/1‑91 (MQ=255)
                                               aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                            <  1:18543/91‑1 (MQ=255)
                                               aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                            >  2:1143606/1‑91 (MQ=255)
                                               aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                            >  2:3589588/1‑91 (MQ=255)
                                                cgtCCATTTCGTAAGCATAGGAAAGCGTTACAAATCGTTTACCGAATTCTGGAAAGGATGAGCCTGGAACGATTGCTAAATGTGTCGAttt                                           <  2:1893339/89‑2 (MQ=255)
                                                cATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTc                                           <  2:820202/91‑1 (MQ=255)
                                                cATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTc                                           >  1:1259877/1‑91 (MQ=255)
                                                cATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTc                                           >  1:2655487/1‑91 (MQ=255)
                                                 aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                          >  2:1677791/1‑91 (MQ=255)
                                                  tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                         >  1:151186/1‑91 (MQ=255)
                                                  tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATCCaa                                         >  1:3556584/1‑91 (MQ=255)
                                                  tCCATTTCATAAGCATAGGAAATCCGTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                         <  2:3421261/91‑1 (MQ=255)
                                                   ccATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAg                                        <  1:326009/91‑1 (MQ=255)
                                                     aTTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTa                                      >  1:61087/1‑91 (MQ=255)
                                                     aTTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTa                                      >  1:1703081/1‑91 (MQ=255)
                                                      tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTaa                                     >  2:383328/1‑91 (MQ=255)
                                                          aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                 >  2:1209798/1‑91 (MQ=255)
                                                          aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                 >  2:3304029/1‑91 (MQ=255)
                                                            aaGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGac                               >  2:2724645/1‑91 (MQ=255)
                                                                 tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                          >  2:1353378/1‑91 (MQ=255)
                                                                 tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                          <  2:3352383/91‑1 (MQ=255)
                                                                 tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                          >  1:812578/1‑91 (MQ=255)
                                                                  aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                         <  2:1535934/91‑1 (MQ=255)
                                                                   ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                        <  2:3204736/91‑1 (MQ=255)
                                                                        tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                   >  2:1163481/1‑91 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACATAAAATTCGaaa                  >  2:2864337/1‑91 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                  >  1:2581551/1‑91 (MQ=255)
                                                                          cTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAt                 >  1:3177447/1‑91 (MQ=255)
                                                                           ttACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATc                >  2:2519169/1‑91 (MQ=255)
                                                                            tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg               >  1:1627159/1‑91 (MQ=255)
                                                                            tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg               >  1:3840346/1‑91 (MQ=255)
                                                                             aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCATGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGt              >  1:1554918/1‑91 (MQ=255)
                                                                              cAAGTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAGATGTGTTGATTCAAGTAAATCGACACAAATTTCGACATCGTc             <  2:2630449/91‑1 (MQ=255)
                                                                              cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTc             >  1:2106512/1‑91 (MQ=255)
                                                                              cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTc             <  1:304409/91‑1 (MQ=255)
                                                                               aaaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTATATCGACACAAAATTCGAAATCGTCa            >  1:3684535/1‑91 (MQ=255)
                                                                                aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt           <  2:6180/91‑1 (MQ=255)
                                                                                  tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCg         >  1:3431599/1‑91 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAAGCGACACAAAATTCGAAATCGTCATCgg        <  1:349293/91‑1 (MQ=255)
                                                                                    cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt       >  1:2682520/1‑91 (MQ=255)
                                                                                    cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt       <  1:2871953/91‑1 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa     <  1:3861613/91‑1 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa     >  1:513517/1‑91 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa     >  1:1534043/1‑91 (MQ=255)
                                                                                      ttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa     >  1:1770783/1‑91 (MQ=255)
                                                                                       tACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGat    <  1:377544/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATAT  >  USA300TCH1516_ALE/1055029‑1055208

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: