BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	115,788	A→G	C438R (TGC→CGC)	hofB ←	HofB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,788	0	A	G	98.0%	449.7 / ‑0.8	147	C438R (TGC→CGC)	hofB	HofB
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (80/64); total (83/64)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 177 total aligned reads displayed.

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

ccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAggcggc                                                                                         >  1:1152145/1‑93 (MQ=255)

 caGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAggcggcg                                                                                        >  1:1284029/1‑92 (MQ=255)

 caGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGcc                                                                                        >  2:4179646/1‑93 (MQ=255)

      ccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtttt                                                                                   <  1:1564299/93‑1 (MQ=255)

      ccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtttt                                                                                   <  2:631907/93‑1 (MQ=255)

           gTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaa                                                                              <  2:911029/93‑1 (MQ=255)

           gTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaa                                                                              >  2:771676/1‑93 (MQ=255)

            ttGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  2:62485/93‑1 (MQ=255)

             tGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAg                                                                            <  1:2404730/93‑1 (MQ=255)

                ttACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGt                                                                         >  2:3921796/1‑93 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGtc                                                                        >  2:2186946/1‑92 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  1:387097/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  1:3165920/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCAAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  2:8164/93‑1 (MQ=255)

                  tcTCGCCATGCGGCATCCCCAGTACGCTGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTAc                                                                       <  1:3011216/92‑1 (MQ=255)

                  aCTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTAc                                                                       >  2:3730443/1‑93 (MQ=255)

                     cGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCa                                                                    <  1:4717928/93‑1 (MQ=255)

                      gCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCAt                                                                   <  1:1102286/93‑1 (MQ=255)

                       tcATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTTAGCCTTGCTCTACGGCCAGGCGGCCGTCTGCTAAAATCGTACGCATa                                                                  >  2:4216285/2‑93 (MQ=255)

                       ccATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATa                                                                  >  2:4309936/1‑93 (MQ=255)

                       ccATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATa                                                                  <  2:262484/93‑1 (MQ=255)

                       ccATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATa                                                                  >  1:3088200/1‑93 (MQ=255)

                        cATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAc                                                                 >  1:2790172/1‑93 (MQ=255)

                         aTGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAcc                                                                <  1:4011708/93‑1 (MQ=255)

                          tGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAccc                                                               <  2:2655841/93‑1 (MQ=255)

                             ggCATCCCGAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                            <  1:1751552/93‑1 (MQ=255)

                             ggCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                            <  2:4608273/93‑1 (MQ=255)

                             ggCACCCCCAGGTCCCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGCCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                            <  2:4462470/93‑1 (MQ=255)

                                aTCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGt                                                         >  1:3573915/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        <  2:4611630/93‑1 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        <  2:1626439/93‑1 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  1:3742762/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGGACGCATACCCGCCTGTc                                                        >  1:4218061/1‑93 (MQ=255)

                                   cccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgt                                                      >  2:1668767/1‑93 (MQ=255)

                                   cccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgt                                                      >  1:2001128/1‑93 (MQ=255)

                                   cccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgt                                                      <  1:2902089/93‑1 (MQ=255)

                                    ccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtg                                                     >  2:2149165/1‑93 (MQ=255)

                                      aGTACGCGGATAAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcg                                                   <  2:2716200/93‑1 (MQ=255)

                                           gcgGATTAACTCTTCAAGGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTgcgtgcgt                                              <  2:1197256/93‑1 (MQ=255)

                                           gcgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTgcgtgcgt                                              <  2:4133397/93‑1 (MQ=255)

                                             ggATTAACTCTTCCAAGGTGGTTAAGCCTTGCTCTACGGGCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCACCTGTCGTGCGTGCGttt                                            >  2:3147409/1‑93 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            <  1:3356454/93‑1 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  1:547942/1‑93 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  2:1205226/1‑93 (MQ=255)

                                              gATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                           >  1:818978/1‑93 (MQ=255)

                                                ttACCTCTTCAAAGGTGGTAAAGCCTTGCTCTGCGCCCAGGCGGCCGTTTTCAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCa                                         <  2:188837/93‑1 (MQ=255)

                                                  aaCTCTTCAAAGGTGGTTAAGCCTTGCTCTCCGGCCAGGCGGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                       <  2:343115/93‑1 (MQ=255)

                                                  aaCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                       >  2:3698646/1‑93 (MQ=255)

                                                    ctctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGa                                     <  1:2665538/93‑1 (MQ=255)

                                                     tctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    <  2:3576237/93‑1 (MQ=255)

                                                     tctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    >  2:2925503/1‑93 (MQ=255)

                                                     tctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    >  2:1610369/1‑93 (MQ=255)

                                                     ctttCAAAGGTGGTTAAGCCTTGCTCTCCGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    <  2:2804272/91‑1 (MQ=255)

                                                      ctTCAAAGGTGGTTAAGCCTTGCTCTACGGCAAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAtt                                   <  2:1384081/93‑1 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCTTTGCTCTACGGCCAGGCGGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  <  2:1841445/93‑1 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  >  1:2152932/1‑93 (MQ=255)

                                                        tCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCa                                 >  1:4610501/1‑93 (MQ=255)

                                                         cAAGGGTGTTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCaa                                <  2:399994/93‑1 (MQ=255)

                                                          aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                               <  1:2505549/93‑1 (MQ=255)

                                                           aaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGACTGTCGTGCGTGCGTTTCCAGCGACTCAACg                              >  2:2741375/1‑93 (MQ=255)

                                                            aGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGt                             >  2:2219161/1‑93 (MQ=255)

                                                            aGGTGGTTAAGCCTTGATCTACGGCCAGGCAGCCGTTTAAAAAAAGCGAACGCATACCCGCCTGTCGTGAGTGCGTTTCCAGCGATTCAACGt                             <  1:4211970/93‑1 (MQ=255)

                                                             ggtggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTc                            <  1:4380627/93‑1 (MQ=255)

                                                              gtggtTAGCCCTTCCCCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                           <  1:1589652/93‑1 (MQ=255)

                                                                ggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                         >  1:618707/1‑93 (MQ=255)

                                                                 gtTAAGCCTTGCTCTACGGCCAGGCGCCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:2088110/93‑1 (MQ=255)

                                                                  ttAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAt                       >  1:1622338/1‑93 (MQ=255)

                                                                    aacccTTGCTCTACGGCCAGGCGGCCTTTTTCAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                     <  2:4741960/90‑1 (MQ=255)

                                                                    aaGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                     >  2:2551505/1‑93 (MQ=255)

                                                                    aaGCCTTGCTCTACGGCCAGGCAGACGTTTTCAAAAAGCGAACGCATAACCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGGAttg                     >  2:2696879/1‑92 (MQ=255)

                                                                     aGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAg                    >  1:3231654/1‑93 (MQ=255)

                                                                      gCCTTGCTCTCCGGCCAGGGGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                   <  2:1245811/93‑1 (MQ=255)

                                                                      gCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                   >  2:4276809/1‑93 (MQ=255)

                                                                      gCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                   <  1:4631526/93‑1 (MQ=255)

                                                                      gCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGc                   >  2:2398673/1‑93 (MQ=255)

                                                                       ccTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCg                  <  1:896393/93‑1 (MQ=255)

                                                                       ccTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCg                  >  2:2982188/1‑93 (MQ=255)

                                                                       ccTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCg                  >  2:1705949/1‑93 (MQ=255)

                                                                        cTTGCTCTACGGCCAGGCGGCTGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCgg                 <  2:3960928/93‑1 (MQ=255)

                                                                         ttGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGa                >  1:1721514/1‑93 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              >  1:1476435/1‑93 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              >  1:2197518/1‑93 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            >  1:10126/1‑93 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            >  1:442716/1‑93 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGGATTAGCGGAAATa            >  2:2053849/1‑93 (MQ=255)

                                                                              ctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATaa           >  2:1147999/1‑93 (MQ=255)

                                                                               tACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGACATAAg          >  2:469352/1‑93 (MQ=255)

                                                                                ccGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGc         <  1:4537432/92‑1 (MQ=255)

                                                                                 cGGCCGGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  1:1937465/93‑1 (MQ=255)

                                                                                 cGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  1:1297727/93‑1 (MQ=255)

                                                                                 cGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  1:4135278/93‑1 (MQ=255)

                                                                                   gCCAGGCGGCCTTTTTAAAAAAGCTTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa      <  2:2991560/93‑1 (MQ=255)

                                                                                   gCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa      <  2:3122355/93‑1 (MQ=255)

                                                                                   gCCAGGCGCCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGa      <  2:3327341/93‑1 (MQ=255)

                                                                                    ccAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc     >  1:3281652/1‑93 (MQ=255)

                                                                                    ccAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGAc     >  2:1825257/1‑93 (MQ=255)

                                                                                     cAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg    >  1:1051530/1‑93 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTgaagaa  >  1:2478749/1‑93 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGaa  <  1:4638873/93‑1 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGaa  <  1:1326772/93‑1 (MQ=255)

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`