BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,355	A→G	H314R (CAC→CGC)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,355	0	A	G	98.0%	467.2 / ‑2.1	149	H314R (CAC→CGC)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (80/66); total (82/67)
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.96e-01

Only 100 of 161 total aligned reads displayed.

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCcgcgt                                                                                        >  2:4376581/1‑92 (MQ=255)

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        >  1:909313/1‑93 (MQ=255)

      gcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGatgatg                                                                                   >  1:4067319/1‑93 (MQ=255)

       cTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGt                                                                                  <  1:4183795/93‑1 (MQ=255)

       cTGTAGCTACACGCCGGATCAAGGTTACACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGt                                                                                  >  2:3781389/1‑93 (MQ=255)

          gaGCTACACGCCGGATCAAGGTTCCACGCGCATTTAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCGCCAGCCACGATGATGTCTg                                                                               <  2:1563067/92‑1 (MQ=255)

           aGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGt                                                                              >  1:4675091/1‑93 (MQ=255)

            gCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGtt                                                                             <  1:3270273/93‑1 (MQ=255)

             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            <  2:2944772/93‑1 (MQ=255)

               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttga                                                                          >  1:4163805/1‑93 (MQ=255)

                cacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgat                                                                         >  2:2689018/1‑93 (MQ=255)

                cacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgat                                                                         <  2:3666208/93‑1 (MQ=255)

                 acGCGGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        <  1:282790/93‑1 (MQ=255)

                 acGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        >  2:1098883/1‑93 (MQ=255)

                 acGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        >  2:1098884/1‑93 (MQ=255)

                 acGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        <  1:2861492/93‑1 (MQ=255)

                 acGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        >  1:48205/1‑93 (MQ=255)

                    ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     >  2:3061312/1‑93 (MQ=255)

                    ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     >  2:3728128/1‑93 (MQ=255)

                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   >  1:4294862/1‑93 (MQ=255)

                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   <  2:1917784/93‑1 (MQ=255)

                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   >  2:3489471/1‑93 (MQ=255)

                       gATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGttt                                                                  <  2:2786104/93‑1 (MQ=255)

                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 >  2:4017437/1‑93 (MQ=255)

                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 <  2:4326852/93‑1 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACTAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                >  2:258372/1‑93 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGCTGTCTGTTTGATTGAGttttc                                                                >  2:3254193/1‑91 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                >  1:1485224/1‑93 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGGGCTGGCATCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                <  1:2560088/93‑1 (MQ=255)

                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               >  1:2041988/1‑93 (MQ=255)

                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               <  1:915073/93‑1 (MQ=255)

                           aaGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAgg                                                              <  2:2438568/93‑1 (MQ=255)

                           aaGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAgg                                                              >  1:3887796/1‑93 (MQ=255)

                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  1:2186674/93‑1 (MQ=255)

                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  1:1252367/93‑1 (MQ=255)

                             ggTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTg                                                            >  1:1055218/1‑93 (MQ=255)

                               ttCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGcc                                                          <  2:87991/93‑1 (MQ=255)

                               ttCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGcc                                                          <  1:3731940/93‑1 (MQ=255)

                                tCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGccc                                                         <  2:3581741/93‑1 (MQ=255)

                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      >  1:4279068/1‑93 (MQ=255)

                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     <  2:4105764/93‑1 (MQ=255)

                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     <  1:2239417/93‑1 (MQ=255)

                                     gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAg                                                    >  1:2897060/1‑93 (MQ=255)

                                       gcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTc                                                  <  1:4551783/93‑1 (MQ=255)

                                        cATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCa                                                 <  1:3238602/93‑1 (MQ=255)

                                         cttGACCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGTTGATGTCTTTTTGATTGATTTTCAGGTGCCCGCCAGTCAg                                                <  2:2536875/92‑1 (MQ=255)

                                         aTTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAg                                                >  2:1391110/1‑93 (MQ=255)

                                          ttGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               >  2:3510549/1‑93 (MQ=255)

                                          ttGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               >  2:4114109/1‑93 (MQ=255)

                                          ttGAACGCGTGCTGGCCTCCGGAACTGGTGCGCTTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               <  1:4211566/93‑1 (MQ=255)

                                            gAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAt                                             <  2:1795617/93‑1 (MQ=255)

                                              aCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAttt                                           >  1:1801643/1‑93 (MQ=255)

                                                gcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCa                                         <  1:1912149/93‑1 (MQ=255)

                                                gcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCa                                         >  1:3077911/1‑93 (MQ=255)

                                                 cgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCaa                                        <  2:2621610/93‑1 (MQ=255)

                                                  gTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCaaa                                       >  2:2963344/1‑93 (MQ=255)

                                                   tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      >  2:3853860/1‑93 (MQ=255)

                                                    gCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACt                                     >  2:1452373/1‑93 (MQ=255)

                                                       ggCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  >  1:2025845/1‑93 (MQ=255)

                                                       ggCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  >  2:550402/1‑93 (MQ=255)

                                                        gCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcg                                 <  1:4138414/93‑1 (MQ=255)

                                                          cTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCa                               >  1:2928550/1‑93 (MQ=255)

                                                          cTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCa                               <  1:728122/93‑1 (MQ=255)

                                                          cTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCa                               <  1:4708474/93‑1 (MQ=255)

                                                             cGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaa                            >  1:2958280/1‑93 (MQ=255)

                                                              ggTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaaa                           <  1:616657/93‑1 (MQ=255)

                                                               gTATTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAag                          <  2:487762/93‑1 (MQ=255)

                                                               gTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAag                          >  1:700377/1‑93 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCCCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTCCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         <  1:488005/93‑1 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         >  1:147072/1‑93 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         <  1:4341139/93‑1 (MQ=255)

                                                                gaCTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         <  1:877466/92‑1 (MQ=255)

                                                                 aCTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAagag                        <  2:4575617/93‑1 (MQ=255)

                                                                  cTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAgaga                       <  2:256238/93‑1 (MQ=255)

                                                                  cTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAgaga                       <  2:2385279/93‑1 (MQ=255)

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      <  2:4445778/93‑1 (MQ=255)

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      <  2:3067008/93‑1 (MQ=255)

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      >  1:1797570/1‑93 (MQ=255)

                                                                    ggTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATc                     >  1:1335346/1‑93 (MQ=255)

                                                                     gTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCg                    <  2:703483/93‑1 (MQ=255)

                                                                      tgtgcgTATTGTCACCAGCCGCGATTATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   <  2:1328796/90‑1 (MQ=255)

                                                                         gcgTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCa                <  2:1492822/93‑1 (MQ=255)

                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             >  1:905864/1‑93 (MQ=255)

                                                                             aTTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATc            >  2:1393934/1‑93 (MQ=255)

                                                                             aTTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATc            >  1:3871809/1‑93 (MQ=255)

                                                                             aTTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATc            <  1:25775/93‑1 (MQ=255)

                                                                              ttGTCCCCAGCCGCGATGATGTCTGTTTGATTGATTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           <  2:2536776/93‑1 (MQ=255)

                                                                              ttGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           >  2:4626401/1‑93 (MQ=255)

                                                                               tGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCt          <  2:840296/93‑1 (MQ=255)

                                                                                gTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTg         >  2:2921954/1‑93 (MQ=255)

                                                                                 tCACCAGCCGCGATGATGTCTGATTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGCCCAAATCCTGa        >  1:3214584/1‑93 (MQ=255)

                                                                                  cACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATGCTGaa       >  1:261469/1‑93 (MQ=255)

                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     >  1:3453659/1‑93 (MQ=255)

                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     <  1:1240104/93‑1 (MQ=255)

                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     >  1:3768725/1‑93 (MQ=255)

                                                                                      aGCCGCGATGATGTCTTTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   <  1:3346367/93‑1 (MQ=255)

                                                                                      aGCCGCGATGATGTCTGTTTGATTGAGTTTCAGTGGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   <  1:682843/93‑1 (MQ=255)

                                                                                      aGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   <  2:499020/93‑1 (MQ=255)

                                                                                      aGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   >  1:1052133/1‑93 (MQ=255)

                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg  >  2:17099/1‑93 (MQ=255)

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

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