Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,119 C→A 5.6% 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.6% 796.8 / 28.7 231G355C (GGT→TGT) patA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (123/95);  new base A (6/7);  total (129/102)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 241 total aligned reads displayed.

ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATG  >  USA300TCH1516_ALE/1055029‑1055209
                                                                                          |                                                                                          
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                            <  2:1622947/91‑1 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                            >  2:2595244/1‑91 (MQ=255)
 tttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCg                                                                                           >  2:808009/1‑91 (MQ=255)
    aGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAAt                                                                                        <  2:660568/91‑1 (MQ=255)
     gTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAAtt                                                                                       >  1:3289274/1‑91 (MQ=255)
       aTTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCt                                                                                     >  2:3131705/1‑91 (MQ=255)
       aTTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCt                                                                                     <  1:2014482/91‑1 (MQ=255)
       aTTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTcg                                                                                     <  2:2807569/91‑2 (MQ=255)
        tttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTg                                                                                    >  2:3334857/1‑91 (MQ=255)
         ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGt                                                                                   <  1:3238911/91‑1 (MQ=255)
         ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTcggt                                                                                   <  2:1583049/91‑4 (MQ=255)
         ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTcggt                                                                                   <  2:1092375/91‑4 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                  >  1:1619903/1‑91 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCAGTa                                                                                  >  1:2388231/1‑91 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTATGTa                                                                                  >  1:3237658/1‑91 (MQ=255)
            aaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaaa                                                                                >  2:1680132/1‑91 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                               <  1:1010195/91‑1 (MQ=255)
               atTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGa                                                                             >  1:3314883/1‑91 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                           >  2:2729502/1‑91 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                           <  1:3384104/91‑1 (MQ=255)
                    gCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGc                                                                        >  2:2745490/1‑91 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                              >  2:2130043/1‑84 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                       >  2:24911/1‑91 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                       >  2:24154/1‑91 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                       <  1:1433357/91‑1 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                       >  1:1066567/1‑91 (MQ=255)
                      aaGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                      <  1:2824141/91‑1 (MQ=255)
                      aaGTCTCTTCATACCTTCTTGTAAACCATCCATTTCATAAGCATTGGAGATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                      <  2:3486452/91‑1 (MQ=255)
                        gTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTgg                                                                    <  2:2815195/91‑1 (MQ=255)
                         tctcTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGa                                                                   <  2:1778635/91‑1 (MQ=255)
                          ctctTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGaa                                                                  >  2:765322/1‑91 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                               <  2:3771706/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                               <  2:1082648/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                               <  1:591837/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                               >  1:117914/1‑91 (MQ=255)
                              tCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAt                                                              >  1:2349098/1‑91 (MQ=255)
                               cATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAtt                                                             >  2:529576/1‑91 (MQ=255)
                                aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTg                                                            <  1:2220731/91‑1 (MQ=255)
                                  aCCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCt                                                          >  2:30642/1‑91 (MQ=255)
                                      tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                      <  1:1109292/91‑1 (MQ=255)
                                      tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCATGGAAAGATTGCTAAAt                                                      >  2:1615555/1‑91 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                     <  1:1345766/91‑1 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                     <  2:3359302/91‑1 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                     <  2:546082/91‑1 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                     <  1:983866/91‑1 (MQ=255)
                                          ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                  <  1:2489549/91‑1 (MQ=255)
                                          ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                  <  2:2080215/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                 <  1:3852256/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                 >  2:2931027/1‑91 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                 <  1:2542074/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                 <  1:2699600/91‑1 (MQ=255)
                                            aaaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTg                                                >  1:726826/1‑91 (MQ=255)
                                             aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                               >  1:1753874/1‑91 (MQ=255)
                                             aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                               >  1:3531463/1‑91 (MQ=255)
                                              aaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAt                                              >  2:718814/1‑91 (MQ=255)
                                              aaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAt                                              >  2:2373427/1‑91 (MQ=255)
                                                 aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                           >  1:851520/1‑91 (MQ=255)
                                                 aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                           <  2:13914/91‑1 (MQ=255)
                                                  tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                      >  1:2377827/1‑79 (MQ=255)
                                                    cATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGCATTCTGCAAAGGATGAGCCTGGAACCATTGCTAAACGTGTTGATTCAAGt                                        >  2:3603029/1‑91 (MQ=255)
                                                      tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                  >  2:1421151/1‑79 (MQ=255)
                                                      tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                  <  1:1421151/79‑1 (MQ=255)
                                                          aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                  <  1:321210/91‑1 (MQ=255)
                                                          aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                  >  1:2897882/1‑91 (MQ=255)
                                                           tAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGa                                 >  2:1467039/1‑91 (MQ=255)
                                                            aaGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGac                                <  2:2181873/91‑1 (MQ=255)
                                                              gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                              >  2:2042762/1‑91 (MQ=255)
                                                              gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                              >  1:411589/1‑91 (MQ=255)
                                                               cATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGAcaca                             <  1:1068882/91‑1 (MQ=255)
                                                                aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                            <  2:2326852/91‑1 (MQ=255)
                                                                 tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                           <  2:3233795/91‑1 (MQ=255)
                                                                    gAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAtt                        <  1:1563448/91‑1 (MQ=255)
                                                                      aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                      >  2:329050/1‑91 (MQ=255)
                                                                       aTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGa                     <  2:1563777/91‑1 (MQ=255)
                                                                        tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                    <  1:2865800/91‑1 (MQ=255)
                                                                        tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                    <  1:328987/91‑1 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                   <  1:1009551/91‑1 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                   <  1:2132165/91‑1 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                   <  1:1680132/91‑1 (MQ=255)
                                                                          cTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACTATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAt                  <  1:3086474/91‑1 (MQ=255)
                                                                          cTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAt                  <  2:2976403/91‑1 (MQ=255)
                                                                            tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg                <  2:437196/91‑1 (MQ=255)
                                                                             aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGt               >  2:2677752/1‑91 (MQ=255)
                                                                             aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGt               >  2:3693456/1‑91 (MQ=255)
                                                                              cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTc              >  1:499250/1‑91 (MQ=255)
                                                                               tcaTCCTTTACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCa             >  1:3427192/3‑91 (MQ=38)
                                                                                aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt            >  2:1974133/1‑91 (MQ=255)
                                                                                  tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATGGTCATCg          >  1:1947933/1‑91 (MQ=255)
                                                                                   ccTTTACCGAATTCTGAAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg         <  2:3702886/91‑1 (MQ=255)
                                                                                     tttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg       >  2:2268864/1‑91 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa      <  2:3801199/91‑1 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa      <  1:2674537/91‑1 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa      <  1:105544/91‑1 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa      <  1:2102200/91‑1 (MQ=255)
                                                                                      ttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa      >  1:1518938/1‑91 (MQ=255)
                                                                                       tACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGat     >  1:334521/1‑91 (MQ=255)
                                                                                        aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata    >  1:3733400/1‑91 (MQ=255)
                                                                                        aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata    >  1:3863406/1‑91 (MQ=255)
                                                                                         ccGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGatat   <  2:1177118/91‑1 (MQ=255)
                                                                                          cGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg  >  1:2398321/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATG  >  USA300TCH1516_ALE/1055029‑1055209

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: