BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,825	G→A	P163P (CCG→CCA)	b2881 →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,825	0	G	A	99.8%	1520.4 / NA	420	P163P (CCG→CCA)	b2881	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (193/226); total (194/226)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 439 total aligned reads displayed.

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

    gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTc                                                                                     <  1:4703731/93‑1 (MQ=255)

    gATGCTGGCTGGCAGCAATACTATCAGGTCACTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCCCCATTGATATCGCTCCAACTTTc                                                                                     <  1:7619536/93‑1 (MQ=255)

      tGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCg                                                                                   <  1:5579387/93‑1 (MQ=255)

       gctggcttGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  1:8574526/93‑1 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  2:2570930/93‑1 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  2:8142104/93‑1 (MQ=255)

         tgcctgccAGCAATACTATCCGGTCATTGAACTGGCGGTTGCACGTAAAAGTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                <  2:68031/90‑1 (MQ=255)

            ctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacga                                                                             <  1:6788163/93‑1 (MQ=255)

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacgac                                                                            >  2:4100716/1‑93 (MQ=255)

              ggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAcc                                                                           <  1:5164290/93‑1 (MQ=255)

                 agcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAg                                                                        >  1:9049439/1‑93 (MQ=255)

                  gcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTaca                                                                       >  2:2283454/1‑91 (MQ=255)

                  gcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGa                                                                       <  2:6544259/93‑1 (MQ=255)

                  gcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGa                                                                       <  2:1875051/93‑1 (MQ=255)

                    aaTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAg                                                                     <  2:7873060/93‑1 (MQ=255)

                     aTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGt                                                                    >  1:4186686/1‑93 (MQ=255)

                     aTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGt                                                                    <  1:1808685/93‑1 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   <  1:9766842/93‑1 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   >  1:270277/1‑93 (MQ=255)

                       aCTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCa                                                                  >  2:8946997/1‑93 (MQ=255)

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 <  1:4831441/93‑1 (MQ=255)

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 >  2:1211382/1‑93 (MQ=255)

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 >  2:6053365/1‑93 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                >  2:953600/1‑93 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  1:2249334/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                >  1:7203286/1‑93 (MQ=255)

                          aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTg                                                               <  2:5791493/93‑1 (MQ=255)

                           tCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTgg                                                              >  1:1965326/1‑93 (MQ=255)

                             aGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTa                                                            <  1:578423/93‑1 (MQ=255)

                             aGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACCTTCCGTGACGACCTAGTAGTCATTGGTa                                                            >  1:5091084/1‑93 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  2:5554374/93‑1 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  2:4704520/93‑1 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  1:7067757/93‑1 (MQ=255)

                               gTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAg                                                          <  1:2879799/93‑1 (MQ=255)

                                 cATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCa                                                        >  1:6260259/1‑93 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       >  1:9292959/1‑93 (MQ=255)

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      >  1:8876758/1‑93 (MQ=255)

                                    tGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTa                                                     <  2:4423650/93‑1 (MQ=255)

                                    tGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTa                                                     >  2:5287513/1‑93 (MQ=255)

                                     gAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  2:7981278/93‑1 (MQ=255)

                                      aaCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATc                                                   >  2:3992154/1‑93 (MQ=255)

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  >  2:2008861/1‑93 (MQ=255)

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  <  2:2862739/93‑1 (MQ=255)

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  >  1:4329294/1‑93 (MQ=255)

                                        cTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCt                                                 >  2:7382081/1‑93 (MQ=255)

                                         tGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATACTa                                                >  2:8223731/1‑93 (MQ=255)

                                          ggCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTaa                                               <  1:3423553/93‑1 (MQ=255)

                                           gCGGGTGCACGTAAAAATAATACGCAGGCCAGCATTGATATAGCTCCAAATTTCCGTGACGACCTAGCAGACATTGGTACGCATTATCCTaaa                                              >  2:8386146/1‑93 (MQ=255)

                                            cGGTTGCACGTAAAAATAATCGGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTACCCTaaaa                                             <  1:5854216/93‑1 (MQ=255)

                                             ggTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAAc                                            <  2:340779/93‑1 (MQ=255)

                                              gTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACt                                           >  2:5024450/1‑93 (MQ=255)

                                               ttGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTg                                          <  2:4530278/93‑1 (MQ=255)

                                                tGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGa                                         <  2:3759863/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        <  2:6600307/93‑1 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTGGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  2:3887954/93‑1 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  1:2501055/93‑1 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  2:4056165/93‑1 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  2:4783551/93‑1 (MQ=255)

                                                    cGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGcc                                     >  2:562997/1‑93 (MQ=255)

                                                    cGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGcc                                     >  2:9107659/1‑93 (MQ=255)

                                                     gTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcg                                    >  1:3924118/1‑93 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCCCCATTGATATCGCCTCAACCTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   <  2:2803297/93‑1 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   >  1:7401215/1‑93 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTGTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   <  1:1393814/93‑1 (MQ=255)

                                                        aaaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGa                                 >  2:2973053/1‑93 (MQ=255)

                                                         aaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaa                                <  1:9870100/93‑1 (MQ=255)

                                                         aaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaa                                <  1:6436586/93‑1 (MQ=255)

                                                          aataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaa                               >  2:752764/1‑93 (MQ=255)

                                                          aataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaa                               >  2:176482/1‑93 (MQ=255)

                                                           ataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaaa                              <  2:6850796/93‑1 (MQ=255)

                                                            taatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             <  2:5420484/93‑1 (MQ=255)

                                                            taatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTAACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAt                             <  2:8749619/93‑1 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:1281517/93‑1 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:315555/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  1:748480/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  1:7986121/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           >  1:8789579/1‑93 (MQ=255)

                                                               tCCTCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATggg                          >  1:3636423/1‑92 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATggg                          >  2:2181587/1‑92 (MQ=255)

                                                                ccGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTg                         <  1:1066933/93‑1 (MQ=255)

                                                                ccGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTg                         >  2:2393484/1‑93 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTAATGCCGCGAAAATGGTGCa                       <  2:7450334/93‑1 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     <  2:6248661/93‑1 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     >  2:7661613/1‑93 (MQ=255)

                                                                        caccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaa                 >  1:4089420/1‑93 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                >  2:6456275/1‑93 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                <  2:93094/93‑1 (MQ=255)

                                                                          ccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAc               >  1:1221306/1‑93 (MQ=255)

                                                                          ccaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAc               <  1:8180243/93‑1 (MQ=255)

                                                                             ttGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCt            >  2:7310081/1‑93 (MQ=255)

                                                                             ttGATATCGCTCCAACGTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAACCTAAACCGGATGCCGCGAAAGTGGTGCAGGGGAAACCCt            <  1:5285787/93‑1 (MQ=255)

                                                                               gATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGc          <  2:1394201/93‑1 (MQ=255)

                                                                               gATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGc          <  1:9511483/93‑1 (MQ=255)

                                                                                atatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCt         <  1:211593/93‑1 (MQ=255)

                                                                                 tatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTa        >  1:2950009/1‑93 (MQ=255)

                                                                                  atCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTAt       >  2:7016940/1‑93 (MQ=255)

                                                                                     gCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGtt    <  1:5327627/93‑1 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  2:3619811/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  2:8763902/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  2:9800/1‑93 (MQ=255)

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

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