Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	431,505	T→C	10.9%	S77P (TCA→CCA)	SAS0387 →	exotoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	431,505	0	T	C	10.9%	309.3 / 5.6	91	S77P (TCA→CCA)	SAS0387	exotoxin
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base C (7/4);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAT  >  BX571857/431411‑431601                                                                                               |                                                                                                 aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCAGCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                <  2:523171/97‑1 (MQ=255)   aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaa                                                                                              <  1:730763/97‑1 (MQ=255)    tGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                             >  2:1078984/1‑97 (MQ=255)     gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACCCCACCTACAGCTAAAGTAGAAGCACCGCCATCAAAAc                                                                                            >  2:1308293/1‑97 (MQ=17)     gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTAGACGCACCGCCACCAAAAc                                                                                            >  2:942503/1‑97 (MQ=14)      caacaaCACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAAcc                                                                                           >  1:926945/1‑97 (MQ=255)       aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCa                                                                                          <  1:229746/97‑1 (MQ=255)         caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                        <  1:1258847/97‑1 (MQ=255)         caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCCACTAAAGTAGAAGCACCGCaaccaaaaccaac                                                                                        >  2:1395340/1‑96 (MQ=17)          aacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGACGCACCACCACCAAAACCaacc                                                                                       >  2:647920/1‑95 (MQ=14)          aacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAc                                                                                       >  1:454381/1‑97 (MQ=255)          aacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAc                                                                                       <  2:580251/97‑1 (MQ=255)           acacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcg                                                                                      <  1:265284/97‑1 (MQ=255)           acacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcg                                                                                      >  1:78743/1‑97 (MQ=255)            cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCCATCACAACCacacac                                                                                     >  2:596722/1‑95 (MQ=18)            cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgc                                                                                     <  1:1163540/97‑1 (MQ=255)              ccATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGa                                                                                   >  1:322746/1‑97 (MQ=255)               cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGac                                                                                  <  1:715617/97‑1 (MQ=255)               cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGac                                                                                  <  2:843799/97‑1 (MQ=255)                 tCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaa                                                                                <  2:1398349/97‑1 (MQ=255)                  cTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaac                                                                               >  1:1175598/1‑97 (MQ=255)                   ttCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTCCAACTAAAGTAGCAGCACCGCCACCAAAACCAAACGCGacaaca                                                                              >  2:302219/1‑97 (MQ=18)                    tCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAacac                                                                             <  1:898615/97‑1 (MQ=255)                    tCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAacac                                                                             <  2:45725/97‑1 (MQ=255)                      aaCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCa                                                                           <  1:400161/97‑1 (MQ=255)                      aaCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCa                                                                           <  2:1424721/97‑1 (MQ=255)                       aCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCAt                                                                          <  1:668927/97‑1 (MQ=255)                        gtAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                         <  2:335427/96‑1 (MQ=255)                        cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                         >  2:1295567/1‑97 (MQ=255)                           aaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCACCACCAAAACCAAACGCGACAACACCATCtcc                                                                      >  2:943915/1‑95 (MQ=21)                           aaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTAGAAACACCACAACCACAACCAACCGCGACAACACCATcacc                                                                      >  2:808081/1‑94 (MQ=17)                             gTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCaa                                                                    >  1:626839/1‑97 (MQ=255)                             gTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACAATCTTCaa                                                                    >  1:322456/1‑97 (MQ=255)                             gTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCACCTTCaa                                                                    >  2:105814/1‑97 (MQ=255)                              tGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAAc                                                                   >  1:48936/1‑97 (MQ=255)                               ggAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACt                                                                  >  2:127432/1‑97 (MQ=255)                               ggAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAACCACCGCAACCACACCCCAACGCGACAACACCATCTCCAAca                                                                  >  2:620545/1‑96 (MQ=18)                                gAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAACACCAAACGCGACAACACCACCTTCAACtc                                                                 >  2:832038/1‑96 (MQ=255)                                gAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTa                                                                 >  1:84036/1‑97 (MQ=255)                                 aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaa                                                                <  2:586689/97‑1 (MQ=255)                                  aGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaaa                                                               >  1:809532/1‑97 (MQ=255)                                     aCCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTg                                                            >  1:104771/1‑97 (MQ=255)                                       cGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCCACAACACCATCTTCAACAAAAATGGa                                                          >  1:1125206/1‑97 (MQ=255)                                        gCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGaa                                                         >  1:1383912/1‑97 (MQ=255)                                         caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAg                                                        <  1:312028/97‑1 (MQ=255)                                          aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                       <  1:1193584/97‑1 (MQ=255)                                          aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                       <  1:906009/97‑1 (MQ=255)                                           acaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCa                                                      >  2:484960/1‑97 (MQ=255)                                             aaGCAGCAAACGCGACAACACCATCTTCAACTAACGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAcc                                                    >  2:1094274/1‑97 (MQ=255)                                             aaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAcc                                                    >  1:1261765/1‑97 (MQ=255)                                              agcagcAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCg                                                   >  1:974267/1‑97 (MQ=255)                                                  gcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaac                                               <  1:358111/97‑1 (MQ=255)                                                  gcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaac                                               <  1:109887/97‑1 (MQ=255)                                                    aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaacaa                                             >  2:278910/1‑97 (MQ=255)                                                    aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaacaa                                             <  1:236319/97‑1 (MQ=255)                                                      aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagc                                           <  1:677867/97‑1 (MQ=255)                                                      aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagc                                           <  2:242544/97‑1 (MQ=255)                                                      aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagc                                           <  1:430286/97‑1 (MQ=255)                                                       cgcgACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                          >  1:1095460/1‑97 (MQ=255)                                                         cgACAACACCATCTTCAACTAAAGTAGAAGCGCCGCAATCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagcagc                                        <  2:141741/97‑1 (MQ=255)                                                          gACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAgcagca                                       >  2:1186193/1‑97 (MQ=255)                                                          gACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAgcagca                                       <  2:224796/97‑1 (MQ=255)                                                           acaacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCaa                                      <  2:917956/97‑1 (MQ=255)                                                            caacaCCACCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCaaa                                     <  2:195780/97‑1 (MQ=255)                                                               cacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgc                                  <  1:646896/97‑1 (MQ=255)                                                                acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                 >  2:574760/1‑97 (MQ=255)                                                                    ccTTCAACTAAAGTAGAAGCACCGCAACCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaa                             <  2:639521/96‑1 (MQ=255)                                                                     cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaac                            <  1:1318519/97‑1 (MQ=255)                                                                      ttCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaaca                           >  2:279483/1‑97 (MQ=255)                                                                       tCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAacac                          <  1:575193/97‑1 (MQ=255)                                                                           cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                      >  1:1220667/1‑97 (MQ=255)                                                                            tAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                     <  1:924861/97‑1 (MQ=255)                                                                            tAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                     <  1:762847/97‑1 (MQ=255)                                                                             aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCtt                    >  1:1215986/1‑97 (MQ=255)                                                                              aaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTc                   <  1:1138214/97‑1 (MQ=255)                                                                               aGTAGAAGCCCCGCAACCAAAACCGAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCg                  <  2:999936/97‑1 (MQ=255)                                                                               aGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCg                  <  1:939399/97‑1 (MQ=255)                                                                                gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                 >  2:249741/1‑97 (MQ=255)                                                                                gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                 >  2:249161/1‑97 (MQ=255)                                                                                 tAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                >  1:1216530/1‑97 (MQ=255)                                                                                 tAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                >  1:187814/1‑97 (MQ=255)                                                                                    aaGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTaa             >  1:465184/1‑97 (MQ=255)                                                                                     aGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt            >  2:25754/1‑97 (MQ=255)                                                                                     aGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt            >  2:36676/1‑97 (MQ=255)                                                                                      gCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATg           >  2:679152/1‑97 (MQ=255)                                                                                       cACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGt          >  2:1050820/1‑97 (MQ=255)                                                                                        aCCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         <  2:1153805/97‑1 (MQ=255)                                                                                        aCCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         >  2:917009/1‑97 (MQ=255)                                                                                        aCCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa         >  2:560944/1‑97 (MQ=255)                                                                                         ccGCAATCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAg        <  2:848697/97‑1 (MQ=255)                                                                                         ccGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAg        <  1:461232/97‑1 (MQ=255)                                                                                            cAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGaca     <  1:162118/97‑1 (MQ=255)                                                                                               tCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt  >  1:13658/1‑97 (MQ=255)                                                                                               |                                                                                                 AAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAT  >  BX571857/431411‑431601

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