BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,626,539	A→C	5.6%	T210P (ACG→CCG)	radA →	DNA repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,539	0	A	C	5.6%	313.5 / 4.5	119	T210P (ACG→CCG)	radA	DNA repair protein
Reads supporting (aligned to +/- strand): ref base A (50/60); new base C (3/5); total (53/66)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

GATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAAC  >  NC_000913/4626394‑4626688

gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAAcgcg                                                                                                                                                   >  1:623116/1‑150 (MQ=255)

  tCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAAcgcgcg                                                                                                                                                 <  2:741687/150‑1 (MQ=255)

       cTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGTCATCGCCGGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGtgtg                                                                                                                                            <  1:136115/150‑1 (MQ=255)

         gATTCCCGAAGAAGAGCAACCGAAGCTGATGGTACTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGc                                                                                                                                          <  2:1611457/150‑1 (MQ=255)

             gCCGAAGTAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCTGAGGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAtt                                                                                                                                      <  2:1003391/150‑1 (MQ=255)

             gCCGAAGAAGAGCAACGGATGCTGAGGGTAATTGACTCGATCCAGGTGATGCATATGGCGGCGGGACAGTCATCGCCTGGGAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAtt                                                                                                                                      <  1:85892/150‑1 (MQ=255)

             gCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAtt                                                                                                                                      >  2:1062174/1‑150 (MQ=255)

               cgacgaAGAGCAACCGAAGCTTATGGTAATTGACTCGATCCACGTGATGCATCTGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGt                                                                                                                                    <  2:1272338/150‑1 (MQ=255)

               cGAAGAAGAGCAACCGAAGCTGATGGTATTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCACACCGCGCGTTTTGGCGaagcg                                                                                                                                    <  1:105701/150‑5 (MQ=255)

               cGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCACACCGCGCGTTTTGGCGaagcg                                                                                                                                    <  1:332481/150‑5 (MQ=255)

               cGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGt                                                                                                                                    >  1:949588/1‑150 (MQ=255)

               cGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAATGt                                                                                                                                    >  1:901211/1‑150 (MQ=255)

               cGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATACAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACCCGCTTCGCCAAACCGCGCGGTGTGGCGAgtgt                                                                                                                                    >  2:98725/1‑150 (MQ=255)

                gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                                                                                   >  1:1491214/1‑150 (MQ=255)

                 aagaagAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCa                                                                                                                                  <  1:1603529/150‑1 (MQ=255)

                  agaagaGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGAGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAt                                                                                                                                 >  1:638549/1‑150 (MQ=255)

                   gaagaGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtg                                                                                                                                >  2:647730/1‑150 (MQ=255)

                    aagaGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGGGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtgg                                                                                                                               <  1:1783241/150‑1 (MQ=255)

                     agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggt                                                                                                                              >  2:1121710/1‑150 (MQ=255)

                     agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGACTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggt                                                                                                                              >  1:135404/1‑150 (MQ=255)

                      gagCCACAGATGCTGCTGGTAATTGACGCGATCCAGGTGCTGCATATGGCGGAGTTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGAGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggtg                                                                                                                             <  1:1435474/150‑1 (MQ=255)

                      gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggtg                                                                                                                             >  1:510333/1‑150 (MQ=255)

                       agCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggtgg                                                                                                                            >  2:839293/1‑150 (MQ=255)

                        gCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTggg                                                                                                                           <  1:1614843/150‑1 (MQ=255)

                          aaCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGc                                                                                                                         <  2:583395/150‑1 (MQ=255)

                             cGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACATTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACg                                                                                                                      <  2:851050/150‑1 (MQ=255)

                             cGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGCGGCACg                                                                                                                      >  2:1502913/1‑150 (MQ=255)

                                 gCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAAc                                                                                                                  >  1:773149/1‑150 (MQ=255)

                                 gCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAGGGTGGGGCACGTAAc                                                                                                                  >  1:142662/1‑150 (MQ=255)

                                    gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCaa                                                                                                               >  1:956452/1‑150 (MQ=255)

                                    gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCaa                                                                                                               >  2:1160954/1‑150 (MQ=255)

                                     aTGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCaaa                                                                                                              <  1:1810621/150‑1 (MQ=255)

                                      tGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAg                                                                                                             <  1:807530/150‑1 (MQ=255)

                                       ggTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGGCATCGCCTGGTAGCGTGGCGCAGGTGCGGGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGa                                                                                                            <  2:804119/150‑1 (MQ=255)

                                          aaTTGACTCGCTCCAGGTGATGCATATCGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATgg                                                                                                         <  2:149505/150‑1 (MQ=255)

                                               aCTCGCTCCAGGTCACCCCTATGGCGGAGGTACAGTCATCGCCTGGTAGCGTGGCGCACGCTCGTGAACCGGCGGCCTATTTGACACTCGGCGCCAAACCGCGCGGTGTGGCGTTTGTCAGGGTGGGGCACGTAACCAAAGCTGGtgccc                                                                                                    <  2:1818486/150‑5 (MQ=255)

                                               aCTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCgc                                                                                                    >  1:1805728/1‑150 (MQ=255)

                                               aCTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCgc                                                                                                    >  1:1798476/1‑150 (MQ=255)

                                                  cGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCgctgg                                                                                                 >  1:1008927/1‑150 (MQ=255)

                                                   gATCCAGGTGATGCATATGGCGGTTGTACAGTCATCTCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGTTGTGGCGATTGTCATGGTGGGGCACGTTACCAAAGATGGTTCgctggc                                                                                                <  2:1541051/150‑1 (MQ=255)

                                                   gATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCgctggc                                                                                                <  2:1808095/150‑1 (MQ=255)

                                                     tCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCgctggctg                                                                                              <  2:407898/150‑1 (MQ=255)

                                                     tCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGGGGCGCAGGGGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTGCgctggctg                                                                                              <  2:961000/150‑1 (MQ=255)

                                                       gaGGTGATGCCTTTGGCGGTTGTACTGTCCTGGCCTGGTAGCGGGGCGCTGGGGCGGGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCtggctggc                                                                                            <  1:1632387/149‑1 (MQ=255)

                                                        aGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGcc                                                                                           <  1:396316/150‑1 (MQ=255)

                                                         ggTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGccc                                                                                          <  2:537296/150‑1 (MQ=255)

                                                           ggATGCATATGGCGGATGTACAGTCATCGCCTGGTAGTGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGa                                                                                        <  1:630626/149‑1 (MQ=255)

                                                            gATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGaa                                                                                       <  1:1257915/150‑1 (MQ=255)

                                                               gCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGt                                                                                    >  2:1848960/1‑150 (MQ=255)

                                                                 atatGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGc                                                                                  <  2:600209/150‑1 (MQ=255)

                                                                    tGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTgg                                                                               <  2:135404/150‑1 (MQ=255)

                                                                     ggCGGATGTCCAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTTTTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGa                                                                              >  1:1499829/1‑150 (MQ=255)

                                                                      gCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGaa                                                                             <  2:1679376/150‑1 (MQ=255)

                                                                      gCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGaa                                                                             <  1:827376/150‑1 (MQ=255)

                                                                           tGTACAGTCCTCGCCTGGTAGCGTGGCGCAGGTCCGGGACACGGCGGCTTATTTGACACGCTTCGCCAAACCGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCACAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTg                                                                        <  2:264535/150‑1 (MQ=255)

                                                                             tACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTa                                                                      <  2:1502690/150‑1 (MQ=255)

                                                                             tACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGAGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTa                                                                      >  1:1465056/1‑150 (MQ=255)

                                                                             tACAGGCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTa                                                                      <  2:543553/150‑1 (MQ=255)

                                                                              aCAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACTCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTAt                                                                     >  2:554046/1‑150 (MQ=255)

                                                                               cAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATc                                                                    >  2:22254/1‑150 (MQ=255)

                                                                               cAGGCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATc                                                                    <  2:1302686/150‑1 (MQ=255)

                                                                                   cATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACt                                                                <  2:519827/150‑1 (MQ=255)

                                                                                   cATCGCCTGGGAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGGGGGGCACGTAACCACAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACt                                                                <  2:120421/150‑1 (MQ=255)

                                                                                    aTCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTg                                                               <  1:965394/150‑1 (MQ=255)

                                                                                      cGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGtt                                                             >  1:1608004/1‑150 (MQ=255)

                                                                                           ggTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGCACACTGTATCGACTGTTCGGTg                                                        >  2:1779039/1‑150 (MQ=255)

                                                                                           ggTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                        >  2:1094568/1‑150 (MQ=255)

                                                                                             tAGCGTGGCGCAGGTGCGTGAATCGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCt                                                      <  2:1267482/150‑1 (MQ=255)

                                                                                              aGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCAGGCTGGCCCGAAAGTGCTGGAACACTGAATCGACTGTTCGGTGCtt                                                     >  1:526589/1‑150 (MQ=255)

                                                                                              aGCGTGGCGCAGGTGCGTGAAACGGCGGCGTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCtt                                                     <  2:1684880/150‑1 (MQ=255)

                                                                                                  tGGCTCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCTCCAAAATGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACAAATGATGGTTCGGTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTgg                                                 <  1:13199/150‑1 (MQ=255)

                                                                                                   ggCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGa                                                <  2:1331629/150‑1 (MQ=255)

                                                                                                       cAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGc                                            <  2:623116/150‑1 (MQ=255)

                                                                                                         ggTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGa                                          <  2:773149/150‑1 (MQ=255)

                                                                                                                gAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGAc                                   <  1:948743/150‑1 (MQ=255)

                                                                                                                gAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGAGATGCCGAc                                   >  1:1096774/1‑150 (MQ=255)

                                                                                                                   aCGGCGGCTTAGTGGACACGCTTCGCCAAAACGCGCGGCGGGGCGATTGGCATGGGGGGGCACGTAACCAAAGATGGTTCGCGGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTcc                                <  1:409063/150‑1 (MQ=255)

                                                                                                                    cggcggCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATTGCGACGCCGACTccc                               >  2:1148965/1‑150 (MQ=255)

                                                                                                                     ggcggcTTATTTGACACGCTTCGCCAAAACGCGCGGTGGGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCg                              <  1:1515826/150‑1 (MQ=255)

                                                                                                                          cTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                         <  2:949588/150‑1 (MQ=255)

                                                                                                                          cTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                         <  1:1062174/150‑1 (MQ=255)

                                                                                                                          cTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                         <  2:1258104/150‑1 (MQ=255)

                                                                                                                          attttttGACAGGCTTCGCCAAATCGCGCGGTGTGGCGTTTGTCATTGTGTGGCACGTAACCAAAGATGGGTCGCTGGCTGGCCCGAAAGTGCTGGAACAATGTATCGAGTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                         <  1:142108/146‑1 (MQ=255)

                                                                                                                           ttATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                        <  2:1636202/150‑1 (MQ=255)

                                                                                                                           ttATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                        >  2:578671/1‑150 (MQ=255)

                                                                                                                            tATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGc                       >  2:1098162/1‑150 (MQ=255)

                                                                                                                             aTTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGTTGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGcc                      >  2:843566/1‑149 (MQ=255)

                                                                                                                                tGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGACGACTCCCGTTTTCGCACCt                   >  2:389514/1‑150 (MQ=255)

                                                                                                                                   cacGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGATCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCAGTGTTCGCACCTTgc                >  1:1648692/1‑150 (MQ=255)

                                                                                                                                   cacGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAAGGTGGGGAACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTgc                >  2:1008268/1‑150 (MQ=255)

                                                                                                                                     cGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTgcgc              >  1:1170486/1‑150 (MQ=255)

                                                                                                                                     acaatCGCCACACCGCGCGTTTTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTgcgc              >  2:834229/5‑150 (MQ=255)

                                                                                                                                      gCTTCGCCAAAACGCGCGGTGTGGCGATTGTGATGGGGGGGCCCGTCCCCACAGATGGTTGGCAGGCACGCCCGAAAGAGCTTGACCAATGTATCGACTGTTAGGTGCGTGTGGGTGGCGGTGCGGAGTCCCGTTTGGGGACCGTgcgcg             >  2:84135/1‑149 (MQ=255)

                                                                                                                                      gCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCa             <  2:1459976/150‑1 (MQ=255)

                                                                                                                                           gCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCAt        >  2:235761/1‑150 (MQ=255)

                                                                                                                                             cAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATaa      >  2:907849/1‑150 (MQ=255)

                                                                                                                                              aaaaCGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATaat     <  1:1082205/150‑2 (MQ=255)

                                                                                                                                               acaCGCGCGGGGTGGCGATTGTCATGGTGGGGCACTTAACCAAAGATGGTGCGCTGGCTGGCCCGAAAGTGCTGGAACACTGAATCGACTGTTCGTTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATaaaa    <  1:1290778/148‑1 (MQ=255)

                                                                                                                                               aaaCGCGCGGGGTGGCGATTGTCATGGTGGTGCACGCAACCAAAGATGTTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATaaaa    <  2:1221791/150‑1 (MQ=255)

                                                                                                                                                 aCGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACATTGCGCAGCCATAAAAAc  >  1:998988/1‑150 (MQ=255)

GATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAAC  >  NC_000913/4626394‑4626688

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